Identification of differentially expressed circulating exosomal lncRNAs in IgA nephropathy patients
BMC Immunology volume 21, Article number: 16 (2020)
Although immunoglobulin A nephropathy (IgAN) is one of the foremost primary glomerular disease, treatment of IgAN is still in infancy. Non-invasive biomarkers are urgently needed for IgAN diagnosis. We investigate the difference in expression profiles of exosomal long non-coding-RNAs (lncRNAs) in plasma from IgAN patients compared with their healthy first-degree relatives, which may reveal novel non-invasive IgAN biomarkers.
We isolated exosomes from the plasma of both IgAN patients and their healthy first-degree relatives. High-throughput RNA sequencing and real-time quantitative polymerase chain reaction (qRT-PCR) was used to validate lncRNA expression profiles. Pathway enrichment analysis was used to predict their nearest protein-coding genes.
lncRNA-G21551 was significantly down-regulated in IgAN patients. Interestingly, the nearest protein-coding gene of lncRNA-G21551 was found to be encoding the low affinity receptor of the Fc segment of immunoglobulin G (FCGR3B).
Exosomal lncRNA-G21551, with FCGR3B as the nearest protein-coding gene, was down-regulated in IgAN patients, indicating its potential to serve as a non-invasive biomarker for IgAN.
In Asian populations, Immunoglobulin A nephropathy (IgAN) is known as the most prevalent of primary chronic glomerular disease. IgAN is mostly characterized by the occurrence of IgA1 deposits in the mesangium of glomeruli; usually occurs in young or middle-aged adults. Among the patients, about 20 ~ 40% of cases eventually progress to end-stage renal disease (ESRD) within 10 ~ 20 years . Although the exact mechanism of IgAN is still largely unknown, it is thought to be an immune related disease with a multi-“hit” pathogenetic process with overproduction of aberrantly glycosylated IgA1(galactose-deficient IgA1, Gd-IgA1) . Till now, renal biopsy is still the standard for IgAN diagnosis, with no effective disease-specific therapies currently available. However, renal biopsy is invasive with limitations in assessing disease activity only at the time of biopsy, which could lead to inconclusive findings and decisions .
Although a number of IgAN-specific biomarkers have been discussed, the most reported is Gd-IgA1, it has not been validated in large multiracial cohorts of IgAN patients . There are still other studies reported on serum and urinary biomarkers, such as sCD89 , the transferrin receptor (TfR) , products of complement system (C3b, C3c, C5b-9) . However, there are limitations of these biomarkers, such as the validation in other populations of IgAN patients or lack sensitivity and specificity towards the disease diagnosis and progression. One study identified circulating sCD89-IgA immune complexes was thought about be associated with IgAN disease progression . However, another study measured sCD89-IgA immune complexes in the 326 IgAN patients and found no association with disease progression. With these contradicting results, it was suggested sCD89-IgA is not a good predictor for assessing the progression of IgAN . To date, there are no IgAN-specific biomarkers yet to replace renal biopsy as a diagnostic or to add valuable information in evaluating progression of IgAN revealing the urgent need to identify effective and non-invasive biomarkers to improve early detection and individualized treatment for IgAN.
Long non-coding RNAs (lncRNAs) are known as a heterogeneous class of transcripts with a of more than 200 bases, with no protein-coding potential . Emerging evidence shows that lncRNAs play important roles in various biological processes, including gene expression, protein folding recruitment during chromatin modifications, X-chromosome inactivation/ and immunoregulation . Recent studies have also reported the link between lncRNAs to various kidney diseases [10,11,12,13]. Li et al. reported that lncRNA MALAT1 expression was substantially increased in diabetic nephropathy (DN) and could reduce pyroptosis of renal tubular epithelial cells, antagonizing the role of miR-23c on the down-regulation of target gene ELAV like RNA binding protein 1 (ELAVL1), which resulted in a better understanding of the pathogenesis of DN and help the development of new therapeutic strategies . lncRNA X-inactive Specific Transcript (Xist) was identified to be up-regulated in membranous nephropathy (MN). Down-regulation of Xist may improve MN by reducing podocytes apoptosis via miR-217/TLR4 pathway . Liao et al. reported that lncRNA RP11-2B6.2 was increased in kidney tissue from lupus nephritis (LN) patients and positively associated with disease activity and the level of type-I interferon (IFN) . In IgAN patients, Zuo et al. reported 167 differentially expressed lncRNAs (55 up-regulated and 112 down-regulated lncRNAs) while compared with controls, which provided new information on the potential role of lncRNAs in IgAN . Subsequently, this explicit change of lncRNA profile indicated the possibility of lncRNAs to be used as a molecular biomarker for IgAN. However, the inherent instability of lncRNAs (easily degraded by RNase in the blood) results in utterly difficult and inconsistent detection of lncRNAs in the blood is, therefore limits their clinical application .
Exosomes are small membrane vesicles with diameters less than 150 nm that are present in nearly all biological fluids (e.g., blood, breast milk, saliva and urine) . Exosomes are known to encapsulate certain proteins, lipids, and RNA, and mediate intercellular communication between various types of cells . Increasing researches are showing that not only can exosomes be diagnostic and prognostic markers for various kinds of malignancies , but they also play an important role in immunoregulation and the pathogenesis of immune related diseases such as rheumatoid arthritis (RA) , systemic lupus erythematosus (SLE)  and IgAN . Exosomes could act as a protective barrier that can protect lncRNAs from extracellular degradation [14, 15], hence could be used as an excellent biomarker to detect significant changes in lncRNA profile in IgAN patients. In recent years, high-throughput RNA sequencing (RNA-seq) has become an attractive choice for the identification of differentially expressed genes (DEG) in various diseases as it has covers most of genome and exact detection of even those low-expressing genes. RNA-seq could achieve the resolution of a “single-base” and capture all transcripts . This study utilized RNA-seq to search for differentially-expressed blood exosomal lncRNAs between IgAN patients and their first-degree healthy relatives. In addition, real-time PCR (qRT-PCR) was also conducted to validate the results of RNA-seq, which may provide clues to identify potential novel lncRNA biomarkers for IgAN.
In this study, we recruited 17 IgAN patients and their healthy controls (their healthy first-degree relatives). The study was divided into two phases: (i) the screening phase was composed of 6 patients and 6 healthy controls; and (ii) the validation phase, with 11 patients and 11 healthy controls. The demographic and clinical data of subjects were summarized in Table 1. There were no significant differences in age, sex, BMI and blood pressure between patients with IgAN and the controls. However, IgAN patients presented higher 24-UPE (1.14 ± 0.51 vs. 0.00 g/24 h), sCr (94 ± 34.02 vs. 67.00 ± 19.42 μmol/L, p = 0.007), UA (422.0 ± 118.0 vs. 341.7 ± 86.16 μmol/L, p = 0.030), BUN (6.01 ± 2.90 vs. 4.64 ± 1.31 mmol/L, p = 0.086) and lower eGFR (84.66 ± 36.76 vs. 125.89 ± 53.90 ml/min/1.73m2,p = 0.014) as compared to healthy controls.
Identification and characterization of plasma exosomes
Exosome Precipitation Solution (ExoQuick-TC, System Biosciences) was used to isolate exosomes from plasma. After exosome isolation and purification, ZETASIZER Nano series-Nano-ZS (Malvern Instruments Ltd., Malvern, UK) was used to determine the hydrodynamic size of the exosomes. As shown in Fig. 1a, the size of exosome was ~ 30–200 nm. Flow cytometry analysis detected high exosomal surface marker proteins CD63 (87.6%) and CD81 (96.6%), respectively, confirming the purity of exosomes (Fig. 1b). TEM was also performed to determine the physical morphology and size of the exosomes, as seen in Fig. 1c, TEM image showed the lipid bilayer membrane at approximately 100 nm.
Identification of differentially-expressed lncRNA profiles
RNA-Seq (Illumina, Sna Diego, TX, USA) was performed to characterize the lncRNA expression profiles of exosomes in plasma samples of patients with IgAN and healthy relatives. As shown in Table 2, 70 lncRNAs were differentially expressed with significant fold-change (|log2(FC)| > 1) and base mean values. Among the 70 lncRNAs, 31 lncRNAs were upregulated and 39 lncRNAs were down-regulated in the IgAN group compared to control. The heatmap of differential expression and hierarchical clustering of lncRNAs in plasma samples of patients with IgAN and their corresponding relatives was demonstrated in Fig. 2a, while the volcano plot of differential expression of lncRNAs was showed in Fig. 2b.
Validation of candidate lncRNAs by qRT-PCR and nearest protein-coding genes
To confirm the results obtained from high-throughput sequencing, ten candidate lncRNAs were selected from both up-regulated and down-regulated groups according to their base mean and |log2(FC)| (Table 2). qRT-PCR was performed on validation cohort (11 IgAN patients with their healthy relatives). Consistent with the sequencing results, lncRNA-G21551 was significantly down-regulated in patients with IgAN compared with their healthy relatives (8.30 ± 1.32 vs 15.90 ± 3.18 for IgAN patients vs control, p = 0.045) (Fig. 3a). Two other candidates, lnc-SPATA31E1–10 and lncRNA-G111779, were also down-regulated in patients with IgAN, but the difference was not statistically significant (6.11 ± 0.89 vs 14.12 ± 3.64, p = 0.055, 6.00 ± 0.83 vs 10.07 ± 1.86, p = 0.066 for IgAN patients vs control respectively) (Fig. 3b, c). The remaining 17 lncRNAs could not be validated by qRT-PCR due to their low abundance. FCGR3B was calculated to be the nearest protein-coding gene of lncRNA-G21551 using pathway enrichment analysis and visualization using the R package clusterProfiler. FCGR3B encodes for the low affinity receptor (FcgR3B receptor) of the Fc segment of immunoglobulin G (IgG) .
In this study, the exosomal lncRNA profiles of IgAN patients were measured and compared to their healthy first-degree relative. Through high-throughput RNA seq, we identified lncRNA-G21551 as a potential diagnostic biomarker for IgAN. We also predicted its potential role in IgAN pathogenesis through pathway enrichment analysis and visualization using R package clusterProfiler. The nearest protein-coding gene of lncRNA-G21551 was identified to be FCGR3B, which encodes for the low affinity receptor (FcgR3B receptor) of the Fc segment of immunoglobulin G (IgG). Until the pathogenesis of IgAN is elucidated, renal biopsy will remain the golden standard for the diagnosis of IgAN. Therefore, our finding may bring forward a clue to find the disease biomarkers by exosomal lncRNA profiles in IgAN. The establishment of differentially expressed exosomal lncRNA profiles in IgAN could be important for illustrating the pathogenesis of IgAN.
lncRNAs-based biomarkers have been reported in a variety of diseases including membranous nephropathy, IgAN and hepatocellular carcinoma [11, 13, 14]. Zuo et al. identified 167 differentially expressed lncRNAs (including 55 upregulated lncRNAs and 112 downregulated lncRNAs) in peripheral blood mononuclear cells (PBMCs) of IgAN, which may aid in the elucidation of a basic pathogenic mechanism . However, unprotected lncRNAs are easily degraded by RNAse, and exosomes act as a protective layer that can prevent extracellular degradation of lncRNAs . Recently, emerging studies on exosomal non-coding RNAs (including lncRNAs and miRNAs) in renal disease have been reported , but few reports has been focused on IgAN. Min et al. reported a significantly difference in urinary exosomal miRNA profiles (including miR-29c, miR-146a and miR-205) between IgAN and healthy controls, which may serve as novel biomarkers for IgAN . However, the differential expression of lncRNAs in plasma exosomes in IgAN patients was not reported, which may contribute to the discovery of potential new biomarkers or pathogenic factors of IgAN from plasma, as the increased levels of circulatory polymeric IgA1 with aberrant O-glycosylation of its hinge region was reported to be the first-hit as of this disease .
In order to rule out the influence of genetic background, healthy first-degree relatives of the patients were used as normal controls. Through RNA-seq and qRT-PCR analysis, a large number of exosomal lncRNAs in the plasma were found differentially expressed between IgAN patients and their healthy first-degree relative, among them, the expression of lncRNA-G21551 were found to be significantly down-regulated in IgAN patients (8.298 ± 1.319 vs 15.896 ± 3.176 for IgAN patients vs control, p < 0.05), therefore, lncRNA-G21551 may serve as a biomarker for IgAN. However, the remaining 17 lncRNAs could not be validated by qRT-PCR due to their low abundance in exosomes. In our study, we used pathway enrichment analysis and visualization using the R package clusterProfiler and found that the nearest protein-coding gene of lncRNA-G21551 was FCGR3B, and hypothesize that lncRNA-G21551 may play a vital role in the pathogenesis of IgAN by regulating the expression of FCGR3B.
It was previously reported that copy number (CN) variation of the FCGR3B gene is associated with susceptibility to systemic lupus nephritis (SLE) and ANCA-associated systemic vasculitis (AASV) [23, 24]. Furthermore, the FcgR3B receptor is primarily expressed on neutrophils; while FCGR3B CNs are correlated with the expression of FcgR3B, functioning for the clearance of immune complex . IgAN is frequently characterized by depositions of IgA (mainly IgA1) or IgA-containing immune complexes in the glomerular mesangial areas or the capillary wall. IgA1 deposits are usually detected along with complement component 3 (C3), and often with IgG or IgM or both in glomeruli . However, a recent study has shown that FCGR3B polymorphisms have significant influence on the incidence and pathological grade of IgAN, suggesting that the impairment of IgG-IC clearance by the FCGR3B gene and subsequent glomerular deposition may also contribute to the glomerular lesions .
However, in this current study, we did not carry out mechanistic study on the direct evidence that could elucidate the interaction between lncRNA-G21551 and FCGR3B and the mechanisms involved in IgA nephropathy. The functional research of lncRNA-G21551 in IgAN may be an interesting new research area and is currently the scope of our research group.
In summary, our study demonstrates a significant difference in plasma exosomal lncRNA expression profiles between IgAN patients and their first-degree relatives, providing novel information on the potential role of exosomal lncRNAs in IgAN. Pathway enrichment analysis and visualization reveals that the FCGR3B gene may be closely associated with the pathogenesis of IgAN. Therefore, the expression level of exosomal lncRNA-G21551 could be utilized as a promising biomarker for IgAN diagnosis.
Participants and sample selection
Patients with primary IgAN from Sun Yat-sen Memorial Hospital, Sun Yat-sen University (Guangzhou, China) between March 2018 and March 2019, and their first-degree relatives were recruited in this study. Inclusion criteria are as follows: 1) biopsy-proven IgAN (within 30 days prior to enrollment); 2) age ≥ 14 years; and 3) an adequate biopsy sample with ≥10 glomeruli. Exclusion criteria are as follows: 1) secondary IgAN (Henoch-Schönlein purpura, systemic lupus erythematosus, liver disease, etc.), 2) an estimated glomerular filtration rate (eGFR) < 30 mL/min/1.73 m2 (calculated by the Chronic Kidney Disease Epidemiology Collaboration [CKD-EPI] creatinine equation ; 3) prior treatment with RAAS (renin-angiotensin-aldosterone system) inhibitor and / or immunosuppressants drugs; 4) presence of diabetes, concomitant infections, severe metabolic syndrome, and malignant tumors. As this is a patient-control matched study, patients’ healthy first-degree relatives (parents, siblings, or children) were chosen as respective controls. For this study, a total of 17 patients and their first-degree relatives were recruited. This study was conducted according to the principles of the Declaration of Helsinki and was approved by the Ethical Review Committee of Sun Yat-sen Memorial Hospital, Sun Yat-sen University (SYSEC-KY-KS-2018-080). Written informed consent was obtained from all participants prior to the study.
Whole blood samples were collected from each participant using anticoagulant EDTA tubes and centrifuged at 3000×g, 4 °C for 10 min. Then, the supernatant was centrifuged at 15000×g under the same conditions. The plasma supernatant was stored immediately at − 80 °C until further analysis. Samples for RNA-Seq were obtained from 6 patients and 6 relatives. For validation, samples from other 11 patients with IgAN and 11 relatives were subjected to qRT-PCR to detect the expression level of ten candidate lncRNAs selected from both up- and down-regulated groups according to their base mean and |log2(FC)|. Demographic and baseline clinical data including gender, age, 24-h urinary protein excretion (UPE) and serum creatinine (sCr), blood urea nitrogen (BUN), Serum uric acid (UA) and eGFR were recorded at the time of kidney biopsy.
Exosome isolation and identification
Exosomes form plasma of IgAN patients and their first-degree relatives were isolated using Exosome Precipitation Solution (ExoQuick-TC, System Biosciences, USA) according to the manufacturer’s instructions. Morphology of the isolated exosomes was then identified with Transmission Electron Microscope (TEM), size distribution analysis. Flow cytometer analysis was then use to confirm the purity of isolated exosomes.
After exosome extraction from serum, total RNAs were extracted from the17 patients and the corresponding relatives using miRNeasy Mini kit (Qiagen, Germany) individually according to the manufacturer’s instructions. Quantification of the total RNAs was performed by the Agilent 2200 TapeStation (Agilent Technologies, CA, USA).
cDNA library construction and high-throughput RNA sequencing
For RNA-seq analysis, total RNA from the exosomes of the 6 IgAN patients and their corresponding relatives was used for library preparation and sequencing, which were performed at RiboBio (Guangzhou, China). Briefly, RNA was fragmented to approximately 200 bp. The individual RNA sample were then subjected to first and second strand cDNA synthesis followed by adaptor ligation and low-cycle PCR enrichment according to the instructions provided with the NEBNext® Ultra™ RNA Library Prep Kit for Illumina (NEB, USA). The purified library products were evaluated using the Agilent 2200 TapeStation and Qubit®2.0 (Life Technologies, USA) and then sequenced (2× 150 bp) using a HiSeq30000.
Sequencing data analysis
To obtain high-quality, clean sequencing data, Fastp (version 0.19.4)  was used to filter low-quality reads, to cut adapters and for quality control of raw FASTQ files to obtain clean reads. The clean reads of each experiment were aligned against the human genome (UCSC/hg19) with HISAT2 (version 2.1.0)  and then subsequently assembled by StringTie (version 1.3.4d)  separately. All assemblies were merged into one transcriptome by TACO .
The newly assembled transcriptome was aligned to GENCODE v27 and Lncipedia v5.2 using GFFCompare (http://github.com/gpertea/gffcompare, version 0.10.1) to find novel transcripts which were assigned “class code” values of ‘i’, ‘u’ or ‘x’. The distances between novel transcripts and reference protein-coding transcripts were calculated by BED Tools (version 1.2.4) . CPAT  and PLEK (version 1.2)  were used to calculate the coding potential of novel transcripts. Salmon (version 1.11.2)  was applied to quantify transcript expression.
To analyze differential gene expression of lncRNA, several correlative packages in R were used. Tximport (version 1.12.3)  was applied to import quantification of transcript expression in R. Then, differentially expressed genes were determined by DESeq2 (version 1.24.0)  using FDR 0.05 as the threshold, and ggplot2  was used for visualization. Nearest protein-coding genes for differential lncRNAs were used to perform pathway enrichment analysis and visualization using the R package clusterProfiler .
Quantitative real-time PCR (qRT-PCR) analysis
qRT-PCR was used to verify the RNA-Seq data. LncRNAs were chosen based on expression level and biological significance. Sixteen μL of total RNA was used to synthesize the first strand of cDNA using PrimeScript™ RT Master Mix (Catalog No. RR036A, Takara, Japan). Real-time PCR was performed using TB Green (Catalog No. RR420A, Takara, Japan) in 96-well plates using the Biorad CFX384 Real-Time System (Bio-Rad, CA). The relative levels of target exosome-packaged lncRNAs were normalized against a synthesized exogenous reference λ polyA+ RNA (Catalog No. 3789, Takara, Japan).
Statistical analysis was performed with IBM SPSS Statistics 22.0 software (SPSS Inc., Chicago, IL, USA). Continuous data are expressed as mean ± standard deviation (S.D.). Data conforming to normal distribution were compared using Student t-test, while those with non-normally distributed were tested using Mann-Whitney U-test. Percentages (%) or frequencies was used for categorical data, and chi-squared test was used for comparison analysis between groups. p < 0.05 was considered statistically significant.
Availability of data and materials
Neither this manuscript nor substantial parts of it are under consideration for publication elsewhere, have been published nor made available elsewhere in a manner that could be construed as a prior or duplicate publication of the same content. The predicted sequence of lncRNA-G21551 was shown in supplementary data.
Immunoglobulin A nephropathy
Fc-Gamma Receptor 3B
End-stage renal disease
Systemic lupus erythematosus
ELAV like RNA binding protein 1
X-inactive Specific Transcrip
Differentially expressed genes
Real-time quantitative polymerase chain reaction
Estimated glomerular filtration rate
Chronic Kidney Disease Epidemiology Collaboration
Urinary protein excretion
Blood urea nitrogen
Transmission Electron Microscope
Peripheral blood mononuclear cell
ANCA-associated systemic vasculitis
- BED Tools:
Browser Extensible Data Tools
Coding-Potential Assessment Tool
Predictor of long non-coding RNAs and messenger RNAs based on an improved k-mer scheme
Transcriptome Assemblies Combined into One
Wyatt RJ, Julian BA. IgA nephropathy. N Engl J Med. 2013;368:2402–14.
Moresco RN, Speeckaert MM, Delanghe JR. Diagnosis and monitoring of IgA nephropathy: the role of biomarkers as an alternative to renal biopsy. Autoimmun Rev. 2015;14:847–53.
Schena FP, Cox SN. Biomarkers and precision medicine in IgA nephropathy. Semin Nephrol. 2018;38:521–30.
Moresco RN, Speeckaert MM, Zmonarski SC, et al. Urinary myeloid IgA fc alpha receptor (CD89) and transglutaminase-2 as new biomarkers for active IgA nephropathy and henoch-Schonlein purpura nephritis. BBA Clin. 2016;5:79–84.
Delanghe SE, Speeckaert MM, Segers H, et al. Soluble transferrin receptor in urine, a new biomarker for IgA nephropathy and Henoch-Schonlein purpura nephritis. Clin Biochem. 2013;46:591–7.
Endo M, Ohi H, Ohsawa I, et al. Glomerular deposition of mannose-binding lectin (MBL) indicates a novel mechanism of complement activation in IgA nephropathy. Nephrol Dial Transplant. 1998;13:1984–90.
van der Boog PJ, De Fijter JW, Van Kooten C, et al. Complexes of IgA with FcalphaRI/CD89 are not specific for primary IgA nephropathy. Kidney Int. 2003;63:514–21.
Jhee JH, Kang HY, Wu M, et al. Circulating CD89-IgA complex does not predict deterioration of kidney function in Korean patients with IgA nephropathy. Clin Chem Lab Med. 2017;56:75–85.
Quinn JJ, Chang HY. Unique features of long non-coding RNA biogenesis and function. Nat Rev Genet. 2016;17:47–62.
Li X, Zeng L, Cao C, et al. Long noncoding RNA MALAT1 regulates renal tubular epithelial pyroptosis by modulated miR-23c targeting of ELAVL1 in diabetic nephropathy. Exp Cell Res. 2017;350:327–35.
Huang YS, Hsieh HY, Shih HM, et al. Urinary Xist is a potential biomarker for membranous nephropathy. Biochem Biophys Res Commun. 2014;452:415–21.
Liao Z, Ye Z, Xue Z, et al. Identification of renal long non-coding RNA RP11-2B6.2 as a positive regulator of type I interferon signaling pathway in lupus nephritis. Front Immunol. 2019;10(975):1–11.
Zuo N, Li Y, Liu N, et al. Differentially expressed long noncoding RNAs and mRNAs in patients with IgA nephropathy. Mol Med Rep. 2017;16:7724–30.
Lee YR, Kim G, Tak WY, et al. Circulating exosomal noncoding RNAs as prognostic biomarkers in human hepatocellular carcinoma. Int J Cancer. 2019;144:1444–52.
Shah R, Patel T, Freedman JE. Circulating extracellular vesicles in human disease. N Engl J Med. 2018;379:958–66.
Wang T, Nasser MI, Shen J, et al. Functions of Exosomes in the triangular relationship between the tumor, inflammation, and immunity in the tumor microenvironment. J Immunol Res. 2019;2019:4197829.
Yoo J, Lee SK, Lim M, Sheen D, Choi EH, Kim SA. Exosomal amyloid a and lymphatic vessel endothelial hyaluronic acid receptor-1 proteins are associated with disease activity in rheumatoid arthritis. Arthritis Res Ther. 2017;19:119.
Sole C, Cortes-Hernandez J, Felip ML, et al. miR-29c in urinary exosomes as predictor of early renal fibrosis in lupus nephritis. Nephrol Dial Transplant. 2015;30:1488–96.
Min QH, Chen XM, Zou YQ, et al. Differential expression of urinary exosomal microRNAs in IgA nephropathy. J Clin Lab Anal. 2018;32:1–9.
Marguerat S, Bahler J. RNA-seq: from technology to biology. Cell Mol Life Sci. 2010;67:569–79.
Takai T. Roles of fc receptors in autoimmunity. Nat Rev Immunol. 2002;2:580–92.
Karpman D, Stahl AL, Arvidsson I. Extracellular vesicles in renal disease. Nat Rev Nephrol. 2017;13:545–62.
Fanciulli M, Norsworthy PJ, Petretto E, et al. FCGR3B copy number variation is associated with susceptibility to systemic, but not organ-specific, autoimmunity. Nat Genet. 2007;39:721–3.
Qi Y, Zhou X, Bu D, et al. Low copy numbers of FCGR3A and FCGR3B associated with Chinese patients with SLE and AASV. Lupus. 2017;26:1383–9.
Davies KA, Robson MG, Peters AM, et al. Defective fc-dependent processing of immune complexes in patients with systemic lupus erythematosus. Arthritis Rheum. 2002;46:1028–38.
Xu G, He Q, Shou Z, et al. NA1/NA2 heterozygote of Fcgr3b is a risk factor for progression of IgA nephropathy in Chinese. J Clin Lab Anal. 2007;21:298–302.
Levey AS, Stevens LA, Schmid CH, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150:604–12.
Chen S, Zhou Y, Chen Y, et al. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:i884–90.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.
Pertea M, Pertea GM, Antonescu CM, et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33:290–5.
Niknafs YS, Pandian B, Iyer HK, et al. TACO produces robust multisample transcriptome assemblies from RNA-seq. Nat Methods. 2017;14:68–70.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.
Wang L, Park HJ, Dasari S, et al. CPAT: coding-potential assessment tool using an alignment-free logistic regression model. Nucleic Acids Res. 2013;41:e74.
Li A, Zhang J, Zhou Z. PLEK: a tool for predicting long non-coding RNAs and messenger RNAs based on an improved k-mer scheme. Bmc Bioinformatics. 2014;15:311.
Patro R, Duggal G, Love MI, et al. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14:417–9.
Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Res. 2015;4:1521.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Ginestet C. Review Reviewed Work(s): ggplot2: Elegant Graphics for Data Analysis by H. Wickham Review by: Cedric Ginestet Source: Journal of the Royal Statistical Society. Series A (Statistics in Society), Vol. 174, No. 1 (JANUARY 2011), p. 245 Published by: Wiley for the Royal Statistical Society Stable URL: https://www.jstor.org/stable/23013414. Accessed: 20-06-2019 10:08 UTC. Journal of the Royal Statistical Society (2011) 1: 245.
Yu G, Wang LG, Han Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16:284–7.
We thank Huayin Health Technology Co., Ltd. (Guangzhou, China) for providing us with the technical support of exosomes identification with Transmission Electron Microscope (TEM). We also thank all patients who participated in this study.
This study was supported by grants of the National Natural Science Foundation of China (No. 81770709, 81370786), Guangzhou Municipal Science and Technology Project (No. 201904010142), and Guangdong Science and Technology Department (No. 2017B030314026).
Ethics approval and consent to participate
This study was conducted according to the Declaration of Helsinki and was approved by the Ethical Review Committee of Sun Yat-sen Memorial Hospital, Sun Yat-sen University (SYSEC-KY-KS-2018-080). Written informed consent was obtained from all participants prior to the study.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Guo, N., Zhou, Q., Huang, X. et al. Identification of differentially expressed circulating exosomal lncRNAs in IgA nephropathy patients. BMC Immunol 21, 16 (2020). https://doi.org/10.1186/s12865-020-00344-1