Clinical characteristics of the study population
The clinical characteristics of the LEMS and control groups are presented in Table 1.
Table 1 Clinical characteristics of patients with LEMS and healthy controls
Subject ID
|
Age, years
|
Sex
|
Cerebellar signs
|
Presence of a tumor
|
Treated
|
LEMS1
|
56
|
Male
|
Negative
|
SCLC
|
Negative
|
LEMS2
|
55
|
Male
|
Negative
|
SCLC
|
Negative
|
LEMS3
|
54
|
Male
|
Negative
|
SCLC
|
Negative
|
N1
|
52
|
Male
|
Negative
|
Negative
|
Negative
|
N2
|
53
|
Male
|
Negative
|
Negative
|
Negative
|
N3
|
54
|
Male
|
Negative
|
Negative
|
Negative
|
LEMS: Lambert-Eaton myasthenic syndrome; N: number, healthy control; SCLC: small-cell lung cancer
Differentially expressed lncRNAs and mRNAs
A total of 34,832 lncRNAs were detected. Among them, 320 lncRNAs were differentially expressed in the LEMS group compared to the control group. The number of differentially upregulated lncRNAs was 197, with ENST00000411554 being the most significantly upregulated at 3.67-fold. The number of differentially downregulated lncRNAs was 123, with ENST00000538934 being the most significantly downregulated (fold change = 4.29).
A total of 168 differentially expressed mRNAs were identified in the LEMS group compared to the control group. The number of differentially upregulated mRNAs was 132, with uc003tgl.2 being the most significantly upregulated (fold change = 3.92). The number of differentially downregulated mRNAs was 36, with the 10.69-fold reduction in NM_004345 being the most significant downregulation (Figure 1). Based on the values of -log10 (P-value) and -log2 (fold change) of differentially expressed lncRNAs and mRNAs, volcano maps were drawn (Figure 2).
Using joint analysis, we have found that 9 different lncRNA and their 11 neighboring mRNAs were differentially expressed in the LEMS group compared to the control group (Table 2). We considered mRNAs within 300 KB upstream or downstream of the lncRNA as their neighboring mRNAs. The results of joint analysis of differentially expressed lincRNAs and neighboring mRNAs are showed in Table 3.
Table 2 Analysis of differentially expressed lincRNAs and neighboring mRNAs between patients with LEMS and healthy controls
lncRNA
|
Gene symbol
|
Fold change
|
P-value
|
Regulation
|
Relative position
|
mRNA
|
Nearby gene symbol
|
T357796
|
G084427
|
2.02
|
0.049
|
Up
|
Downstream
|
NM_003289
|
TPM2
|
ENST00000427872
|
RP11-439A17.4
|
1.97
|
0.028
|
Up
|
Upstream
|
NM_001017986
|
FCGR1B
|
NR_024397
|
NUTM2A-AS1
|
1.76
|
0.041
|
Up
|
Upstream
|
NM_004897
|
MINPP1
|
ENST00000452525
|
AC022201.4
|
1.64
|
0.045
|
Up
|
Upstream
|
NM_022173
|
TIA1
|
NR_034139
|
LOC338963
|
1.63
|
0.029
|
Up
|
Downstream
|
NM_004644
|
AP3B2
|
NR_034139
|
LOC338963
|
1.63
|
0.029
|
Up
|
Upstream
|
NM_001007122
|
FSD2
|
ENST00000538934
|
RP11-807H22.6
|
4.29
|
0.032
|
Down
|
Downstream
|
NM_000802
|
FOLR1
|
ENST00000538934
|
RP11-807H22.6
|
4.29
|
0.032
|
Down
|
Upstream
|
ENST00000442948
|
FOLR3
|
ENST00000422305
|
RP11-399E6.1
|
1.76
|
0.0169
|
Down
|
Upstream
|
NM_001031694
|
SCMH1
|
T153963
|
G035671
|
1.62
|
0.045
|
Down
|
Downstream
|
NM_001080466
|
BTBD17
|
NR_120479
|
LOC101930452
|
1.51
|
0.050
|
Down
|
Downstream
|
NM_002258
|
KLRB1
|
Validation of microarray results by qRT-PCR
Among the mRNAs differentially expressed in LEMS patients, AP3B2 and the protein it synthesizes are known to be abnormally expressed in some nervous system diseases. Therefore, we focused on AP3B2 and its neighboring lncRNA LOC338963 (NR_034139) for further analysis. We used qRT-PCR to validate the microarray data using primers in the University of California, Santa Cruz database that were synthesized by Shanghai Xianjun Biotechnology Co., Ltd (Shanghai, China) (Table 3). LOC338963 (NR_034139) and AP3B2 were upregulated by 2.07- and 2.15-fold, respectively, in the LEMS group compared to the control group (Figure 3).
Table 3 PCR primers for LOC338963 (NR_034139) and AP3B2
Gene name
|
Primers
|
LOC338963
|
F: 5' GGCTACTCGGGGATAAAGACTG 3'
R: 5' GTCCGTCGCTTCGCTCACT 3'
|
AP3B2
|
F: 5' CACCTGTCGGAGTGACCACATT 3'
R: 5' GTGCCAATCACCATTTTCTCG 3'
|
b-actin (H)
|
F: 5' GTGGCCGAGGACTTTGATTG 3'
R: 5' CCTGTAACAACGCATCTCATATT 3'
|
F: forward; R: reverse
We further verified the differential expression of LOC338963 (NR_034139) and AP3B2 using qRT-PCR on blood samples taken from patients with either paraneoplastic LEMS (nine patients), non-tumor LEMS (eight), or SCLC (four) (Table 4). LOC338963 (NR_034139) and AP3B2 were significantly upregulated in the paraneoplastic LEMS group compared to both the non-tumor LEMS and SCLC groups (Figure 4). LOC338963 (NR_034139) and AP3B2 were expressed at 2.34- and 5.06-fold, and 4.36- and 14.97-fold higher levels in patients with paraneoplastic LEMS than in the comparative groups, respectively.
Table 4 Clinical characteristics of patients with P-LEMS, N-LEMS and SCLC
|
P-LEMS
|
N-LEMS
|
SCLC
|
Number of patients
|
9
|
8
|
4
|
Age (average ± SD)
|
65y±10
|
48y±17
|
63y±6
|
Sex (female/male)
|
5/4
|
4/4
|
1/3
|
Tumor Classification/number
|
SCLC/4
|
|
SCLC/4
|
|
liver cancer/3
|
|
|
|
carcinoma of stomach/1
|
|
|
|
thymic tumor/1
|
|
|
LEMS: Lambert-Eaton myasthenic syndrome; P-LEMS: paraneoplastic LEMS; N-LEMS: non-tumor LEMS; SCLC: small cell carcinoma of lung
GO and KEGG analysis
GO analysis of differentially expressed mRNAs revealed the following categories: immune system processes, tissue development, cellular defense response, natural killer cell-mediated immunity, and other functional pathways. Figures 5 and 6 show the top ten mRNAs according to the enrichment scores for Molecular Function (blue), Biological Process (red), and Cellular Component (green).
To infer the pathways in which the mRNAs participated, we performed KEGG enrichment analysis (Table 5 and Figure 7-8). The number of downregulated mRNAs was relatively small, so there were fewer enrichment pathways compared to those related to the upregulated mRNAs. KEGG analysis indicated that lncRNAs were mostly related to 12 biological information pathways, mainly focusing on antigen processing and presentation, natural killer cell-mediated cytotoxicity, longevity regulation, oxytocin signaling pathway, primary immune deficiency, glutathione synapse, salt resistance, and endocytosis.
GO analysis of AP3B2 and its interacting proteins showed that the main Biological Process enrichment terms were anterograde synaptic vesicle transport (Enrichment Score =15.80410035) and establishment of synaptic vesicle localization (Enrichment Score =12.86966623). The Cellular Component enrichment terms were AP-type membrane coat adaptor complex (Enrichment Score = 14.29073004) and axon cytoplasm (Enrichment Score = 14.19449914). The Molecular Function enrichment term was transporter activity (Enrichment Score = 1.468521083). KEGG analysis of AP3B2 and its interacting proteins revealed that AP3B2 was enriched into Lysosome pathways (Figure 9).
Table 5 Significantly different KEGG pathways in patients with LEMS versus healthy control
Pathways
|
Genes
|
P-value
|
Antigen processing and presentation
|
CD8B//KLRC1//KLRC2//KLRC3//KLRD1
|
0.00
|
Natural killer cell-mediated cytotoxicity
|
KLRC1//KLRC2//KLRC3//KLRD1//KLRK1
|
0.00
|
Malaria
|
CD40LG//KLRB1//KLRK1
|
0.00
|
Glutamatergic synapse
|
ADCY1//HOMER3//PLA2G4C//SLC1A7
|
0.01
|
Longevity regulatory pathway
|
ADCY1//IRS1//PRKAB2
|
0.01
|
Oxytocin signaling pathway
|
ADCY1//CD38//PLA2G4C//PRKAB2
|
0.01
|
Longevity regulatory pathway
|
ADCY1//IRS1//PRKAB2
|
0.02
|
Primary immunodeficiency
|
CD40LG//CD8B
|
0.02
|
Hematopoietic cell lineage
|
CD38//CD8B//FCGR1A
|
0.02
|
Graft-versus-host disease
|
KLRC1//KLRD1
|
0.03
|
Antifolate resistance
|
FOLR1//FOLR3//IZUMO1R
|
1.42
|
Endocytosis
|
FOLR1//FOLR3//IZUMO1R
|
0.01
|