1
Physiological Genomics: RESEARCH ARTICLE │ General Interest
2 3
Transcriptional Profiling Reveals Gland-specific Differential Expression in
4
the Three Major Salivary Glands of the Adult Mouse
5 6
Xin Gao,1,3¶ Maria S. Oei,1¶ Catherine E. Ovitt,4 Murat Sincan,2* and James E. Melvin1*
7 8
1
Secretory Mechanisms and Dysfunctions Section, 2Office of the Clinical Director, National
9
Institute of Dental and Craniofacial Research, National Institutes of Health, Bethesda, MD 20892
10
3
11
20742
12
4
13
Rochester, NY 14642
14 15
Running Head: Salivary Gland-specific Differences in Gene Expression
Joint Institute for Food Safety and Applied Nutrition, University of Maryland, College Park, MD
Center for Oral Biology and Department of Biomedical Genetics, University of Rochester,
16 17
¶
18
*
19
Address for correspondence: James E. Melvin, D.D.S., Ph.D.
These authors contributed equally to this paper Co-Senior Authors: M. Sincan and J.E. Melvin
20
10 Center Drive
21
Building 10/Room 1N-117B
22
Bethesda, MD 20892
23
Email:
[email protected]
24
Tel: (301) 402-1706
25 26
1
Downloaded from www.physiology.org/journal/physiolgenomics by ${individualUser.givenNames} ${individualUser.surname} (216.165.126.118) on February 28, 2018. Copyright © 2018 American Physiological Society. All rights reserved.
27
ABSTRACT
28 29
Gao X, Oei MS, Ovitt CE, Sincan M, Melvin JE. Transcriptional Profiling Reveals
30
Gland-specific Differential Expression in the Three Major Salivary Glands of the Adult Mouse.
31
Physiol Genomics––RNA-seq was used to better understand the molecular nature of the biological
32
differences among the three major exocrine salivary glands in mammals. Transcriptional profiling
33
found that the adult murine parotid, submandibular and sublingual salivary glands express greater
34
than 14,300 protein-coding genes, and nearly 2,000 of these genes were differentially expressed.
35
Principle component analysis of the differentially expressed genes revealed three distinct clusters
36
according to gland type. The three salivary gland transcriptomes were dominated by a relatively
37
few number of highly expressed genes (6.3%) that accounted for more than 90% of transcriptional
38
output. Of the 912 transcription factors expressed in the major salivary glands, greater than 90%
39
of them were detected in all three glands, while expression for about 2% of them was enriched in
40
an individual gland. Expression of these unique transcription factors correlated with sublingual
41
and parotid specific subsets of both highly expressed and differentially expressed genes. Gene
42
ontology analyses revealed that the highly expressed genes common to all glands were associated
43
with global functions, while many of the genes expressed in a single gland play a major role in the
44
function of that gland. In summary, transcriptional profiling of the three murine major salivary
45
glands identified a limited number of highly expressed genes, differentially expressed genes and
46
unique transcription factors that represent the transcriptional signatures underlying gland-specific
47
biological properties.
48 49
Keywords: RNA-seq; Exocrine Gland; Salivary Gland; Transcription Factors; Gene Expression
50 51
2
Downloaded from www.physiology.org/journal/physiolgenomics by ${individualUser.givenNames} ${individualUser.surname} (216.165.126.118) on February 28, 2018. Copyright © 2018 American Physiological Society. All rights reserved.
52
INTRODUCTION
53 54
THE PRIMARY function of salivary glands is to secrete saliva, a fluid mixture of water,
55
electrolytes and proteins. The saliva secreted by mammals originates largely from three paired
56
exocrine glands, the parotid, submandibular and sublingual salivary glands. Saliva lubricates and
57
moistens the oral cavity to facilitate speech, swallowing and taste, initiates digestion, and inhibits
58
dental decay and other opportunistic infections (1, 46). Although the three major salivary glands
59
are similar in many respects, several functional characteristics are gland-specific. For example,
60
saliva contains more than a thousand proteins, some of which are highly enriched or exclusively
61
expressed by an individual salivary gland (10, 54). Moreover, although the essential ion and water
62
transport proteins have been identified (16, 33), the differences in the fluid secretion mechanism
63
among the different salivary glands have not been fully elucidated (23).
64
Salivary glands are mainly composed of two cell types, i.e. acinar cells that are responsible
65
for secreting most of the proteins, fluid and electrolytes found in saliva, and duct cells that
66
primarily reabsorb the NaCl secreted by the acinar cells. Some of the gland-specific functional
67
properties are associated with the distinctive cell types found in each salivary gland (2, 10, 33, 41,
68
51). In most mammals, parotid acinar cells are exclusively serous in nature, while sublingual acinar
69
cells are primarily mucous cells capped by a few serous cells. The acinar cells of the rodent
70
submandibular gland are seromucous and their duct cells display sexual dimorphism (41).
71
The genetic and molecular basis for the functional differences among the three major
72
salivary glands remains unclear, but can likely be associated with the regulation of salivary gland-
73
specific gene expression by different groups of transcription factors. Transcription factors are the
74
central regulators of gene expression. Approximately 1,600 transcription factors comprise greater
75
than 7% of the protein-coding genes in the mouse genome, which makes this family the single
76
largest gene family (20). Mammalian genes are typically flanked by distinct binding sites for
77
several transcription factors that work in unison to regulate expression of each gene. While
78
globally expressed transcription factors often regulate housekeeping functions (28), others may
79
control specific biological tasks, like organ development (25, 29, 47). Aberrant expression of
80
transcription factors also frequently correlates with tumorigenesis (6, 31). However, expression
81
data are lacking for the specific transcription factors that maintain the adult gland-specific
82
phenotype for each of the three major salivary glands. 3
Downloaded from www.physiology.org/journal/physiolgenomics by ${individualUser.givenNames} ${individualUser.surname} (216.165.126.118) on February 28, 2018. Copyright © 2018 American Physiological Society. All rights reserved.
83
The current study evaluated the transcriptional profiles of the three murine salivary glands
84
using RNA-seq to better understand the molecular and genetic nature of their well-recognized
85
biological differences. We found that the expression patterns of gland-specific protein-coding
86
genes agree with the adult male mouse submandibular gland transcriptome (12) and the salivary
87
gland secretion literature, validating our RNA-seq data. Transcriptome variation across the three
88
murine salivary glands correlated with the differential co-expression of a limited number of
89
transcription factors and salivary gland-specific protein-coding genes.
90 91
MATERIALS AND METHODS
92 93
General Methods. Chemicals were purchased from Sigma-Aldrich (St Louis, MO, USA)
94
unless otherwise stated. Experiments were performed using 12-weeks old male and female
95
C57BL/6J mice (Jackson Laboratory, ME) that were housed in pathogen-free, micro-isolator cages
96
with free access to laboratory chow and water with a 12-hour light/dark cycle. Animal procedures
97
were approved by the Animal Care and Use Committee of the National Institute of Dental and
98
Craniofacial Research, National Institutes of Health (ASP 16-802).
99
RNA preparation and next-generation sequencing. Six adult mice (3 male and 3 female)
100
were euthanized by exposure to CO2, followed by cervical dislocation, and the parotid,
101
submandibular and sublingual glands immediately removed and placed in RNAlater Stabilization
102
Solution (Thermo Fisher Scientific, Waltham, MA, catalog #AM7020) at 4°C. The submandibular
103
and sublingual glands are encapsulated, and thus, were easily separated from the surrounding
104
tissue, while the parotid gland is more diffuse and was removed with the aid of a dissecting
105
microscope to exclude surrounding connective tissue. A single parotid and submandibular gland
106
were used from each animal. The sublingual glands are considerably smaller than the other two
107
major salivary glands, thus the left and right glands were pooled from each animal. The 18 gland
108
samples were shipped overnight to Otogenetics Corporation (Atlanta, GA). Total RNA was
109
extracted using the RNeasy Micro Kit (Qiagen, Valencia, CA, catalog #74004). The integrity and
110
purity of total RNA were assessed using Agilent Bioanalyzer or TapeStation and OD260/280.
111
cDNA was generated from high quality total RNA using the SMARTer PCR cDNA Synthesis Kit
112
with modified oligo(dT) primers (Clontech Laboratories, Inc., Mountain View, CA, catalog
113
#634926), and adapters were removed by digestion with RsaI. The resulting 1-2 μg of cDNA were 4
Downloaded from www.physiology.org/journal/physiolgenomics by ${individualUser.givenNames} ${individualUser.surname} (216.165.126.118) on February 28, 2018. Copyright © 2018 American Physiological Society. All rights reserved.
114
fragmented using Bioruptor (Diagenode Inc., Denville, NJ) and profiled using Agilent Bioanalyzer
115
or TapeStation. Illumina libraries were made from qualified fragmented cDNA using the
116
SPRIworks HT Reagent Kit (Beckman Coulter, Inc., Indianapolis, IN, catalog #B06938). The
117
quality, quantity, and size distribution of the Illumina libraries were determined using Agilent
118
Bioanalyzer or Tapestation. The libraries were then submitted for Illumina HiSeq2500 sequencing
119
according to the standard operation with v1 or v2 chemistry, and 100-106 nucleotide (nt) paired
120
end reads were generated and checked for data quality using FastQC (Babraham Institute,
121
Cambridge, UK). Each of the 18 sample libraries were sequenced at approximately 40 million
122
reads.
123
Analysis and Comparison of RNA-Seq Data.
The data generated by Otogenetics
124
Corporation were subjected to analysis using the platform provided by DNAnexus (Mountain
125
View, CA) and then: (1) mapped against the mouse mm10 reference genome with STAR 2.4.0j;
126
(2) measured for expression levels of genes and transcripts with Cufflinks 2.2.1; and (3) compared
127
across the three major salivary glands with Cufflinks/Cuffdiff 2.2.1 for significance (q500-fold more
171
abundant in the SLG and SMG than in the PG, while Mucin 19/Muc19 was expressed at a 2,000-
172
fold greater level in the SLG than in the PG and SMG.
173
To further validate our RNA-seq data, we compared the adult male SMG transcriptome
174
generated in the present study with the adult male SMG transcriptome previously reported (12) 6
Downloaded from www.physiology.org/journal/physiolgenomics by ${individualUser.givenNames} ${individualUser.surname} (216.165.126.118) on February 28, 2018. Copyright © 2018 American Physiological Society. All rights reserved.
175
using the same age and mouse background (12-weeks old, C57BL/6). Gluck et al. identified 246
176
protein-coding genes enriched in the adult male SMG transcriptome, 220 of which were detected
177
in our male SMG transcriptome (Supplemental Table 3). The 90% overlap in detected protein-
178
coding genes in the two male SMG transcriptomes confirms the reproducibility of the RNA-seq
179
method.
180
Highly Expressed Genes in the Three Major Murine Salivary Glands. Among the
181
expressed protein-coding genes, we selected the ones most highly expressed for further
182
examination, as these genes generally hold keys to understanding the genetic and biochemical
183
foundation of a given organ system. To identify the highly-expressed genes (HEGs) in each gland,
184
the sorted FPKMs were individually graphed on a logarithmic scale and a log FPKM cutoff value
185
of 4 was used, as revealed by the turning points on the graphs (blue lines; Supplemental Figure 2).
186
We observed in each of the major salivary glands about 700 HEGs with log FPKM values greater
187
than 4 (Figure 1). Allowing that the three major salivary glands have similar biological functions,
188
it was therefore not surprising that 685 of the 901 total HEGs (about 76%) in the PG, SLG and
189
SMG presented with high expression levels in at least two of the three salivary glands (Figure 1;
190
list of HEGs given in Supplemental Table 4). In contrast, 7.8%, 11.0% and 5.2% of the HEGs (PG,
191
SLG and SMG, respectively) displayed enriched expression in a single gland, indicating that the
192
latter group of HEGs contributes to the unique transcriptomic signature of each gland. The 901
193
different HEGs in the three salivary glands represent only about 6.3% of the total number of
194
expressed genes (14,371), but these highly expressed genes dominated expression, contributing
195
96.6%, 87.8%, and 93.4% of the transcriptional outputs in PG, SLG, and SMG, respectively (total
196
HEG FPKMs divided by total FPKMs for each gland).
197
Gene Ontology (GO) enrichment analyses (refinement P value