G. Gasparin

Putative Regulatory Factors Associated with Intramuscular Fat Content

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.

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

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

SNP and INDEL detection in a QTL region on chicken chromosome 2 associated with muscle deposition

Genetic improvement is important for the poultry industry, contributing to increased efficiency of meat production and quality. Because breast muscle is the most valuable part of the chicken carcass, knowledge of polymorphisms influencing this trait can help breeding programs. Therefore, the complete genome of 18 chickens from two different experimental lines (broiler and layer) from EMBRAPA was sequenced, and SNPs and INDELs were detected in a QTL region for breast muscle deposition on chicken chromosome 2 between microsatellite markers MCW0185 and MCW0264 (105,849-112,649 kb). Initially, 94,674 unique SNPs and 10,448 unique INDELs were identified in the target region. After quality filtration, 77% of the SNPs (85,765) and 60% of the INDELs (7828) were retained. The studied region contains 66 genes, and functional annotation of the filtered variants identified 517 SNPs and three INDELs in exonic regions. Of these, 357 SNPs were classified as synonymous, 153 as non-synonymous, three as stopgain, four INDELs as frameshift and three INDELs as non-frameshift. These exonic mutations were identified in 37 of the 66 genes from the target region, three of which are related to muscle development (DTNA, RB1CC1 and MOS). Fifteen non-tolerated SNPs were detected in several genes (MEP1B, PRKDC, NSMAF, TRAPPC8, SDR16C5, CHD7, ST18 and RB1CC1). These loss-of-function and exonic variants present in genes related to muscle development can be considered candidate variants for further studies in chickens. Further association studies should be performed with these candidate mutations as should validation in commercial populations to allow a better explanation of QTL effects.

Variant discovery in a QTL region on chromosome 3 associated with fatness in chickens

Abdominal fat content is an economically important trait in commercially bred chickens. Although many quantitative trait loci (QTL) related to fat deposition have been detected, the resolution for these regions is low and functional variants are still unknown. The current study was conducted aiming at increasing resolution for a region previously shown to have a QTL associated with fat deposition, to detect novel variants from this region and to annotate those variants to delineate potentially functional ones as candidates for future studies. To achieve this, 18 chickens from a parental generation used in a reciprocal cross between broiler and layer lines were sequenced using the Illumina next-generation platform with an initial coverage of 18X/chicken. The discovery of genetic variants was performed in a QTL region located on chromosome 3 between microsatellite markers LEI0161 and ADL0371 (33,595,706-42,632,651 bp). A total of 136,054 unique SNPs and 15,496 unique INDELs were detected in this region, and after quality filtering, 123,985 SNPs and 11,298 INDELs were retained. Of these variants, 386 SNPs and 15 INDELs were located in coding regions of genes related to important metabolic pathways. Loss-of-function variants were identified in several genes, and six of those, namely LOC771163, EGLN1, GNPAT, FAM120B, THBS2 and GGPS1, were related to fat deposition. Therefore, these loss-of-function variants are candidate mutations for conducting further studies on this important trait in chickens.

Identification of KCNJ11 as a functional candidate gene for bovine meat tenderness

The potassium inwardly rectifying channel, subfamily J, member 11 (KCNJ11) gene was investigated as a candidate for meat tenderness based on the effects reported on muscle for KCNJ11 gene knockout in rat models and its position in a quantitative trait locus (QTL) for meat tenderness in the bovine genome. Sequence variations in the KCNJ11 gene were described by sequencing six amplified fragments, covering almost the entire gene. We identified single nucleotide polymorphisms (SNP) and validated them by different approaches, taking advantage of simultaneous projects that are being developed with the same Nelore population. By sequencing the KCNJ11 in Nelore steers representing extreme phenotypes for Warner-Bratzler shear force (WBSF), it was possible to identify 22 SNPs. We validated two of the identified markers by genotyping the whole population (n = 460). Analysis of association between genotypes and WBSF values revealed a significant additive effect of a SNP at different meat aging times (P ≤ 0.05). In addition, an association between the expression levels of KCNJ11 and WBSF was found, with lower expression levels of KCNJ11 associated with more tender meat (P ≤ 0.05). The results showed that the KCNJ11 gene is a candidate mapped to a QTL for meat tenderness previously identified on BTA15 and may be useful to identify animals with genetic potential to produce tender meat. The effect of KCNJ11 observed on muscle is potentially due to changes in activity of KATP channels, which in turn influence the flow of potassium in the intracellular space, allowing establishment of the membrane potential necessary for muscle contraction.

Genome-wide characterization of genetic variants and putative regions under selection in meat and egg-type chicken lines

BackgroundMeat and egg-type chickens have been selected for several generations for different traits. Artificial and natural selection for different phenotypes can change frequency of genetic variants, leaving particular genomic footprints throghtout the genome. Thus, the aims of this study were to sequence 28 chickens from two Brazilian lines (meat and white egg-type) and use this information to characterize genome-wide genetic variations, identify putative regions under selection using Fst method, and find putative pathways under selection.ResultsA total of 13.93 million SNPs and 1.36 million INDELs were identified, with more variants detected from the broiler (meat-type) line. Although most were located in non-coding regions, we identified 7255 intolerant non-synonymous SNPs, 512 stopgain/loss SNPs, 1381 frameshift and 1094 non-frameshift INDELs that may alter protein functions. Genes harboring intolerant non-synonymous SNPs affected metabolic pathways related mainly to reproduction and endocrine systems in the white-egg layer line, and lipid metabolism and metabolic diseases in the broiler line. Fst analysis in sliding windows, using SNPs and INDELs separately, identified over 300 putative regions of selection overlapping with more than 250 genes. For the first time in chicken, INDEL variants were considered for selection signature analysis, showing high level of correlation in results between SNP and INDEL data. The putative regions of selection signatures revealed interesting candidate genes and pathways related to important phenotypic traits in chicken, such as lipid metabolism, growth, reproduction, and cardiac development.ConclusionsIn this study, Fst method was applied to identify high confidence putative regions under selection, providing novel insights into selection footprints that can help elucidate the functional mechanisms underlying different phenotypic traits relevant to meat and egg-type chicken lines. In addition, we generated a large catalog of line-specific and common genetic variants from a Brazilian broiler and a white egg layer line that can be used for genomic studies involving association analysis with phenotypes of economic interest to the poultry industry.

Gene expression profiling by high throughput sequencing to determine signatures for the bovine receptive uterus at early gestation

The uterus plays a central role among the reproductive tissues in the context of early embryo-maternal communication and a successful pregnancy depends on a complex series of endometrial molecular and cellular events. The factors responsible for the initial interaction between maternal and embryonic tissues, leading to the establishment of pregnancy, remain poorly understood. In this context, Illuminas next-generation sequencing technology has been used to discover the uterine transcriptome signature that is favourable for ongoing pregnancy. More specifically, the present report documents on a retrospective in vivo study in which data on pregnancy outcome were linked to uterine gene expression signatures on day 6 (bovine model). Using the RNA-Seq method, 14.654 reference genes were effectively analysed for differential expression between pregnant and non-pregnant uterine tissue. Transcriptome data revealed that 216 genes were differently expressed when comparing uterine tissue from pregnant and non-pregnant cows. All read sequences were deposited in the Sequence Read Archive (SRA) of the NCBI (http://www.ncbi.nlm.nih.gov/sra). An overview of the gene expression data has been deposited in NCBIs Gene Expression Omnibus (GEO) and is accessible through GEO Series accession number GSE65117. This allows the research community to enhance reproducibility and allows for new discoveries by comparing datasets of signatures linked to receptivity and/or pregnancy success. The resulting information can serve as tool to identify valuable and urgently needed biomarkers for scoring maternal receptivity and even for accurate detection of early pregnancy, which is a matter of cross-species interest. Beyond gene expression analysis as a marker tool, the RNA-Seq information on pregnant uterine tissue can be used to gain novel mechanistic insights, such as by identifying alternative splicing events, allele-specific expression, and rare and novel transcripts that might be involved in the onset of maternal receptivity. This concept is unique and provides a new approach towards strategies that are highly needed to improve efficiency of fertility treatments.

Allele- and parent-of-origin-specific effects on expression of the KCNJ11 gene: A candidate for meat tenderness in cattle

In contrast to the Mendelian inheritance model, parental alleles can contribute unequally to gene expression, which may result in phenotypic variance among individuals and bias in the predicted additive effect of molecular markers associated with production traits. Given the need to understand the effects of allelic variation and parent-of-origin effects on the expression of genes with a commercial interest in cattle, we analyzed the expression of KCNJ11 (potassium inwardly rectifying channel, subfamily J, member 11), which was previously described as a functional candidate gene for meat tenderness. Allele-specific and parent-of-origin-dependent expression of this gene were assessed in bovine muscle using the rs379610823 single nucleotide polymorphism as a reference. Biallelic expression was observed; however, the T allele was expressed at significantly higher levels than the C allele. Furthermore, increased expression of KCNJ11 was found in animals harboring the maternal T allele. This study is the first to describe the differential allelic expression of bovine KCNJ11. Our findings are important for understanding the mechanisms that underlie the pattern of KCNJ11 expression and its potential impact on the phenotypic variation of meat tenderness in Nelore beef cattle. This reinforces the need for further investigation of allelic- and parent-of-origin expression deviation in genetic markers eligible for the selection of target traits.

Gene Co-expression Analysis Indicates Potential Pathways and Regulators of Beef Tenderness in Nellore Cattle

Beef tenderness, a complex trait affected by many factors, is economically important to beef quality, industry, and consumer’s palatability. In this study, RNA-Seq was used in network analysis to better understand the biological processes that lead to differences in beef tenderness. Skeletal muscle transcriptional profiles from 24 Nellore steers, selected by extreme estimated breeding values (EBVs) for shear force after 14 days of aging, were analyzed and 22 differentially expressed transcripts were identified. Among these were genes encoding ribosomal proteins, glutathione transporter ATP-binding cassette, sub-family C (CFTR/MRP), member 4 (ABCC4), and synaptotagmin IV (SYT4). Complementary co-expression analyses using Partial Correlation with Information Theory (PCIT), Phenotypic Impact Factor (PIF) and the Regulatory Impact Factor (RIF) methods identified candidate regulators and related pathways. The PCIT analysis identified ubiquitin specific peptidase 2 (USP2), growth factor receptor-bound protein 10 (GBR10), anoctamin 1 (ANO1), and transmembrane BAX inhibitor motif containing 4 (TMBIM4) as the most differentially hubbed (DH) transcripts. The transcripts that had a significant correlation with USP2, GBR10, ANO1, and TMBIM4 enriched for proteasome KEGG pathway. RIF analysis identified microRNAs as candidate regulators of variation in tenderness, including bta-mir-133a-2 and bta-mir-22. Both microRNAs have target genes present in the calcium signaling pathway and apoptosis. PIF analysis identified myoglobin (MB), enolase 3 (ENO3), and carbonic anhydrase 3 (CA3) as potentially having fundamental roles in tenderness. Pathways identified in our study impacted in beef tenderness included: calcium signaling, apoptosis, and proteolysis. These findings underscore some of the complex molecular mechanisms that control beef tenderness in Nellore cattle.