Research Papers: Pathology:
Functional networks of aging markers in the glomeruli of IgA nephropathy: a new therapeutic opportunity
Metrics: PDF 2638 views | HTML 2539 views | ?
Hong Jiang1,2,3*, Ludan Liang1,2,3*, Jing Qin5,6, Yingying Lu1,2,3, Bingjue Li1,2,3, Yucheng Wang1,2,3, Chuan Lin1,2,3, Qin Zhou1,2,3, Shi Feng1,2,3, Shun H. Yip6, Feng Xu6, EnYin Lai7, Junwen Wang4,6 and Jianghua Chen1,2,3
1 Kidney Disease Center, The First Affiliated Hospital, College of Medicine, Zhejiang University, P.R. China
2 Kidney Disease Immunology Laboratory, The Third Grade Laboratory, State Administration Of Traditional Chinese Medicine Of China, Hangzhou, P.R. China
3 Key Laboratory Of Multiple Organ Transplantation, Ministry Of Health, Key Laboratory Of Nephropathy, Zhejiang, P.R. China
4 Centre for Genomic Sciences, LKS Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
5 School of Life Sciences, The Chinese University of Hong Kong, Hong Kong SAR, China
6 School of Biomedical Sciences, LKS Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
7 Department of Physiology, Zhejiang University School of Medicine, Hangzhou, China
* These authors have contributed equally to this work
Jianghua Chen, email:
Hong Jiang, email:
Junwen Wang, email:
Keywords: IgAN, glomerular, RNA sequence, aging, networks, Pathology Section
Received: December 12, 2015 Accepted: April 11, 2016 Published: April 26, 2016
IgA nephropathy(IgAN) is the most common primary glomerular disease in China. Primary infections always occur before IgAN. However, the pathology of IgAN is still unclear. Previously we found that LL37, a protein secreted by senescent cells, was specific for the progression of IgAN, and also played a role in the neutrophil function. So we hypothesized that the infiltration of neutrophils, inflammation factors, and aging markers , which were modulated by functional networks, induced the immune response and renal injury. RNA-Sequencing (RNA-seq) can be used to study the whole transcriptome and detect splicing variants that are expressed in a specific cell type or tissue. We separate glomerulus from the renal biopsy tissues. After RNA extraction, the sequences were analyzed with Illumina HiSeq 2000/2500. 381 genes with differential expression between the IgAN patients and the healthy controls were identified. Only PLAU, JUN, and FOS were related to DNA damage, telomere dysfunction-induced aging markers, neutrophil function and IgA nephropathy. The networks showed the possibility of these genes being connected. We conclude that DNA damage and telomere dysfunction could play important roles in IgA nephropathy. In addition, neutrophils are also important factors in this disease. The networks of these markers showed the mechanism pathways that are involved in the duration of the occurrence and progression of IgA nephropathy and might be a new therapeutic opportunity for disease treatment.
IgAN is the most common primary glomerular disease in China. The morbidity of IgA nephropathy is 40-50% in primary glomerular disease [1-3]. The characteristic pathological change is the proliferation of mesangial cells and the deposition of the IgA immune complex. There are always primary infections before IgAN [4, 5], but the detailed mechanism is still unclear. In addition, there is no effective treatment for IgA nephropathy. We have discovered a group of specific proteins that were secreted by senescent cells and were induced by telomere dysfunction . This group of proteins is a good predictor of not only human aging but also chronic disease, such as liver fibrosis and myelodysplastic syndrome (MDS). The follow-up studies also showed that among this group of proteins, LL37 is specific for the progression of IgAN . We found higher LL37 expression in renal tissues from IgAN patients, which is also confirmed by IF staining . Interestingly, cells with positive LL37 staining were considered to be neutrophils according to their nuclear shape, which led to the hypothesis that neutrophils may participate in the pathology of IgAN. LL37 not only is related to DNA damage and telomere dysfunction but also plays an important role in the function of the neutrophils . LL37 was very important in the SVV as a neutrophil extracellular traps(NETs) functional factor . LL37 played a novel role in the protection of NETs against degradation by bacterial nucleases . These data suggest some connective links between telomere dysfunction and NETs. These connective links might relate to the pathologic mechanism of IgAN. In this study, we hypothesized that the infiltration of neutrophils, inflammation factors, and aging markers, such as LL37, which were modulated by functional networks, induced the immune response and renal injury.With rapid technological progression, the elucidation of the transcriptome complexity of an organism and understanding the underlying functions of various differentially expressed genes have become a major focus of post-genome research. In addition to the classical microarray approach for profiling transcripts, the recent development of RNA-Sequencing (RNA-seq) has revolutionized the study of whole transcriptomes, providing a potentially unbiased measure of all transcripts and splicing variants that are expressed in a specific cell type [11,12] or tissue [11, 12]. However, the analysis and interpretation of the huge amount of data generated from RNA-Seq pose a practical challenge and demand accurate and easily automated bioinformatics tools for processing datasets . In this study, we detect differentially expressed RNA transcripts between healthy and IgAN patients in aging, inflammation, NETs, and IgAN databases and found the corresponding network, decreasing the cost for the mechanisms that are found, not only in fees, but also in labor. And we could present the bioinformatics proof for the following study of the mechanism’s verification.
Demographic characteristics of the healthy controls and the IgAN patients
Using RNA-Seq, we analyzed six human glomeruli, three are from the renal biopsy tissues of the transplant donors; these are the healthy controls with pathologic evidence. Other three are from IgAN patients. Table 1 shows the clinical characteristics of the healthy controls and the IgAN patients. The gender, age, urine protein concentration, urine RBC counts, serum ALB, SCr, serum BUN, eGFR (MDRD formula), SBP, DBP, and pathologic diagnosis are listed. IgA nephropathy patients showed higher proteinuria and urine RBC levels and lower eGFR levels. These diagnoses are supported by a pathologic analysis of kidney biopsy tissues.
Table 1: Clinical characteristics of 3 healthy controls and 3 IgA nephropathy patients.
Their gender, age, urine protein concentration, urine RBC counts, serum ALB, SCr, serum BUN, eGFR(MDRD formula), SBP, DBP, and pathologic diagnosis are listed above. IgA nephropathy patients showed more proteinuria, urine RBC, and worse eGFR. All those diagnosis were supported by pathologic analysis of kidney biopsy tissues.
Gene expression analysis by RNA-Seq
To investigate the genes that are changed in IgAN patients, we performed RNA-seq expression screen in the six renal glomeruli. Our bioinformatics analysis of the RNA-seq data revealed a total of 381 genes that were differentially expressed between the IgAN patients and the healthy controls (q-value < 0.05). Among them, 229 genes were upregulated, and 152 genes were down regulated in the IgA nephropathy patients compared to the healthy controls (Table S1, Suppl. Figure 1). Thus, our study focused on the hypothesis that modulation by functional networks was induced by the infiltration of aging neutrophils, which caused IgA immune complex deposition and renal injury. We next used the DAVID program to perform enrichment analysis between the differentially expressed genes and gene sets curated from Pubmed search and manual curation. We are particularly interested in genes involved in aging, inflammation, NETs and IgAN.
Aging transcript expression analysis by RNA-seq
A total of 21 aging transcripts showed differential expression between the IgAN patients and the healthy controls (q-value < 0.05). Of these transcripts, 15 were upregulated, and 6 were downregulated (Figure 1). Red and blue in the heat map depict higher and lower gene expression, respectively. The color intensity indicates the magnitude of the expression differences.
Figure 1: The glomerular expression levels of 21 transcripts correlated with aging. The red and blue colors in the heat map depict higher and lower gene expression, respectively. The color intensity indicates the magnitude of the expression differences. A fold change of 2 and a q-value of 0.05 were used as cutoffs. Official names according to the NCBI were given, test values of IgAN and the control are listed, the fold changes indicate the relative alteration of the control value to the IgAN value, and the log2 value of the fold change was used. A positive log2 (fold change) indicates upregulation in the disease group, and a negative log2 (fold change) indicates downregulation.
Inflammation transcript expression analysis by RNA-seq
A total of 25 inflammation transcripts showed differential expression between the IgAN patients and the healthy controls (q-value < 0.05). Of these transcripts, 15 were upregulated, and 10 were downregulated (Figure 2). Red and blue in the heat map depict higher and lower gene expression, respectively. The color intensity indicates the magnitude of the expression differences.
Figure 2: The glomerular expression levels of 25 transcripts correlated with inflammation. Refer to Figure 1 for description of the figure.
IgA nephropathy transcripts in a glomeruli analysis by RNA-seq
A total of 16 IgAN transcripts showed differential expression between the IgAN patients and the healthy controls (q-value < 0.05). Of these transcripts, 8 were upregulated, and 8 were downregulated (Figure 3). Red and blue in the heat map depict higher and lower gene expression, respectively. The color intensity indicates the magnitude of the expression differences.
Figure 3: The glomerular expression levels of 16 transcripts that were previously reported as being related to IgA nephropathy. Refer to Figure 1 for description of the figure.
NETs transcripts in a glomeruli analysis by RNA-seq
A total of 9 NETs transcripts showed differential expression between the IgAN patients and the healthy controls (q-value < 0.05). Of these transcripts, 2 were upregulated, while 7 were downregulated (Figure 4A). Red and blue in the heat map depict higher and lower gene expression, respectively. The color intensity indicates the magnitude of the expression differences.
We conducted NETs visualization with kidney biopsy tissues from IgA nephropathy patients and healthy donors (Figure 4B, 4C). No neutrophil extracellular traps (NETs) were deposited in the kidney biopsies from either the IgA nephropathy patients or the healthy donors. NETs were identified by co-localization with extracellular DNA (blue), the H2A-H2B-complex (green) and LL37 (red). Almost no co-localization of autoantigen was found in the kidney biopsy sections from both populations.
Figure 4: Expression levels of 16 transcripts in glomeruli showed correlation with IgAN. A. The glomerular expression levels of 9 transcripts correlated with NETs. A fold change of 2 and a q-value of 0.05 were used as cutoffs. Official names according to the NCBI were given, test values of IgAN and the control are listed, the fold change was used to indicate the relative alteration of the control value to the IgAN value, and the log2 value of the fold change was used. A positive log2 (fold change) indicates upregulation in the disease group, and a negative log2 (fold change) indicates downregulation. B. and C. No neutrophil extracellular traps (NETs) were deposited in kidney biopsies from either the IgA nephropathy patients or the healthy donors. NETs were identified by co-localization with extracellular DNA (blue), the H2A-H2B-complex (green) and LL37 (red). Almost no co-localization of autoantigens was found in the kidney biopsy sections from both populations. IgA: IgA nephropathy patients HD: healthy donars .
Overlap of aging, inflammation, and IgAN related genes
By overlap, we found three genes, JUN, FOS, and PLAU were simultaneously participated in aging and inflammation, and were also previously reported to be related to IgA nephropathy. Two of the three genes, JUN and FOS, could also increase the expression of LL-37 and are therefore indirectly related to NETs (Figure 5A). Figure 5B shows the network of differentially expressed genes related to aging, inflammation and IgAN. They participate in RAAS activation, complement and coagulation cascades.
Figure 5: Three transcripts of three genes, JUN, FOS, PLAU, highly expressed in IgAN patients, relatively, were simultaneously participated in aging and inflammation and IgAN. A. Three transcripts of the genes JUN, FOS, and PLAU were highly expressed in the IgAN patients and were simultaneously included in aging-, inflammation- and IgAN-related genes, respectively. Two of these genes, JUN and FOS could also increase the expression of LL-37. B. Network analysis includes differentially expressed genes related with aging(diamond), inflammation(circle) and IgAN(square). Red means up regulation and green means down regulation. Color depth indicates the log2(fold change) value. Lines in grey shows relationship between genes. Arrow means promotion and bar means inhibition. Genes in purple frame were key genes found in our study.
Putative mechanisms of aging neutrophil-related IgAN
Both JUN and FOS were components of the AP-1 transcription factor. The increased expression of AP-1 and NF-kappa B can both activate the RAAS system, resulting in IgA nephropathy. PLAU is a part in the complement and coagulation cascades, and therefore correlates with IgAN progression. Meanwhile, over-produced AIgA1 could increase JUN expression and thus promoted the disease. Matrix deposition as a characteristic of IgAN could also induce PLAU production in a protective role against disease (Figure 6).
Figure 6: Putative mechanism of of aging neutrophil-related IgAN. The RAAS system is activated by JUN, FOS, JUND(components of AP-1 transcription factor) and NF-kappa B, and subsequently induces IgAN. PLAU as a component of complement and coagulation cascades, can also promotes IgAN. Meanwhile, AIgA1 can increase JUN expression. Matrix deposition, an important phenomenon in IgAN can promote PLAU expression as a protective role.
Owning to different cell types and functions, genes differentially expressed between kidney interstitial, tubules and glomerulus. Because IgAN is a glomerular disease, the gene expression variance of the glomerulus should be more obvious and more related to the disease itself. The glomeruli were separated from the renal biopsy tissues, which allowed us to more precisely identify the differentially expressed transcripts specifically to the glomerulus of IgAN. The glomeruli separating techniques were stable and commonly used in mice researches [21, 22]. Because human samples are difficult to obtain, there are no studies on the human glomerulus. In addition, a human glomeruli dissection technique was recently developed. Koitabashi et al. set up a method to isolate human glomeruli from kidney biopsy tissues with further proteomic research  in 2011. Bruschi M et al. performed similar studies on LN patients in 2014. Most of these studies focus on method setup and protein inspection; RNA-seq of isolated human glomeruli is rare.
Aging is a biological process characterized by cell senescence and telomere dysfunction . When the cell is stimulated to grow, while the cell cycle is arrested, then the cell becomes senescent (geroconversion) . Irreversible cell cycle arrest happens in response to the telomere dysfunction . Our preliminary study found that LL37 was an aging biomarker induced by DNA damage and telomere dysfunction . LL37 is not only a marker related to IgAN  but also an important factor in NETs. Neutrophil infiltration is important in IgAN [27-30]. We hypothesized that the infiltration of the neutrophils, inflammation factors, and aging markers, such as LL37, which were modulated by functional networks, induced the immune response and renal injury. We analyzed the transcriptome and the differentially expressed genes and the genes are significantly enriched in process involved in aging, inflammation, and NETs. RNA-Sequencing techniques [31-33] provide a good way to verify this hypothesis at a low cost. Although previous study found LL37 closely related to IgAN, we did not observe differential RNA expression of LL37 between the healthy controls and the IgAN patients in this transcriptome, maybe because of limited sample size. However, some other aging markers were related to the IgAN patients. Among these 21 aging markers, 3 important aging markers were related to LL37: APOA1, JUN, and FOS. In addition, there was no difference in NETs between the IgAN patients and the healthy controls, suggesting that NETs were not specifically induced in IgAN patients. NETs were specific in some other immune diseases, such as SVV  but not IgAN.
PLAU, JUN, and FOS are also related to neutrophil infiltration [34-36]. Both JUN and FOS are components of the AP-1 transcription factor. JUND, though not included in the overlap group, is also a subunit of AP-1. The increased expression of AP-1 activates the RAAS system , subsequently resulting in IgA nephropathy. PLAU is a component of complement and coagulation cascades . And complement cascade is closely related with IgAN . Meanwhile, over-produced AIgA1 could increase JUN expression and thus promote disease onset . Matrix deposition, as one of the characteristics of IgAN, can also induce the production of PLAU as a protective role .
RNA-Sequence techniques have vastly helped in revising the initial hypothesis, saving money and labor. And we could present the bioinformatics proof for the following study of the mechanism’s verification.
There were several limitations in this study. First, because of the difficulty to obtain consent forms, the sample size is limited, only 6 samples were available for analysis with complete and qualified consent forms for a biopsy of kidney tissues. Second, the RNA from the glomerulus is only enough for RNA-Sequencing; with more RNA available, other analyses, such as MicroRNA sequencing, could be simultaneously performed, and the networks could be much more detailed. In this study, the functional networks were induced by the infiltration of aging neutrophils, which caused IgA immune complex deposition and renal injury. This study could help a lot in the follow-up study for detail mechanism’s verification. The networks that were correlated with these markers showed the mechanistic pathways that were involved in the occurrence and progression of IgA nephropathy, which might offer a new therapeutic opportunity for disease treatment.
MATERIALS AND METHODS
Studies on human samples were conducted according to the declaration of Helsinki. Three live-donor renal transplant donors supplied the renal biopsy tissues that were used in the transplantations; these tissues were the healthy controls with pathologic evidence. Three IgAN patients supplied the renal biopsy tissues. These renal tissues were collected from Jun. 2011 to Oct. 2011. The dissection procedures were previously published . The glomeruli were dissected at 4°C in an albumin-enriched physiological salt solution (0.1%) using a stereoscopic microscope and sharpened forceps. Tubuli were removed, with the exception of the region of the thick ascending limb of Henle’s loop, which touches the glomerulus. The dissection time was limited to 120 min after the biopsies, and the glomeruli were dissected under a microscope (Leica S8AP0).
RNA quantification and qualification
RNA degradation and contamination were monitored on 1% agarose gels.
RNA purities were checked using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA).
RNA concentrations were measured using the Qubit® RNA Assay Kit in the Qubit® 2.0 Fluorometer (Life Technologies, CA, USA).
RNA integrities were assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA).
Library preparation for sequencing
A total of 3 μg of RNA per sample was used as the input material for the RNA sample preparations. Sequencing libraries were generated using the NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA) following the manufacturer’s recommendations, and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from the total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under an elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5×). First-strand cDNA was synthesized using a random hexamer primer and M-MuLV Reverse Transcriptase (RNaseH-). Second-strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. The remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After the adenylation of the 3’ ends of the DNA fragments, the NEBNext Adaptor with a hairpin loop structure was ligated to prepare for hybridization. To select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with the AMPure XP system (Beckman Coulter, Beverly, USA). Then, 3 μl of USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37°C for 15 min followed by 5 min at 95 °C before PCR. Then, PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. Finally, the PCR products were purified (AMPure XP system), and the library quality was assessed on the Agilent Bioanalyzer 2100 system.
Clustering and sequencing (novogene experimental department)
The clustering of the index-coded samples was performed on a cBot Cluster Generation System using the TruSeq SR Cluster Kit v3-cBot-HS (Illumina) according to the manufacturer’s instructions. After the cluster generation, the library preparations were sequenced on an Illumina HiSeq 2000/2500 platform, and 150 bp/100 bp/50 bp paired/single-end reads were generated.
The quality of raw sequencing reads was checked with FastQC . Sequencing adapter sequences at the two ends of raw reads were cleaned by Trimmomatic . Then, the cleaned reads were mapped to the human genome hg19 by STAR . The gene expression was quantified, and the differentially expressed genes (DEGs) were identified by Cuffdiff 2  with the q-value < 0.05. Functional annotation and pathway enrichment analysis were performed with DAVID . Polysearch2 , web-server that supports text mining by integrating millions of scientific literatures from Medline and Pubmed as well as other databases to ensure a more accurate prediction, was used to search for inflammation/aging-related gene relationships. In the study, we used the gene as the query to find inflammation/aging related genes. In the advanced search options, the synonym list was default from the tool. For example, the filter words was customized using inflammation related keywords (inflammation, inflammatory, inflamed, pro-inflammatory, anti-inflammatory) to extract inflammation-related genes. A relevance score was calculated from the associations after the extraction, and later converted to Z score. Associated genes with z score above 1 were selected for subsequent analysis. IgA nephropathy related genes and proteins were searched through the PubMed database(Table S2).
Visualization of NETs
Double immunofluorescence staining was conducted on frozen sections of kidney biopsies from IgA nephropathy patients(n = 3) and healthy donors(n = 3). Slides were blocked and permeabilized with phosphate buffered saline containing 2% BSA (Sigma, Germany) and 0.2% Triton X-100 (Sigma, Germany) for 15 min. The incubation of primary antibodies, mouse monoclonal anti-DNA/Histone 1 (1:500; Millipore, America) and rabbit anti-LL37  was carried out overnight at 4 Celsius degree. Then, the samples were incubated with Alexa-Fluor-488-labeled goat-anti-mouse antibody (1:500; Invitrogen, UK) and Alexa-Fluor-633-labeled goat-anti-rabbit antibody (1:500; Invitrogen, UK). The mounted samples were viewed with an ICS SP8 confocal laser scanning microscope, and the intensity of NETs and LL37 was calculated using LAS-AF software.
This study was supported by the National Fund Committee of China (2012CB517603, 2012AA02A512, 81470938, 81200546, and 81300619) and 2014-2015 Dr. Cheng Yu Tung Fellowships.
ConflictS of Interest
There is no conflict of interest.
1. Cai GY and Chen XM. Immunoglobulin A nephropathy in China: progress and challenges. American journal of nephrology. 2009; 30:268-273.
2. Zhou FD, Zhao MH, Zou WZ, Liu G and Wang H. The changing spectrum of primary glomerular diseases within 15 years: a survey of 3331 patients in a single Chinese centre. Nephrology, dialysis, transplantation. 2009; 24:870-876.
3. Zhang L, Wang F, Wang L, Wang W, Liu B, Liu J, Chen M, He Q, Liao Y, Yu X, Chen N, Zhang J-e, Hu Z, Liu F, Hong D, Ma L, et al. Prevalence of chronic kidney disease in China: a cross-sectional survey. The Lancet. 2012; 379:815-822.
4. Wyatt RJ and Julian BA. IgA nephropathy. The New England journal of medicine. 2013; 368:2402-2414.
5. Zand L and Fervenza FC. Does tonsillectomy have a role in the treatment of patients with immunoglobulin A nephropathy? Nephrology, dialysis, transplantation. 2014; 29:1456-1459.
6. Jiang H, Schiffer E, Song Z, Wang J, Zurbig P, Thedieck K, Moes S, Bantel H, Saal N, Jantos J, Brecht M, Jeno P, Hall MN, Hager K, Manns MP, Hecker H, et al. Proteins induced by telomere dysfunction and DNA damage represent biomarkers of human aging and disease. Proceedings of the National Academy of Sciences of the United States of America. 2008; 105:11299-11304.
7. Lu YY, Yang X, Chen WQ, Ju ZY, Shou ZF, Jin J, Zhang XH, Chen JH and Jiang H. Proteins induced by telomere dysfunction are associated with human IgA nephropathy. Journal of Zhejiang University Science B. 2014; 15:566-574.
8. Brinkmann V, Reichard U, Goosmann C, Fauler B, Uhlemann Y, Weiss DS, Weinrauch Y, Zychlinsky A. Neutrophil Extracellular Traps Kill Bacteria. Science. 2004; 303:1532-5.
9. Kessenbrock K, Krumbholz M, Schonermarck U, Back W, Gross WL, Werb Z, Grone HJ, Brinkmann V and Jenne DE. Netting neutrophils in autoimmune small-vessel vasculitis. Nature medicine. 2009; 15:623-625.
10. Merres J, Hoss J, Albrecht LJ, Kress E, Soehnlein O, Jansen S, Pufe T, Tauber SC and Brandenburg LO. Role of the cathelicidin-related antimicrobial peptide in inflammation and mortality in a mouse model of bacterial meningitis. Journal of innate immunity. 2014; 6:205-218.
11. Mortazavi A, Williams BA, McCue K, Schaeffer L and Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature methods. 2008; 5:621-628.
12. Farkas MH, Au ED, Sousa ME and Pierce EA. RNA-Seq: Improving Our Understanding of Retinal Biology and Disease. Cold Spring Harbor perspectives in medicine. 2015; 5.
13. Steijger T, Abril JF, Engstrom PG, Kokocinski F, Hubbard TJ, Guigo R, Harrow J, Bertone P and Consortium R. Assessment of transcript reconstruction methods for RNA-seq. Nature methods. 2013; 10:1177-1184.
14. Andrews S. FastQC: a quality control tool for high throughput sequence data. Available online at: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/. 2010.
15. Bolger AM, Lohse M and Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30:2114-2120.
16. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M and Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013; 29:15-21.
17. Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL and Pachter L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013; 31:46-53.
18. Huang da W, Sherman BT and Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic acids research. 2009; 37:1-13.
19. Liu Y, Liang Y and Wishart D. PolySearch2: a significantly improved text-mining system for discovering associations between human diseases, genes, drugs, metabolites, toxins and more. Nucleic acids research. 2015; 43:W535-542.
20. Dorschner RA, Pestonjamasp VK, Tamakuwala S, Ohtake T, Rudisill J, Nizet V, Agerberth B, Gudmundsson GH, Gallo RL. Cutaneous Injury Induces the Release of Cathelicidin AntiMicrobial Peptides Active Against Group A Streptococcus. Journal of Investigative Dermatology. 2001; 117:91-97.
21. Lai E. Angiotensin II-Nitric Oxide interaction in glomerular arterioles. Clinical and Experimental Pharmacology and Physiology. 2005; 32:410-414.
22. Ilatovskaya DV and Staruschenko A. Single-channel analysis of TRPC channels in the podocytes of freshly isolated Glomeruli. Methods in molecular biology. 2013; 998:355-369.
23. Koitabashi K, Okamoto K, Arirto M, Sato T, Nagai K, Kurokawa MS, Suematsu N, Yasuda T, Kimura K and Kato T. Micro-sieving: isolation of whole glomeruli from a single renal needle biopsy sample. Nephron Clinical practice. 2011; 117:c225-229.
24. Lopez-Otin C, Blasco MA, Partridge L, Serrano M and Kroemer G. The hallmarks of aging. Cell. 2013; 153:1194-1217.
25. Blagosklonny MV. Cell cycle arrest is not yet senescence, which is not just cell cycle arrest: terminology for TOR-driven aging. Aging(Albany NY). 2012; 4:159-165. doi: 10.18632/aging.100443.
26. Alexandra Bernadotte VMM, Irina M. Spivak. Markers of cellular senescence. Telomere shortening as a marker of cellular senescence. Aging(Albany NY). 2016; 8:3-11. doi: 10.18632/aging.100871.
27. Lai HCCJYGJMCYH. Differential Effects of FMLP-Activated Neutrophils from Patients with IgA Nephropathy Enhanced Endothelin 1 Production of Glomerular Mesangial Cells. Nephron. 2001; 2001:274-279.
28. Imaizumi T, Aizawa T, Segawa C, Shimada M, Tsuruga K, Kawaguchi S, Matsumiya T, Yoshida H, Joh K and Tanaka H. Toll-like receptor 3 signaling contributes to the expression of a neutrophil chemoattractant, CXCL1 in human mesangial cells. Clinical and experimental nephrology. 2014.
29. Brown HJ, Lock HR, Wolfs TG, Buurman WA, Sacks SH and Robson MG. Toll-like receptor 4 ligation on intrinsic renal cells contributes to the induction of antibody-mediated glomerulonephritis via CXCL1 and CXCL2. Journal of the American Society of Nephrology. 2007; 18:1732-1739.
30. Sadik CD, Kim ND and Luster AD. Neutrophils cascading their way to inflammation. Trends in immunology. 2011; 32:452-460.
31. Shalek AK, Satija R, Shuga J, Trombetta JJ, Gennert D, Lu D, Chen P, Gertner RS, Gaublomme JT, Yosef N, Schwartz S, Fowler B, Weaver S, Wang J, Wang X, Ding R, et al. Single-cell RNA-seq reveals dynamic paracrine control of cellular variation. Nature. 2014; 510:363-369.
32. Ozsolak F and Milos PM. RNA sequencing: advances, challenges and opportunities. Nature reviews Genetics. 2011; 12:87-98.
33. Kolodziejczyk AA, Kim JK, Svensson V, Marioni JC and Teichmann SA. The technology and biology of single-cell RNA sequencing. Molecular cell. 2015; 58:610-620.
34. Cao W, Xu J, Zhou ZM, Wang GB, Hou FF and Nie J. Advanced oxidation protein products activate intrarenal renin-angiotensin system via a CD36-mediated, redox-dependent pathway. Antioxidants & redox signaling. 2013; 18:19-35.
35. De Borst MH, Prakash J, Melenhorst WB, van den Heuvel MC, Kok RJ, Navis G and van Goor H. Glomerular and tubular induction of the transcription factor c-Jun in human renal disease. The Journal of pathology. 2007; 213:219-228.
36. Maria Pia Rastaldi ST, Franco Ferrario, Antonio Indaco, Hequn Zou, Pietro Napodano and D’amico aG. Transforming growth factor-b, endothelin-1, and c-fos expression in necrotizing/crescentic IgA glomerulonephritis. Nephrology, dialysis, transplantation. 1998; (1998) 13:1668-1674.
37. Z.-G. SHEN J-LG, D.-S. LI. Screening of differentially expressed genes related to severe sepsis induced by multiple trauma with DNA microarray. European review for medical and pharmacological sciences. 2014; 2014; 18: 734-739.
38. Daha MR and van Kooten C. Role of complement in IgA nephropathy. Journal of nephrology. 2016; 29:1-4.
39. STEVEN C. DIVEN CRC, DIANNE K. HAMMOND, PAUL H. WEIGEL, JANET A. OKA, and GOLDBLUM aRM. IgA induced activation of human mesangial cells: Independent of FcR1 (CD 89). Kidney international. 1998; 54: 837-847.
40. Lee JH, Oh MH, Park JS, Na GJ, Gil HW, Yang JO, Lee EY and Hong SY. Urokinase, urokinase receptor, and plasminogen activator inhibitor-1 expression on podocytes in immunoglobulin A glomerulonephritis. The Korean journal of internal medicine. 2014; 29:176-182.
All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 4.0 License.