Kacper Żukowski

Comprehensive analysis of the whole transcriptomes from two different pig breeds using RNA-Seq method

Next-generation sequencing RNA-Seq technology is a powerful tool that creates new possibilities for whole-transcriptome analysis. In our study, the RNA-Seq method was applied to analyze global changes in transcriptome from muscle tissue (m. semimembranosus) in two pig breeds (Pietrain and Polish Landrace, PL). The breeds differ in terms of muscularity, growth rate and reproduction traits. Using three different approaches (deseq, cufflinks and edger) and taking into account the most restrictive criteria, 35 genes differentially expressed between Pietrain and PL pigs were identified. In both breeds, the most abundant were transcripts encoding ribosomal and cytoskeletal proteins (TPM3, TCAP, TMOD4, TPM2, TNNC1) and calcium-binding proteins involved in muscle contraction, calcium-mediated signaling or cation transport (CASQ1, MLC2V, SLC25A4, MYL3). In PL pigs, we identified up-regulation of several genes that play crucial roles in reproduction: female gamete generation (BDP1, PTPN21, USP9X), fertilization (EGFR) and embryonic development (CPEB4). In the Pietrain breed, only seven genes were over-expressed (CISH, SPP1, TUBA8, ATP6V1C2, IGKC, predicted LOC100510960 and LOC100626400), and they play important roles in, for example, negative regulation of apoptosis, immune response, cell-cell signaling, cell growth and migration as well as the metabolic process. The functions of the majority of selected genes were consistent with phenotypic variation in investigated breeds; thus, we proposed a new panel of candidate genes that can be associated with economically important pig traits.

Genome-wide RNA-Seq analysis of breast muscles of two broiler chicken groups differing in shear force

In this study, a whole transcriptome analysis of breast muscles was conducted in broiler chicken groups differing in shear force. Shear force is a determinant of tenderness, which in turn is one of the most important parameters of meat quality in chickens. In our analysis, a total of 11,560 transcripts and 9824 genes per sample were identified. In chickens with more tender meat, up-regulation of 19 genes and down-regulation of 49 genes was observed. The up-regulated gene group included the ASB2 gene, which is probably involved in the meat conversion process, as its product results in the degradation of filamins, proteins which form muscle fibres. In the down-regulated gene group, genes which play a role in lipogenesis (THRSP, PLIN1) and in collagen synthesis (P4HA3, LEPREL4, PCOLCE2, COL16A1, COL20A1, VWA1) were detected. Their presence suggests the involvement of the extracellular matrix in the determination of meat tenderness. Thus, our study identified a pool of genes that may participate in the tenderisation process in broiler chickens.

Transcriptome profiling of Arabian horse blood during training regimens

BackgroundArabian horses are believed to be one of the oldest and most influential horse breeds in the world. Blood is the main tissue involved in maintaining body homeostasis, and it is considered a marker of the processes taking place in the other tissues. Thus, the aim of our study was to identify the genetic basis of changes occurring in the blood of Arabian horses subjected to a training regimen and to compare the global gene expression profiles between different training periods (T1: after a slow canter phase that is considered a conditioning phase, T2: after an intense gallop phase, and T3: at the end of the racing season) and between trained and untrained horses (T0). RNA sequencing was performed on 37 samples with a 75-bp single-end run on a HiScanSQ platform (Illumina), and differentially expressed genes (DEGs) were identified based on DESeq2 (v1.11.25) software.ResultsAn increase in the number of DEGs between subsequent training periods was observed, and the highest amount of DEGs (440) was detected between untrained horses (T0) and horses at the end of the racing season (T3). The comparisons of the T2 vs. T3 transcriptomes and the T0 vs. T3 transcriptomes showed a significant gain of up-regulated genes during long-term exercise (up-regulation of 266 and 389 DEGs in the T3 period compared to T2 and T0, respectively). Forty differentially expressed genes were detected between the T1 and T2 periods, and 296 between T2 and T3. Functional annotation showed that the most abundant genes up-regulated in exercise were involved in pathways regulating cell cycle (PI3K-Akt signalling pathway), cell communication (cAMP-dependent pathway), proliferation, differentiation and apoptosis, as well as immunity processes (Jak-STAT signalling pathway).ConclusionsWe investigated whether training causes permanent transcriptome changes in horse blood as a reflection of adaptation to conditioning and the maintenance of fitness to compete in flat races. The present study identified the overrepresented molecular pathways and genes that are essential for maintaining body homeostasis during long-term exercise in Arabian horses. Selected DEGs should be further investigated as markers that are potentially associated with racing performance in Arabian horses.

Transcript variants of a region on SSC15 rich in QTLs associated with meat quality in pigs

Abstract A high meat percentage in the porcine carcass has been achieved as a result of selection, but it has contributed to a deterioration of pork quality. The level of intramuscular fat has significantly declined, the pork has lost its tenderness and drip loss in meat has substantially increased, which has led to a deterioration of meat flavour and its technological suitability. The recovery of good pork quality could be supported by the development of genetic markers enabling faster breeding progress. This study presents a method by using RNA-seq data that identifies new variants for a chromosome region rich in QTLs for pork quality and selects gene candidates for these traits. This work included two pig breeds: the Polish Landrace (PL) and Puławska (PUL), which differ in meat quality and fat content. The transcriptome profile was estimated for semimembranosus and longissimus dorsi muscles. Into variant calling analysis, transcripts of both muscles encoded by genes located in a region between microsatellites SW964 and SW906 (43-135.9 Mbp) in SSC15 were included. In total, 439 transcripts were searched, 2,800 gene variants were identified and 6 mutations with a high effect belonging to the frameshift variants were found (ENSSSCG00000015976, ENSSSCG00000027516, WRN and XIRP2). Moreover, several interesting significant missense variants in PDLIM3, PLCD4 and SARAF genes were detected. These genes are recommended as candidates for meat quality; however they require further investigation in an association study.

Effects of Different Sources of Fat in the Diet of Pigs on the Liver Transcriptome Estimated by RNA-Seq

Abstract In this study, we have attempted to analyse the impact of dietary fats on the liver transcriptome in pigs. Four nutritional groups were created. The animals’ diets differed among groups in terms of the presence of corn dried distillers’ grains with solubles (DDGS) (group I - no DDGS, groups II, III, IV - 20% DDGS) as well as the type of fat (rapeseed oil - groups I and II, beef tallow - group III, coconut oil - group IV) used. Using the RNA-Seq method we identified 39 differentially expressed genes (DEGs) as a result of Cuffdiff analysis of the differences among all groups. Analysis of these genes with Panther Gene Classification System revealed that among identified DEGs, genes responsible for lipid and fatty acids metabolism are overrepresented as well as the genes engaged in oxidoreductase and catalytic activity. The article presents for the first time the RNAseq analysis of the liver transcriptome in pigs fed with different sources of fats. The results may be useful for the elaboration of new therapies for cardiovascular diseases in humans as well as for the preparation of new nutrition strategies in animals.

Deep sequencing of a QTL-rich region spanning 128-136Mbp of pig chromosome 15

The present study shows the characterization of the chromosome 15 (SSC15) region that is highly rich in quantitative traits loci (QTLs) associated with pork quality, growth performance, fat and meat carcass contents. The analytic method that was utilized included targeted enrichment DNA sequencing and RNA hybridisation probes. The research included two pig breeds (Puławska and Polish Landrace) that are significantly different in terms of carcass and meat quality features. Filtered sequences were aligned to the Sscrofa10.2 assembly genome with the STAR aligner and GATK HaplotypeCaller was used for identified gene variants in SSC15 region. In Puławska pigs, which were characterized by high meat quality, mutations were predominantly observed in non-coding regions such as introns and intergenic regions. The highest over 50% frequencies of alternate alleles were identified in the introns of TNS1, VIL1 and USP37 genes. In the upstream gene regions of the Polish Landrace pigs, were observed more mutations than in the upstream gene regions of Puławska. The present study showed interesting gene variant panel that could be analyzed in the further association studies in order to understand the impact on important productive pig traits.

The effect of QTL-rich region polymorphisms identified by targeted DNA-seq on pig production traits.

The aim of the present study was to analyse the effect of PLCD4, PECR, FN1 and PNKD mutations on pig productive traits and tested the usefulness of targeted enrichment DNA sequencing method as tool for preselection of genetic markers. The potential genetic markers for pig productive traits were identified by using targeted enrichment DNA sequencing of chromosome 15 region that is QTL-rich. The selected mutations were genotyped by using HRM, Sanger sequencing and PCR-ACRS methods. The association study was performed by using GLM model in SAS program and included over 500 pigs of 5 populations maintained in Poland. The variation (C/T) of PLCD4 gene affected feed conversion, intramuscular fat and water exudation. The PNKD mutations were associated with texture parameters measured after cooking. In turn, the variation rs792423408 (C/T) in the FN1 gene influenced toughness measured in semimembranosus muscle and growth traits that was observed particularly in Duroc breed. Summarizing, the investigated gene variants delivered valuable information that could be used during developing SNP microarray for genomic estimated breeding value procedure in pigs. Moreover, the study showed that the TEDNA-seq method could be used to preselect the molecular markers associated with pig traits. However, the further association study that included large number animal populations is necessary.

Examining the Genetic Background of Porcine Muscle Growth and Development Based on Transcriptome and miRNAome Data

Recently, selection in pigs has been focused on improving the lean meat content in carcasses; this focus has been most evident in breeds constituting a paternal component in breeding. Such sire-breeds are used to improve the meat quantity of cross-breed pig lines. However, even in one breed, a significant variation in the meatiness level can be observed. In the present study, the comprehensive analysis of genes and microRNA expression profiles in porcine muscle tissue was applied to identify the genetic background of meat content. The comparison was performed between whole gene expression and miRNA profiles of muscle tissue collected from two sire-line pig breeds (Pietrain, Hampshire). The RNA-seq approach allowed the identification of 627 and 416 differentially expressed genes (DEGs) between pig groups differing in terms of loin weight between Pietrain and Hampshire breeds, respectively. The comparison of miRNA profiles showed differential expression of 57 microRNAs for Hampshire and 34 miRNAs for Pietrain pigs. Next, 43 genes and 18 miRNAs were selected as differentially expressed in both breeds and potentially related to muscle development. According to Gene Ontology analysis, identified DEGs and microRNAs were involved in the regulation of the cell cycle, fatty acid biosynthesis and regulation of the actin cytoskeleton. The most deregulated pathways dependent on muscle mass were the Hippo signalling pathway connected with the TGF-β signalling pathway and controlling organ size via the regulation of ubiquitin-mediated proteolysis, cell proliferation and apoptosis. The identified target genes were also involved in pathways such as the FoxO signalling pathway, signalling pathways regulating pluripotency of stem cells and the PI3K-Akt signalling pathway. The obtained results indicate molecular mechanisms controlling porcine muscle growth and development. Identified genes (SOX2, SIRT1, KLF4, PAX6 and genes belonging to the transforming growth factor beta superfamily) could be considered candidate genes for determining muscle mass in pigs.

Transcriptomic hallmarks of bone remodelling revealed by RNA-Seq profiling in blood of Arabian horses during racing training regime.

The impact of exercises on young developing organisms is still of interest to researchers. Similarly like Thoroughbreds, Arabian horses competing at the race track. The high percent of lameness and loss of days in training are often the result of weakness in the condition of the musculoskeletal system. The objective of the presented study was to identify by RNA-Seq method, the possible skeletal system originating transcriptomic profile in peripheral blood of Arabian horses undergoing race training. Obtained results showed that one of the most significantly deregulated pathway involved in bone homeostasis was those involved in osteoclast differentiation. Among the significantly expressed molecules, we recognized twelve genes potentially involved in the metabolism of the skeletal system: BGLAP, CTSK, TYROBP, PDLIM7, SLC9B2, TWSG1, NOTCH2, IL6ST, VAV3, NFATc1, CLEC5A, TXLNG. The panel of identified genes should be evaluated as candidate biomarkers for bone homeostasis indicators of Arabians performing on race tracks to assess bone remodelling states during training for race track competitions.

Transcriptomic Changes in Broiler Chicken Hypothalamus during Growth and Development

The hypothalamus plays an overarching role that is reflected in the physiological processes observed in the entire organism. The hypothalamus regulates selected metabolic processes and activities of the autonomic nervous system. The avian hypothalamus due to the structural complexity is not well described and has a slightly different function than the mammalian hypothalamus that is the subject of numerous studies. The present study evaluated activities of hypothalamic genes in fast-growing chickens during development (at the 1st day and 3rd and 6th weeks after hatching). The hypothalamic transcriptomes for 3- and 6-week-old cockerels were analysed using an RNA sequencing method in next-generation sequencing (NGS) technology. The differentially expressed gene analysis was conducted using DESeq2 software. In younger 22-day-old cockerels, 389 genes showed higher expression (fold change > 1.5) than that in 45-day-old birds. These genes played a role in several biological processes because they encoded proteins involved in integrin signalling, regulation of hormone levels, camera-type eye development, and blood vessel development. Moreover, surprisingly in the hypothalamus of 3-week-old cockerels, transcripts were identified for proteins involved in both anorexigenic (POMC, NMU) and orexigenic (PMCH, ALDH1A1, LPL, and GHRH) pathways. The RNA-seq results were confirmed by qPCR methods. In summary, the intensive growth of 3-week-old chickens was reflected in hypothalamic activities because the genes associated with the somatotropin axis and regulation of satiety centre showed increased expression.