Subscriber access provided by REGIS UNIV
Omics Technologies Applied to Agriculture and Food
Comparative Analysis of Long Noncoding RNAs Expressed During Intramuscular Adipocytes Adipogenesis in Fat-type and Lean-type Pigs Yunmei Sun, Xiaochang Chen, Jin Qin, Shuge Liu, Rui Zhao, Taiyong Yu, Guiyan Chu, Gongshe Yang, and Pang Wei-jun J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b04243 • Publication Date (Web): 19 Oct 2018 Downloaded from http://pubs.acs.org on October 20, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 30
Journal of Agricultural and Food Chemistry
1
Comparative Analysis of Long Noncoding RNAs Expressed During Intramuscular
2
Adipocytes Adipogenesis in Fat-type and Lean-type Pigs
3
Yunmei Sun#, Xiaochang Chen#, Jin Qin, Shuge Liu, Rui Zhao, Taiyong Yu, Guiyan Chu,
4
Gongshe Yang*, Weijun Pang*
5
Laboratory of Animal Fat Deposition & Muscle Development, College of Animal Science
6
and Technology, Northwest A&F University, Yangling Shaanxi 712100, China.
7
#These
8
* Correspondence:
9
Weijun Pang, Ph.D.
authors contributed equally to this work.
10
Address: No. 22 Xinong Road, Yangling, Shaanxi Province 712100, China
11
Tel: 86-29-87091017
12
Fax: 86-29-87092430
13
E-mail: pwj1226@nwsuaf.edu.cn
14
Gongshe Yang, Ph.D.
15
Address: No. 22 Xinong Road, Yangling, Shaanxi Province 712100, China
16
Tel: 86-29-87091017
17
Fax: 86-29-87092312
18
E-mail: gsyang999@hotmail.com
19
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
20
ABSTRACT: The meat quality of local breed pigs is more tender and juicier than the
21
imported varieties. The important reason is that the intramuscular fat content is high. Even
22
though modest sequence conservation and evolution, the expression pattern and function
23
of long non-coding RNAs (lncRNAs) seem to be conserved. In spite of that, analysis of
24
lncRNAs associated with intramuscular fat development remains unknown to us in porcine.
25
Here, we systematically investigated lncRNAs of intramuscular adipocytes of fat local
26
Bamei pig and lean Large White pig to consider the function of lncRNAs on intramuscular
27
fat development. We selected three piglets of both breeds separately to isolate
28
intramuscular preadipocytes and performed RNA sequencing across four stages (0, 2, 4
29
and 8 d) during the intramuscular preadipocytes differentiation, and identified 1932
30
lncRNAs (760 novel). In addition, we have screened lnc_000414 closely related to fat
31
synthesis. This lncRNA function as an inhibitor in the proliferation of porcine
32
intramuscular adipocytes. These novel findings will provide new targets for improving
33
pork quality and making pig breeding better.
34
KEYWORDS: pig, lncRNA-seq, intramuscular adipocyte, adipogenesis, lnc_000414
35
ACS Paragon Plus Environment
Page 2 of 30
Page 3 of 30
Journal of Agricultural and Food Chemistry
36
INTRODUCTION
37
Now, with the improvement of living standards, people pay more attention to meat
38
quality. Pork, as a major source of meat, is closely related to human health. Its quality is
39
strongly influenced by fat distribution and deposition1. Therefore, how to improve meat
40
production and quality becomes an important issue for porcine breeding scientists.
41
Coincidently, more and more studies demonstrated that intramuscular fat (IMF) was
42
closely connected with the quality of meat, affecting properties such as flavor, water-
43
holding capacity and tenderness2. Thus, in the past few years, researches on IMF have
44
increased frequently in the field of porcine fat deposition3, suggesting it is essential to study
45
candidate genes controlling IMF deposition.
46
The regulation processes of proliferation, differentiation and tissue development in
47
porcine IMF are actually a result of the interaction between inheritance, epigenetics and
48
signaling molecules4-6. And most importantly, long noncoding RNA (lncRNA) plays a
49
major role in these processes. LncRNAs is a type of non-coding RNA longer than 200 nt7,
50
and previous studies have shown that it could involve in many biological processes,
51
including transcription level and post-transcription level, although it has only weak ability
52
to encode proteins. They were associated with many life processes, including cell
53
proliferation, differentiation and apoptosis by regulating the expression of their target
54
genes8-10. There were some evidences demonstrated that lncRNAs play critical roles during
55
adipogenesis in vitro as well as in vivo. As a matter of fact, more and more studies have
56
demonstrated that lncRNAs could regulate adipogenesis. For instance, PU.1 AS lncRNA
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
57
elevates adipogenic differentiation of preadipocyte through inhibiting the translation of
58
PU.1 mRNA in different animal models
59
are thought to be potential therapy targets for weight loss, suggesting that the associated
60
lncRNAs 11may have an effect on the formation of brown and beige adipocytes, and
61
contribute to the treatment of obesity13.
62
11,12.
Nevertheless, brown and beige adipocytes
Bamei pig, a Chinese native fat-type pig breed, is known for good meat quality and 14.
63
relative high lean meat percentage
Thus, it is a useful model for studying muscle
64
development and meat yield of indigenous pig breed. Although, there are also some reports
65
on differentially expressed lncRNAs between fat-type pig and lean-type pig breeds15. There
66
are also more new lncRNAs to consider in other breeds especially indigenous fat-type pig
67
breed. In this study, we obtained the lncRNAs expression variation in preadipocytes
68
between Bamei pig and Large White pig at different differentiation stages by lncRNA high-
69
throughput sequencing, and aimed to explore the lncRNAs involved in intramuscular fat
70
development of Chinese indigenous pigs. Our findings will be useful to expand the
71
knowledge of regulating meat quality by lncRNAs and be beneficial for the genetic
72
improvement of porcine meat quality traits.
73
MATERIALS AND METHODS
74
Animals and samples
75
The male Chinese native Bamei piglets (three-day-old, 1.2 kg) were purchased from
76
Bamei original breeding filed of Huzhu, Gansu, and the male Large White piglets (three-
77
day-old, 1.5 kg) were from the experimental station of Northwest A&F University. These
ACS Paragon Plus Environment
Page 4 of 30
Page 5 of 30
Journal of Agricultural and Food Chemistry
78
piglets lived in similar rearing conditions.
79
Intramuscular preadipocytes were isolated from longissimus dorsi muscle (LDM) of
80
3-day-old piglets as described previously16. In briefly, we digested the muscle tissues by
81
Dulbeccos modified Eagles medium (DMEM; Gibco, Carlsbad, CA) with 0.1%
82
collagenase I (270 U/mg; Gibco, Carlsbad, CA) at 37°C for 1.5h. The digest samples were
83
sequentially filtered through 70 and 200 mesh filters to separate the cells. Next, we used
84
DMEM/F12 medium to wash the cells thrice. Cells were seeded in dishes containing
85
DMEM/F12 medium added 10% fetal bovine serum (Gibco, Australia). After 1 h, we
86
rinsed off unattached cells with DMEM/F12 medium. When density of cells reaches to
87
approximately100 %, they were induced to differentiate with DMEM/F12 containing 10%
88
FBS, 5μg/ml (872 nM) insulin, 1 μM dexamethasone and 0.5 mM isobutyl methylxanthine
89
(IBMX, Sigma–Aldrich, St. Louis, MO) for 2 days. The medium with 10% FBS and
90
5μg/ml insulin was used to retain differentiation. Finally, we harvested the cells at0, 2, 4
91
and 8 d for subsequent research.
92
In our study, all the experimental procedures with pigs are based on the Experimental
93
Animal Manage Committee of Northwest A&F University (2011-31101684).
94
RNA isolation, library preparation, and sequencing
95
We used Trizol reagent (Invitrogen, USA) to extract the total RNA from all cell
96
samples, and the method was followed manufacturer’s instructions. The quality of RNA
97
and concentration was detected by the NanoDrop 2000 (Thermo Scientific, Waltham, MA,
98
USA). After that, the qualified RNA was stored at -80°C for later use.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
99
Each cell sample was extracted with 3μg of RNA for preparation of sequencing
100
samples. Firstly, we used Epicentre Ribo-zero™ rRNA Removal Kit (Epicentre, USA) and
101
ethanol to remove ribosomal RNA and free residue separately Subsequently, sequencing
102
libraries were generated using the rRNA-depleted RNA by NEBNext® Ultra™ Directional
103
RNA Library Prep Kit for Illumina® (NEB, USA). Briefly, fragmentation Buffer was
104
added to the obtained mRNA to make the fragment into a short fragment. Then, using the
105
fragmented mRNA as a template, the first strand of cDNA was synthesized with random
106
hexamers and buffer, dNTPs, RNase H and DNA polymerase I synthesizes the second
107
strand of cDNA, which is purified by QiaQuick PCR kit and eluted with EB buffer. The
108
end-repair, base A is added, the sequencing linker is added, and the fragment of interest is
109
recovered by agarose gel electrophoresis, and PCR is performed. Amplification was
110
performed to complete the entire library preparation work, and the constructed library was
111
sequenced using Illumina HiSeq2000.
112
The remaining overhangs were converted to blunt ends by exonuclease/polymerase
113
activity. After adenylation at the 3' end of the DNA fragment, a NEBNext adaptor with a
114
hairpin loop structure was ligated to prepare for hybridization. The cDNA fragment of 150
115
to 200 bp was preferentially selected, and a cDNA fragment library was purified using an
116
AMPure XP system (Beckman Coulter, Beverly, USA). The appropriately sized cDNA
117
was then ligated at 37 °C for 15 minutes following 95 °C for 5 minutes by USER Enzyme
118
(NEB, USA). Subsequently, PCR was performed by using Phusion High-Fidelity DNA
119
polymerase, Universal PCR primers and Index (X) Primer.
ACS Paragon Plus Environment
Page 6 of 30
Page 7 of 30
120
Journal of Agricultural and Food Chemistry
Quality control
121
Firstly, raw data is measured by each internal script in the fastq format. Moreover,
122
with the strict screening criteria, the clean data from raw data were collected.
123
Synchronously, the related indicators of clean data were calculated, including Q20, Q30
124
and GC content17. And the high quality clean data was used for further downstream
125
analyzes.
126
Transcriptome assembly
127
The Cufflinks (v2.1.1) and Scripture (beta2) were used for assembling the mapped
128
reads for samples. In order to determine exons connectivity, two different approaches were
129
used to splice reads in these two methods.
130
quantify the expression level of the expression data by using a probabilistic model in a
131
given locus. And Scripture could distinguish expressed loci and used spliced reads to
132
assemble expressed segments though a statistical segmentation model. We used a variety
133
of methods to predict the coding ability of transcripts including PFAM database19, Coding-
134
Non-Coding Index (CNCI) software20 , and Coding Potential Calculator (CPC)18, if the
135
transcript cannot pass the test, it will be excluded. On one hand, all parameters were set as
136
default in Scripture analysis, and on another hand, Cufflinks analysis was based on the
137
program named ‘min-frags-per-transfrag=0’ and ‘--library-type’. Besides these two
138
parameters, other parameters were run with default.
139
Quantification of gene expression level
140
In briefly, Cufflinks could assemble and
Cuffdiff (v2.1.1) was used to calculate FPKMs of both lncRNAs and coding genes in
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
141
each sample. Gene FPKMs were counted by summing the FPKMs of transcripts in each
142
group. FPKM computed the number of fragments in each exon isoform. This was based on
143
the length of the fragments and reads count mapped to this fragment.
144
Differential expression analysis
145
Cuffdiff provided statistical routines for determining differential expression in digital
146
transcript or gene expression data using a model which was based on the negative binomial
147
distribution. Transcripts were considered to be differentially expressed, if P-adjust