scholarly journals Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differentially expressed transcripts in sheep tail fat

PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e12454
Author(s):  
Zehu Yuan ◽  
Ling Ge ◽  
Jingyi Sun ◽  
Weibo Zhang ◽  
Shanhe Wang ◽  
...  

Background Nowadays, both customers and producers prefer thin-tailed fat sheep. To effectively breed for this phenotype, it is important to identify candidate genes and uncover the genetic mechanism related to tail fat deposition in sheep. Accumulating evidence suggesting that post-transcriptional modification events of precursor-messenger RNA (pre-mRNA), including alternative splicing (AS) and alternative polyadenylation (APA), may regulate tail fat deposition in sheep. Differentially expressed transcripts (DETs) analysis is a way to identify candidate genes related to tail fat deposition. However, due to the technological limitation, post-transcriptional modification events in the tail fat of sheep and DETs between thin-tailed and fat-tailed sheep remains unclear. Methods In the present study, we applied pooled PacBio isoform sequencing (Iso-Seq) to generate transcriptomic data of tail fat tissue from six sheep (three thin-tailed sheep and three fat-tailed sheep). By comparing with reference genome, potential gene loci and novel transcripts were identified. Post-transcriptional modification events, including AS and APA, and lncRNA in sheep tail fat were uncovered using pooled Iso-Seq data. Combining Iso-Seq data with six RNA-sequencing (RNA-Seq) data, DETs between thin- and fat-tailed sheep were identified. Protein protein interaction (PPI) network, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were implemented to investigate the potential functions of DETs. Results In the present study, we revealed the transcriptomic complexity of the tail fat of sheep, result in 9,001 potential novel gene loci, 17,834 AS events, 5,791 APA events, and 3,764 lncRNAs. Combining Iso-Seq data with RNA-Seq data, we identified hundreds of DETs between thin- and fat-tailed sheep. Among them, 21 differentially expressed lncRNAs, such as ENSOART00020036299, ENSOART00020033641, ENSOART00020024562, ENSOART00020003848 and 9.53.1 may regulate tail fat deposition. Many novel transcripts were identified as DETs, including 15.527.13 (DGAT2), 13.624.23 (ACSS2), 11.689.28 (ACLY), 11.689.18 (ACLY), 11.689.14 (ACLY), 11.660.12 (ACLY), 22.289.6 (SCD), 22.289.3 (SCD) and 22.289.14 (SCD). Most of the identified DETs have been enriched in GO and KEGG pathways related to extracellular matrix (ECM). Our result revealed the transcriptome complexity and identified many candidate transcripts in tail fat, which could enhance the understanding of molecular mechanisms behind tail fat deposition.

2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Kai Xing ◽  
Kejun Wang ◽  
Hong Ao ◽  
Shaokang Chen ◽  
Zhen Tan ◽  
...  

Abstract Fatness traits are important in pigs because of their implications for fattening efficiency, meat quality, reproductive performance and immunity. Songliao black pigs and Landrace pigs show important differences in production and meat quality traits, including fatness and muscle growth. Therefore, we used a high-throughput massively parallel RNA-seq approach to identify genes differentially expressed in backfat tissue between these two breeds (six pigs in each). An average of 37.87 million reads were obtained from the 12 samples. After statistical analysis of gene expression data by edgeR, a total of 877 differentially expressed genes were detected between the two pig breeds, 205 with higher expression and 672 with lower expression in Songliao pigs. Candidate genes (LCN2, CES3, DGKB, OLR1, LEP, PGM1, PCK1, ACACB, FADS1, FADS2, MOGAT2, SREBF1, PPARGC1B) with known effects on fatness traits were included among the DEGs. A total of 1071 lncRNAs were identified, and 85 of these lncRNAs were differentially expressed, including 53 up-regulated and 32 down-regulated lncRNAs, respectively. The differentially expressed genes and lncRNAs involved in glucagon signaling pathway, glycolysis/gluconeogenesis, insulin signaling pathway, MAPK signaling pathway and so on. Integrated analysis potential trans-regulating or cis-regulating relation between DEGs and DE lncRNAs, suggested lncRNA MSTRG.2479.1 might regulate the expressed level of VLDLR affecting porcine fat metabolism. These results provide a number of candidate genes and lncRNAs potentially involved in porcine fat deposition and provide a basis for future research on the molecular mechanisms underlying in fat deposition.


2018 ◽  
Vol 143 (1) ◽  
pp. 56-66
Author(s):  
Yu Bai ◽  
Ying Zhou ◽  
Xiaoqing Tang ◽  
Yu Wang ◽  
Fangquan Wang ◽  
...  

The appropriate timing of bolting and flowering is one of the keys to the reproductive success of Isatis indigotica. Several flowering regulatory pathways have been reported in plant species, but we know little about flowering regulatory in I. indigotica. In the present study, we performed RNA-seq and annotated I. indigotica transcriptome using RNA from five tissues (leaves, roots, flowers, fruit, and stems). Illumina sequencing generated 149,907,857 high-quality clean reads and 124,508 unigenes were assembled from the sequenced reads. Of these unigenes, 88,064 were functionally annotated by BLAST searches against the public protein databases. Functional classification and annotation assigned 55,991 and 23,072 unigenes to 52 gene ontology (GO) terms and 25 clusters of orthologous group (COG) categories, respectively. A total of 19,927 unigenes were assigned to 124 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and 80 candidate genes related to plant circadian rhythm were identified. We also identified a number of differentially expressed genes (DEG) and 91 potential bolting and flowering-related genes from the RNA-seq data. This study is the first to identify bolting and flowering-related genes based on transcriptome sequencing and assembly in I. indigotica. The results provide foundations for the exploration of flowering pathways in I. indigotica and investigations of the molecular mechanisms of bolting and flowering in Brassicaceae plants.


2019 ◽  
Author(s):  
Fuping Zhang ◽  
Liangting Tang ◽  
Xueqin Ran ◽  
Ning Mao ◽  
Yiqi Ruan ◽  
...  

AbstractBackground/AimsLitter size is one of the most important reproductive traits in pig breeding, which is affected by multiple genes and the environment. Ovaries are the most important reproductive organs and have a profound impact on the reproduction efficiency. Therefore, genetic differences in the ovaries may contribute to the observed differences in litter size. Although QTLs and candidate genes have been reported to affect the litter size in many pig breeds, however, the findings cannot elucidate the marked differences of the reproductive traits between breeds. The aim of present work is to elucidate the mechanisms of the differences for the reproductive traits and identify candidate genes associated with litter size in Xiang pig breed.MethodsThe changes in ovary transcriptome and alternative splicing were investigated at estrus between Xiang pigs with large and small litter size by RNA-seq technology. The RNA-seq results were confirmed by RT-qPCR method.ResultsWe detected 16,219 - 16,285 expressed genes and 12 types of alternative splicing (AS) events in Xiang pig samples. A total of 762 differentially expressed genes were identified by XL (Xiang pig group with larger litter size) vs XS (Xiang pig group with small litter size) sample comparisons. A total of 34 genes were upregulated and 728 genes were downregulated in XL ovary samples compared with the XS samples. Alternative splicing (AS) rates in XL samples were slightly lower than that observed in XS samples. Most of differentially expressed genes were differentially regulated on AS level. Eleven candidate genes were potentially identified to be related to Xiang pig fecundity and litter size, which may be closely related to the gonad development, oocyte maturation or embryo quality.ConclusionThe significant changes in the expression of the protein-coding genes and the level of alternative splicing in estrus ovarian transcriptome between XL and XS groups probably are the molecular mechanisms of phenotypic variation in litter size.


Agriculture ◽  
2021 ◽  
Vol 12 (1) ◽  
pp. 15
Author(s):  
Tianyou Bai ◽  
Benmeng Liang ◽  
Yuhetian Zhao ◽  
Jiangang Han ◽  
Yabin Pu ◽  
...  

The hair follicle is a complex mini-organ in the skin that undergoes organ induction, morphogenesis, and regeneration. However, the accurate molecular mechanism of skin and hair diversity regulation is still elusive. The sheep is an animal model that can be used to further explore the mechanisms of skin and hair diversity. In this study, we carried out a transcriptomic analysis of the mRNA expression in the skin of Xinji fine-wool sheep at different growth stages (3 and 12 months old) and 12-month-old Tan sheep and explored the transcripts’ relationship with hair follicle growth. A total of 1327 mRNAs and 67 transcription factors were identified to be differentially expressed in the different breeds and during different periods of skin development. The differentially expressed genes were enriched in GO terms represented by system development, multicellular organism development, animal organ development, and skin development, and three KEGG pathways typified those governing differences in skin structure. Combining protein–protein interaction networks of skin development (GO:0043588) and functional annotation, nine important candidate genes, namely, LAMA5, OVOL1, SRF, DHCR24, NGFR, SMO, CDSN, HOXC13, and KDF1, and many core genes with minor effects were confirmed to be associated with hair follicle development. Furthermore, members of the zf-C2H2 and homeobox transcription factor families, which were identified to play a crucial role in producing finer and denser wool, were mainly upregulated in 12-month-old Xinji fine-wool sheep when compared with expression in 12-month-old Tan sheep and 3-month-old Xinji fine-wool sheep. This study revealed the major–minor gene interactions in the developmental pathway and provided ideas for an in-depth understanding of the genetic structure and gene regulation in the skin/hair growth process.


2018 ◽  
Vol 2018 ◽  
pp. 1-8 ◽  
Author(s):  
Songbai Yang ◽  
Xiaolong Zhou ◽  
Yue Pei ◽  
Han Wang ◽  
Ke He ◽  
...  

Estrus is an important factor for the fecundity of sows, and it is involved in ovulation and hormone secretion in ovaries. To better understand the molecular mechanisms of porcine estrus, the expression patterns of ovarian mRNA at proestrus and estrus stages were analyzed using RNA sequencing technology. A total of 2,167 differentially expressed genes (DEGs) were identified (P≤0.05, log2  Ratio≥1), of which 784 were upregulated and 1,383 were downregulated in the estrus compared with the proestrus group. Gene Ontology (GO) enrichment indicated that these DEGs were mainly involved in the cellular process, single-organism process, cell and cell part, and binding and metabolic process. In addition, a pathway analysis showed that these DEGs were significantly enriched in 33 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including cell adhesion molecules, ECM-receptor interaction, and cytokine-cytokine receptor interaction. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) confirmed the differential expression of 10 selected DEGs. Many of the novel candidate genes identified in this study will be valuable for understanding the molecular mechanisms of the sow estrous cycle.


2012 ◽  
Vol 111 (suppl_1) ◽  
Author(s):  
Emma L Robinson ◽  
Syed Haider ◽  
Hillary Hei ◽  
Richard T Lee ◽  
Roger S Foo

Heart failure comprises of clinically distinct inciting causes but a consistent pattern of change in myocardial gene expression supports the hypothesis that unifying biochemical mechanisms underlie disease progression. The recent RNA-seq revolution has enabled whole transcriptome profiling, using deep-sequencing technologies. Up to 70% of the genome is now known to be transcribed into RNA, a significant proportion of which is long non-coding RNAs (lncRNAs), defined as polyribonucleotides of ≥200 nucleotides. This project aims to discover whether the myocardium expression of lncRNAs changes in the failing heart. Paired end RNA-seq from a 300-400bp library of ‘stretched’ mouse myocyte total RNA was carried out to generate 76-mer sequence reads. Mechanically stretching myocytes with equibiaxial stretch apparatus mimics pathological hypertrophy in the heart. Transcripts were assembled and aligned to reference genome mm9 (UCSC), abundance determined and differential expression of novel transcripts and alternative splice variants were compared with that of control (non-stretched) mouse myocytes. Five novel transcripts have been identified in our RNA-seq that are differentially expressed in stretched myocytes compared with non-stretched. These are regions of the genome that are currently unannotated and potentially are transcribed into non-coding RNAs. Roles of known lncRNAs include control of gene expression, either by direct interaction with complementary regions of the genome or association with chromatin remodelling complexes which act on the epigenome.Changes in expression of genes which contribute to the deterioration of the failing heart could be due to the actions of these novel lncRNAs, immediately suggesting a target for new pharmaceuticals. Changes in the expression of these novel transcripts will be validated in a larger sample size of stretched myocytes vs non-stretched myocytes as well as in the hearts of transverse aortic constriction (TAC) mice vs Sham (surgical procedure without the aortic banding). In vivo investigations will then be carried out, using siLNA antisense technology to silence novel lncRNAs in mice.


2018 ◽  
Author(s):  
Peter A. Combs ◽  
Joshua J. Krupp ◽  
Neil M. Khosla ◽  
Dennis Bua ◽  
Dmitri A. Petrov ◽  
...  

AbstractPheromones known as cuticular hydrocarbons are a major component of reproductive isolation in Drosophila. Individuals from morphologically similar sister species produce different sets of hydrocarbons that allow potential mates to identify them as a suitable partner. In order to explore the molecular mechanisms underlying speciation, we performed RNA-seq in F1 hybrids to measure tissue-specific cis-regulatory divergence between the sister species D. simulans and D. sechellia. By focusing on cis-regulatory changes specific to female oenocytes, we rapidly identified a small number of candidate genes. We found that one of these, the fatty acid elongase eloF, broadly affects both the complement of hydrocarbons present on D. sechellia females and the propensity of D. simulans males to mate with those females. In addition, knockdown of eloF in the more distantly related D. melanogaster led to a similar shift in hydrocarbons as well as lower interspecific mate discrimination by D. simulans males. Thus, cis-regulatory changes in eloF appear to be a major driver in the sexual isolation of D. simulans from multiple other species. More generally, our RNA-seq approach proved to be far more efficient than QTL mapping in identifying candidate genes; the same framework can be used to pinpoint cis-regulatory drivers of divergence in a wide range of traits differing between any interfertile species.


2020 ◽  
Author(s):  
Dawei Zhang ◽  
Wenjing Wu ◽  
Xin Huang ◽  
Ke Xu ◽  
Cheng Zheng ◽  
...  

Abstract Background: Chinese domestic pig breeds are reputed for pork quality, but their low ratio of lean-to-fat carcass weight decreases production efficiency. A better understanding of the genetic regulation network of SC fat tissue is necessary for the rational selection of Chinese domestic pig breeds. In the present study, SC adipocytes were isolated from Jiaxing Black pigs (a Chinese indigenous pig breed with redundant SC fat deposition) and Large White pigs (a lean-type pig breed with relatively low SC fat deposition) and the expression profiles of mRNAs and lncRNAs were compared by RNA-seq analysis to identify biomarkers correlated with the differences of SC fat deposition between the two breeds.Results: A total of 3,371 differentially expressed genes (DEGs) and 1,182 differentially expressed lncRNAs (DELs) were identified in SC adipocytes between Jiaxing Black (JX) and Large White (LW) pigs, which included 797 upregulated mRNAs, 2,574 downregulated mRNAs, 461 upregulated lncRNAs and 721 downregulated lncRNAs. Gene Ontology and KEGG pathway analyses revealed that the DEGs and DELs were mainly involved in the immune response, cell fate determination, PI3K-Akt signaling pathway and MAPK signaling pathway, which are known to be related to adipogenesis and lipid metabolism. The expression levels of DEGs and DELs according to the RNA-seq data were verified by quantitative PCR, which showed 81.8% consistency. The differences in MAPK pathway activity between JX and LW pigs was confirmed by western blot analysis, with <100-fold elevated p38 phosphorylation in JX pigs.Conclusions: This study offers a detailed characterization of mRNAs and lncRNAs in fat- and lean-type pig breeds. The activity of the MAPK signaling pathway was found to be associated with subcutaneous adipogenesis. These results greatly enhance our understanding of the molecular mechanisms regulating SC fat deposition in pigs.


2021 ◽  
Author(s):  
Zhihui Wang ◽  
Liying Yan ◽  
Yuning Chen ◽  
Xin Wang ◽  
Dongxin Huai ◽  
...  

Abstract Seed weight is a major target of peanut breeding as an important component of seed yield. However, relatively little is known about QTLs and candidate genes associated with seed weight in peanut. In this study, three major QTLs on chromosomes A05, B02 and B06 were determined by applying NGS-based QTL-seq approach for a RIL population. These three QTL regions have been successfully narrowed down through newly developed SNP and SSR markers based on traditional QTL mapping. Among these three QTL regions, qSWB06.3 exhibited stable expression with large contribution to phenotypic variance across all environments. Furthermore, RNA-seq were applied for early, middle and late stages of seed development, and differentially expression genes (DEGs) were identified in ubiquitin-proteasome pathway, serine/threonine protein pathway and signal transduction of hormones and transcription factors. Notably, DEGs at early stage were majorly related to regulating cell division, whereas DEGs at middle and late stages were mainly associated with cell expansion during seed development. Through integrating SNP variation, gene expression and functional annotation, candidate genes related to seed weight in qSWB06.3 were predicted and distinct expression pattern of those genes were exhibited using qRT-PCR. In addition, KASP-markers in qSWB06.3 were successfully validated in diverse peanut varieties and the alleles of parent Zhonghua16 in qSWB06.3 was associated with high seed weight. This suggested that qSWB06.3 was reliable and the markers in qSWB06.3 could be deployed in marker-assisted breeding to enhance seed weight. This study provided insights into the understanding of genetic and molecular mechanisms of seed weight in peanut.


2020 ◽  
Vol 34 ◽  
pp. 205873842097630
Author(s):  
Li Jiang ◽  
Mengmeng Zhang ◽  
Sixue Wang ◽  
Yuzhen Xiao ◽  
Jingni Wu ◽  
...  

The current study intended to explore the interaction of the long non-coding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) under the background of competitive endogenous RNA (ceRNA) network in endometriosis (EMs). The differentially expressed miRNAs (DEmiRs), differentially expressed lncRNA (DELs), and differentially expressed genes (DEGs) between EMs ectopic (EC) and eutopic (EU) endometrium based on three RNA-sequencing datasets (GSE105765, GSE121406, and GSE105764) were identified, which were used for the construction of ceRNA network. Then, DEGs in the ceRNA network were performed with Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and protein-protein interaction (PPI) analysis. Besides, the DEmiRs in the ceRNA network were validated in GSE124010. And the target DELs and DEGs of verified DEmiRs were validated in GSE86534. The correlation of verified DEmiRs, DEGs, and DELs was explored. Moreover, gene set enrichment analysis (GSEA) was applied to investigate the function of verified DEmiRs, DEGs, and DELs. Overall, 1352 DEGs and 595 DELs from GSE105764, along with 27 overlapped DEmiRs between GSE105765 and GSE121406, were obtained. Subsequently, a ceRNA network, including 11 upregulated and 16 downregulated DEmiRs, 7 upregulated and 13 downregulated DELs, 48 upregulated and 46 downregulated DEGs, was constructed. The GO and KEGG pathway analysis showed that this ceRNA network probably was associated with inflammation-related pathways. Furthermore, hsa-miR-182-5p and its target DELs (LINC01018 and SMIM25) and DEGs (BNC2, CHL1, HMCN1, PRDM16) were successfully verified in the validation analysis. Besides, hsa-miR-182-5p was significantly negatively correlated with these target DELs and DEGs. The GSEA analysis implied that high expression of LINC01018, SMIM25, and CHL1, and low expression of hsa-miR-182-5p would activate inflammation-related pathways in endometriosis EU samples. LINC01018 and SMIM25 might sponge hsa-miR-182-5p to upregulate downstream genes such as CHL1 to promote the development of endometriosis.


Sign in / Sign up

Export Citation Format

Share Document