Comparative Analysis of Long Noncoding RNAs Expressed During

2 days ago - ... gas emissions, climate change. The world faces a final opportunity to limit the extreme impacts of global warming, a key internati...
2 downloads 0 Views 1002KB Size
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