- Methodology article
- Open Access
Recommendations for mRNA analysis of micro-dissected glomerular tufts from paraffin-embedded human kidney biopsy samples
BMC Molecular Biology volume 19, Article number: 2 (2018)
Glomeruli are excellent pre-determined natural structures for laser micro-dissection. Compartment-specific glomerular gene expression analysis of formalin-fixed paraffin-embedded renal biopsies could improve research applications. The major challenge for such studies is to obtain good-quality RNA from small amounts of starting material, as applicable for the analysis of glomerular compartments. In this work, we provide data and recommendations for an optimized workflow of glomerular mRNA analysis.
With a proper resolution of the camera and screen provided by the next generation of micro-dissection systems, we are able to separate parietal epithelial cells from glomerular tufts. Selected compartment-specific transcripts (WT1 and GLEPP1 for glomerular tuft as well as PAX2 for parietal epithelial cells) seem to be reliable discriminators for these micro-dissected glomerular substructures. Using the phenol–chloroform extraction and hemalaun-stained sections (2 µm), high amounts of Bowman’s capsule transections (> 300) reveal sufficient RNA concentrations (> 300 ng mRNA) for further analysis. For comparison, in unstained sections from a number of 60 glomerular transections upwards, a minimum amount of 157 ng mRNA with a reasonable mRNA purity [A260/A280 ratio of 1.5 (1.4/1.7) median (25th/75th percentiles)] was reversely transcribed into cDNA. Comparing the effect of input RNA (20, 60, 150 and 300 micro-dissected glomerular transections), transcript expression of POLR2A significantly correlated when 60 and 150 laser micro-dissected glomerular transections were used for analysis. There was a lower inter-assay coefficient of variability for ADAMTS13, when at least 60 glomerular transections were used. According to the algorithms of geNormPlus and NormFinder, PGK1 and PPIA are more stable glomerular reference transcripts compared to GUSB, GAPDH, POLR2A, RPLPO, TBP, B2M, ACTB, 18SrRNA and HMBS.
Our approach implements compartment-specific glomerular mRNA expression analysis into research applications, even regarding glomerular substructures like parietal epithelial cells. We recommend using of at least 60 micro-dissected unstained glomerular or 300 hemalaun-stained Bowman’s capsule transections to obtain sufficient input mRNA for reproducible results. Hereby, the range of RNA concentrations in 60 micro-dissected glomeruli is low and appropriate normalization of Cq values using our suggested reference transcripts (PGK1 and PPIA) allows compensation with respect to different amounts of RNA purity and quantity.
Laser capture micro-dissection (LCM) enables the isolation of single cells or small cell groups from histological sections under direct microscopic control . Combined with quantitative PCR, LCM is a very powerful approach for studying individual transcript expression profiles in discrete cell populations . Especially for renal biopsies, micro-dissection of its different compartments—glomeruli, tubulointerstitium and arteries/arterioles—is a worthwhile exercise despite its small amounts of degraded RNA. About 17 years ago, glomerular micro-dissection was established in frozen sections obtained from a rat model . The method itself [4,5,6] and its potential role in diagnostic work-up of renal biopsies [7, 8] were reviewed several years ago. Hereby, general remarks for the major challenge for such approaches to obtain RNA with sufficient quantity and proper quality from small amounts of starting material are given, but no data for experimental procedures in detail. Especially, it has been shown that the majority of studies were performed with frozen tissue. This limited experience with formalin-fixed paraffin-embedded (FFPE) studies is due to the fact that frozen tissue is preferred, which is most often available from animal studies. However, often this is hardly an option for the large-scale analysis of human FFPE tissue from renal biopsies.
A comprehensive study for RNA quality and quantity was recent performed for the first time in renal FFPE tissue . Furthermore, six studies [9,10,11,12,13] including ours  have performed successful mRNA gene expression analysis from micro-dissected FFPE glomerular transections using different methodological approaches. Another study described a selection of reference transcripts for mRNA quantification from these specimen . However, there was no clear-cut recommendation for glomerular reference mRNA transcripts and renal biopsies stored in RNAlater . Most often, GAPDH (glyceraldehyde-3-phosphate dehydrogenase) is recommended as a glomerular reference transcript, although there is still no systematic study upon which the preference for GAPDH can be based .
Therefore, one aim of this study was to identify stably-expressed reference transcripts by proper normalization strategies and compare our approach for quantitative mRNA analysis of micro-dissected glomerular compartments with others. Another focus was to illustrate the possibilities for laser micro-dissection of glomerular substructures in detail. In addition, we studied the effect of the amount of micro-dissected glomerular transections for further RNA purity and quantity control, especially the net result for relative quantification. Hereby, we focused on two candidate reference transcripts (GAPDH and POLR2A [polymerase (RNA) II (DNA directed) polypeptide A)] and one target transcript [ADAMTS13 (a disintegrin and metalloproteinase with thrombospondin type 1 motives 13)] due to our interest in thrombotic microangiopathy.
Patients and biopsies
Overall, FFPE native and transplant kidney biopsies from 24 patients were selected from the archive of the Institute of Pathology, Hannover Medical School and Institute of Pathology, Department of Nephropathology, Friedrich-Alexander University Erlangen-Nuremberg. (i) Nine native kidney biopsies with acute tubular injury (ATN) or interstitial inflammation without any glomerular disease were selected for the practicability of the micro-dissection of parietal epithelial cells (PECs, patient cohort I). (ii) In order to determine reference transcripts, we used an additional fifteen patients selected by a high number (> 30) of glomeruli per biopsy, including a wide range of native and transplant biopsies (patient cohort II). These patients had different glomerular, tubulo-interstitial and hypertension-related renal diseases. Eight of these patients were analyzed for the effect of input RNA. Patient characteristics for all biopsies are detailed in Table 1. Clinical data included blood pressure, proteinuria (none, non-nephrotic, nephrotic), hematuria and serum-creatinine. There were significantly more hypertensive patients in cohort II compared to cohort I (p = 0.0234; Mann–Whitney U test). We could not find any significant differences for other clinical data between different patient groups.
Micro-dissection and mRNA analysis
All biopsies were cut at a thickness of 2 µm. In order to test the practicability of micro-dissection of PECs (patient cohort I), biopsies were deparaffinized for 30 s in xylol, stained for 2 s in hemalaun (Mayers Hemalaun, Merck, Germany, Cataloge number 109249) and washed in sterile aqua ad injectabilia (Ampuwa, Plastipur, Fresenius, Germany, Cataloge number 1088813). After the drying of several consecutive sections, non-sclerotic glomerular tufts were micro-dissected from PECs using a next generation micro-dissection system (MMI® CellCut Plus® Laser-Micro-dissection System, Molecular Machines and Industries (MMI), Eching, Germany and MMI® CellCamera DXA285cF) with the following settings: energy: 70%, focus: 50%, cutting speed: 40%.
Several non-deparaffinized and unstained consecutive transections were applied for the patient cohort II. Micro-dissection was performed by another micro-dissection system (Zeiss/Palm MicroBeam AxioVert 200M, Camera Hitachi HV-D30, Crystal Laser systems FTSS 355-50) with the same laser settings. For reference transcript analysis, 60 glomerular transections from all fifteen patient samples were used. In order to determine the effect of input RNA for valid mRNA expression analysis, 20, 60, 150 and 300 glomerular transections from eight patient samples out of cohort II were used.
RNA was quantified as previously described [17, 18]. Briefly, RNA was isolated by phenol–chloroform extraction and dissolved in 11 µl DEPC water. Subsequently, 1 µl was used for measuring the RNA amount and A260/A280 ratio by a spectrophotometer (Synergy, Biotek Instruments, Bad Friedrichshall, Germany). The remaining 10 µl RNA solution was reversely transcribed using the Multi-Scribe-based High Capacity Kit (ThermoFisher Scientific previously Applied Biosystems, Foster City, CA, USA). For quantitative mRNA analysis, we used reverse transcription single TaqMan real-time quantitative PCR (RT-qPCR) with TaqMan pre-designed assays after cDNA synthesis and preamplification (ThermoFisher Scientific). The TaqMan assay identity numbers and amplicon sizes of ADAMTS13 and selected candidate reference transcripts [PGK1 (Phosphoglyceratekinase 1), POLR2A, GAPDH, GUSB (Beta glucuronidase), PPIA (Peptidylprolyl isomerase A Alias: Cyclophilin A), ACTB (Beta-Actin), eukaryotic 18S rRNA, B2M (Beta-2-microglobulin), HPRT1 (Hypoxanthine phosphoribosyltransferase 1), RPLPO (ribosomal protein lateral stalk subunit P0), TBP (TATA-box binding protein) and HMBS (hydroxymethylbilane synthase)] as well as compartment markers WT1 (Wilms tumor 1) for podocytes/PECs, PAX2 (Paired box 2) for PECs/tubulointerstitium and GLEPP1 (Alias: PTPRO (Protein tyrosine phosphatase, receptor type O) for podocytes/PECs are detailed in Table 2. No-reverse transcription controls and no-template controls were included during each amplification step. Off-scale Cq values (Cq > 35) and missing Cq values were replaced by Cq = 36 .
Immunohistochemical staining was performed exemplarily in cases with ATN and mild interstitial inflammation (patient cohort I). All staining procedures were undertaken manually with the ABC detection system (peroxidase) from Vector Laboratories PK-6100 (Burlingame, CA, USA). Treating sections with 3% hydrogen peroxide for 10 min quenched endogenous peroxidase activity. Antigen retrieval for PAX2, WT1 and GLEPP1 was performed with Tris buffer pH 6.0 for 2.5 min at 120 °C in a pressure cooker. Secondary antibody was incubated for 30 min at room temperature at a dilution of 1:200. Diaminobenzidine (DAB) was used for a brown reaction product at the site of bound primary antibody.
In detail, for PAX2 immunostaining a rabbit polyclonal antibody (71-600, ThermoFisher, Foster City, CA, USA) was used at a dilution of 1:500. WT1 immunostaining was conducted with overnight incubation (4 °C) with a rabbit polyclonal antibody (RB-9267, ThermoFisher, Foster City, CA, USA) at a dilution of 1:10,000. GLEPP1 immunostaining was conducted with a rabbit polyclonal antibody (HPA034525, Sigma-Aldrich, Munich, Germany) at a concentration of 1:5000 with overnight incubation (4 °C).
Non-parametric Wilcoxon matched-pairs signed rank test was used for comparison of all data except for comparison of clinical data between cohort I and II (Mann–Whitney U test). Median values as well as 25th/75th percentiles are shown in the box plots or given in the text. Furthermore, absolute data or Spearman rank correlations are shown. Graphical and statistical analyses were performed using GraphPad Prism Version 5.0 (GraphPad Software Inc., La Jolla, CA, USA).
The NormFinder and geNormPlus approach was used as recently described . Briefly, NormFinder considers the intra- and inter-group variations independently. The stability value for each gene is a measurement of the estimated systematic error when using this gene for normalization. The commercially-available geNormPlus algorithm in qBaseplus, version 3.0 (http://www.biogazelle.com/qbaseplus Biogazelle, Zwijnaarde, Belgium) is based on calculating the average pairwise variations of each candidate reference gene compared to all candidate reference genes . For the identification of the optimal number of reference genes, geNormPlus provides the following algorithm: differences in the ratio of the average pairwise variation of normalization factors (geometric mean Cq value)  of consecutive candidate reference genes starting with the two most stable genes (n) and the addition of the next most stable gene (n + 1) until all genes have been added suggest the number of optimal reference genes. As a general guideline, the benefit of using an additional reference gene is limited as soon as the difference in the ratio (Vn/n + 1) drops below the 0.15 threshold .
According to our results of NormFinder and geNormPlus, we used the geometric mean of PGK1 and PPIA for normalization of glomerular ADAMTS13, GAPDH and POLR2A transcript expression. By contrast, for comparison of expression levels in PEC, tubulointerstitium and glomeruli we did not perform a NormFinder and geNormPlus analysis. Therefore, the geometric mean of GAPDH and PGK1 was used for normalization of cohort I. This decision was based on the fact, that GAPDH has been used for normalization in tubulointerstitial studies [23, 24].
Visualization of hemalaun-stained glomerular substructures and subsequent compartment-specific transcript analysis
One main goal was to push the visual resolution of glomerular structures to the limits, especially to prove the reliability of micro-dissection for the isolation of PECs. First, we compared non-deparaffinization (Fig. 1A, B) with deparaffinization by xylol followed by hemalaun staining (Fig. 1C, D). Non-deparaffinization is an excellent option for the micro-dissection of glomerular capillaries and to identify sclerosed glomeruli (Fig. 1E, F) or scarred tubulointerstitium. There was no difference in visualization comparing 2 and 3 µm sections, although non-deparaffinized 2 µm sections were much easier to dissect with the laser beam. Without deparaffinization (Fig. 1A, B), it was virtually impossible to differentiate PECs. However, following deparaffinization and a short hemalaun stain, it was feasible to detect nuclei of PECs and endothelial cells (Fig. 1C, D). For comparison, HE (Fig. 1G) or PAS (Fig. 1H, I) stains are given for the samples above (Fig. 1D–F).
Moreover, the high-resolution camera and touch screen from the MMI system enabled micro-dissecting PECs, sparing the deparaffinized and hemalaun-stained glomerular tuft, especially in cases with a wide Bowman’s space (Fig. 2). For further characterization of micro-dissected non-injured glomerular substructures, we used a set of selected mRNA transcripts: WT1 (podocytes), GLEPP1 (podocytes) and PAX2 (PECs). Significantly different expression levels of WT1, GLEPP1 and PAX2 indicated an appropriate separation of substructures (Fig. 3A–C). The geometric mean of PGK1 and GAPDH was used as normalization factor for relative expression analysis of WT1, GLEPP1 and PAX2. For clarification of mRNA expression levels, immunostaining of WT1, GLEPP1 and PAX2 is provided. Both, WT1 (Fig. 3A, D) and GLEPP1 (Fig. 3B, E) were strongly expressed in podocytes and to a lower extend also in PECs, but not in the tubulointerstitium. PAX2 was expressed in PECs and tubular epithelial cells (Fig. 3C, F).
RNA quantities from these hemalaun-stained glomerular tufts and micro-dissected PECs in association with the number of micro-dissected glomerular transections are shown in Fig. 4. Interestingly, there were no significant differences in RNA concentration between glomerular tufts [42.4 (18.8/97.8) ng/µl]; Fig. 4a) and Bowman’s capsules [43.8 (10.7/83.3) ng/µl]; Fig. 4b) in case of 335 (66/484) micro-dissected glomerular transections. Among patients with less micro-dissected glomerular transections (about 60 or less), RNA concentrations below 10 ng/µl in PECs and a broad range between 6 and 26 ng/µl in glomerular tufts were found. There was also one sample with no detectable amplification signal in reference and target transcripts in case of 65 micro-dissected Bowman’s capsule transections. However, high amounts of Bowman’s capsule transections (> 300) revealed sufficient RNA concentrations (> 30 ng/µl; Fig. 4b).
Reference transcripts for unstained glomerular transections
Two established algorithms (NormFinder and geNormPlus) were used to define stable glomerular reference transcripts in micro-dissected glomerular transections from fifteen unstained samples. We included twelve different mRNA candidate reference transcripts (Fig. 5 and Table 2). In one sample, HMBS was not detected and thus this sample was excluded from further stability analysis. Both, NormFinder (Fig. 5a) and geNormPlus (Fig. 5b) analysis recommended PGK1, PPIA, HPRT and GUSB as the top four glomerular reference transcripts (black squares). However, depending on the algorithm, the order of stability of these four reference transcripts varied. According to geNormPlus, the optimal number of reference transcripts was two (Fig. 5c), indicated by the ratio of the average pairwise variation (V) of the best two (V2) and best three (V3) reference transcripts, which was below 0.15.
For the comparison of different studies using human FFPE renal biopsies for glomerular micro-dissection, we performed a literature search in PubMed looking for articles published until March 2017. There are five studies using different reference transcripts or genes (GAPDH, ACTB and Eukaryotic 18S rRNA) [9,10,11,12,13]. Detailed differences between these studies and our study are highlighted in Table 3. GAPDH, ACTB and Eukaryotic 18S rRNA were not recommended as reference transcripts in our analysis of unstained glomerular transections.
mRNA expression in dependence of different amounts of unstained micro-dissected glomerular transections
In order to determine the mRNA minimum required for the detection of reliable mRNA expression in unstained glomerular transections, we used ADAMTS13, GAPDH and POLR2A as target transcripts and the geometric mean of PGK1 and PPIA as normalization factor. To this end, we compared different amounts of micro-dissected glomerular transections: 20, 60, 150 and 300. The total amount of isolated RNA was between 16.7 (13.1/24.1) ng/µl for 20 micro-dissected glomerular transections and 67.4 (53.6/86.3) ng/µl for 300 micro-dissected glomerular transections (Fig. 6a). As expected, there was a significant increase in RNA concentration from 60 glomerular transections upwards (Fig. 6a). However, the amount of RNA was not proportional to the number of glomerular transections. With our protocol, we obtained mRNA with a purity of A260/A280 ratio between 0.9 and 2.3 from all samples (Fig. 6b). There were no significant differences regarding mRNA purity between different amounts of micro-dissected glomerular transections (Fig. 6b). Comparing the effect of input mRNA, in groups with 20 glomerular transections five Cq values of ADAMTS13 had to be replaced by Cq = 36 resulting in a relative expression of < 0.01 in those five cases. Relative expression of ADAMTS13 displayed a significant difference between 150 and 300 glomerular transections (Fig. 6c). In groups with at least 60 glomerular transections the relative expression of ADAMTS13 displayed a lower inter-assay coefficient of variability compared to groups with 20 glomerular transections (Fig. 6d). For comparison expression levels of GAPDH and POLR2A are given (Fig. 6e, f). There were no significant differences for all three transcripts ADAMTS13, GAPDH and POLR2A comparing expression levels between 60 and 150 glomerular transections. However, for 300 glomerular transections there was a significant difference compared either to 60 glomerular transections (GAPDH) or to 150 glomerular transections (POLR2A). The lowest interquartile range for ADAMTS13, GAPDH and POLR2A were in samples with 150 glomerular transections. Expression differences between ADAMTS13, GAPDH and POLR2A revealed similar significance levels, when 60 and 150 or 60 and 300 glomerular transection counts were compared (data not shown).
Furthermore, we tested the reproducibility of mRNA transcript levels between different amounts of micro-dissected glomerular transections. A minimum of 60 micro-dissected glomerular transections revealed a significant correlation for POLR2A expression compared to 150 micro-dissected glomerular transections (Fig. 7a) as well as between 150 and 300 micro-dissected glomerular transections (Spearman r 0.833; p-value 0.010; data not shown). We could not find any further significant correlations comparing different amounts of glomerular transections (even between 150 and 300 transections) regarding the relative transcript expression of GAPDH and ADAMTS13 (Fig. 7b, c).
Both, RNA purity and quantity are decisive for successful mRNA expression studies. Our results confirm the previous finding, that the purity and quantity of RNA recovered from FFPE tissue are independent factors . The variability of RNA purity and quantity may be due to different pre-analytical handling of the biopsy such as time to formalin fixation or rapid versus ordinary embedding [26, 27]. Moreover, different mRNA specimens seem to be affected differently by degradation . Our goal was to document the RNA purity and quantity for quantitative mRNA analysis of micro-dissected glomerular compartments and to evaluate the approach according to the MIQE guidelines .
Despite the certain degradation of our RNA samples, we had a similar range for RNA purity in between different samples as well as between different amounts of RNA, which seems to be valuable for further analyses . As shown in Table 3, one recent study performed a comprehensive analysis for RNA purity, quality and quantity . Beside this study no other documents the A260/A280 ratio, which might be due to the fact that it is difficult or even impossible to obtain reliable absorption values at A260 in micro-dissected FFPE tissue . However, this does not necessarily mean that the RNA purity results are faulty. Landolt et al. have shown a mean A260/A280 ratio between 1.1 and 1.9 in 100 micro-dissected glomerular transections from FFPE biopsies , which was in the same range as in our samples. Most importantly, despite low A260/A280 ratios we could amplify mRNA transcripts, which demonstrates reasonable mRNA quality.
Another measure for RNA integrity is to calculate the ratio of 28S and 18S ribosomal RNA peak area . However, it has recently been postulated that this manner of RNA integrity assessment from degraded FFPE samples is not a sensitive measure of RNA quality and samples with low RNA integrity numbers are still usable for RT-qPCR, provided that the amplicon length is sufficiently short (< 90 nucleotides) . Therefore, we did not calculate RNA integrity.
Furthermore, studies have shown that the mean RNA fragment size—expressed as a DV200 value (percentage of RNA fragments > 200 nucleotides in size) and determined by a Fragment Analyzer (Advanced Analytical, Ankeny, USA) or Bioanalyzer (Agilent 2100, Agilent Technologies, Santa Clara, USA)—is a more reliable predictor for the successful RNA amplification of micro-dissected FFPE samples as compared to the conventional RNA integrity approach from regular tissue samples [9, 34,35,36]. In micro-dissected glomerular transections, there was a broad range for the DV200 values, from 26 to 89% in dependence of different RNA isolation methods . However, the RNA quantity from micro-dissected glomerular transections is even below the pre-determined standard range of the Fragment Analyzer . Therefore, due to the limited amount of RNA, we decided not to measure the DV200 value. Furthermore, in our opinion it is more important to insert all available micro-dissected RNA for reverse transcription and PCR itself. In case of 60 micro-dissected glomerular transections, we used a RNA concentration of 230 (182/277) ng/10 µl, which was a higher amount than requested for lowest RNA quality according to the stepwise approach of Landolt et al. . In general, they suggest that if the DV200 is of high quality (> 70% of fragments > 200 nucleotides), a smaller amount of RNA (about 20 ng) is required; for medium quality (50–70%), an intermediate amount of RNA (about 20–40 ng); and for low quality (30–50%), a high amount of RNA (40–100 ng) is required .
In order to document RNA yields and transection amounts, we used two different approaches. (I) For non-deparaffinized glomerular transections, we conducted a comprehensive comparative study of different transection amounts, which revealed increasing but no proportional RNA concentrations. The total RNA isolated from 60 micro-dissected glomerular transections is higher than the range described in the study of Landolt et al. , which obtained a range of 138 ng (± 43 ng) RNA per 100 micro-dissected glomerular transections using the ExpressArt kit (Amsbio, Abingdon, UK) for RNA extraction. Independent of the RNA concentration, we detected no significant differences for the relative expression of ADAMTS13 and POLR2A between 20, 60 and 150 micro-dissected glomerular transections. For comparison, another non-renal study demonstrated similar mRNA (Epidermal growth factor receptor, Human epidermal growth factor receptor 2, Mouse double minute 2 homolog) expression levels within the range of 106 to 100 micro-dissected cells in FFPE tissue . Especially, strongly-expressed transcripts like POLR2A provide reproducible relative expression levels between 60 and 150 glomerular transections. Moreover, there was a lower inter-assay coefficient of variability for ADAMTS13, when at least 60 unstained glomerular transections were used. Samples with 20 micro-dissected glomerular transections revealed a non-detectable amplification signal for ADAMTS13 in four and for GAPDH in three out of eight analyzed cases. In another sample with 20 glomerular transections there was an off-scale Cq value for ADAMTS13 expression. Therefore, we conclude that 60 glomerular transections with a median RNA concentration of [23 (18/28) ng/µl] are necessary and sufficient for mRNA expression analysis. This number of transections can be reasonably achieved in most routine biopsies.
(II) For the amount of deparaffinized and hemalaun-stained Bowman’s capsule transections, we can only provide a rough estimate, because we performed no comprehensive analysis for different amounts of Bowman’s capsule transections. In all samples the RNA concentration was higher than the detection limit of spectrophotometers (10 ng/µl) . About 300 hemalaun-stained micro-dissected Bowman’s capsule transections even revealed a higher amount of RNA compared to the mean concentration of 60 unstained glomerular transections. Therefore, we recommend micro-dissecting at least 300 hemalaun-stained Bowman’s capsule transections from FFPE renal biopsies. Interestingly, individual samples of our cohort with 180 Bowman’s capsule transections displayed higher RNA concentrations compared to 360 and 450 transections. Although samples were treated according to our diagnostic standard operation procedure, this unusual progression might be due to differences in pre-analytical steps, like the time period between harvesting and fixation of the tissue in formalin, which is crucial since any delay reduces the RNA quality and quantity .
In general, one has to keep in mind that concentration estimation can also be affected by phenol contamination, given that phenol absorbs strongly at 230 nm and ~ 270 nm. Indeed, in many samples we had high absorbance at 230 nm (data not shown). Therefore, the A260/A230 ratio was strongly reduced in our samples. However, there is no consensus on the acceptable lower limit of this ratio and the most important factor is the amount of phenol that is transferred to the downstream reaction, rather than the absorbance ratio. It is important to remember that absorbance ratios also depend on RNA concentration, especially in case of laser micro-dissection when RNA is at very low concentrations. Vice versa, absorbance around 260 nm can affect RNA concentration. Therefore, spectrophotometric analysis (e.g. Synergy, used in our study) of FFPE tissue yield different, consistently higher RNA concentrations compared to fluorometric analysis (e.g. Qubit) .
Furthermore, different RNA extraction methods yield different amounts of RNA . We used the phenol–chloroform extraction because in case of highly-degraded RNA like FFPE autopsy tissue the phenol–chloroform extraction is still a valid alternative . However, in whole non-renal tissue analysis, the High Pure kit (Roche, Mannheim, Germany) performed better than the conventional phenol–chloroform extraction method . In micro-dissected glomerular FFPE transections, the High Pure kit and ExpressArt kit (Amsbio, Abingdon, UK) recovered the largest amounts of RNA . Unfortunately, phenol–chloroform extraction was not used in their study . Until a study compares the High Pure kit or ExpressArt kit to the phenol–chloroform extraction method in micro-dissected glomerular FFPE transections, both methods can be considered adequate.
Operational issues for small amounts of starting material
To date, studies on PECs are limited because the small amounts of PECs could not be analyzed by quantitative real-time PCR for mRNAs. Our approach for laser micro-dissection can overcome this limitation. In contrast to all other studies shown in Table 3, we used the method of cDNA amplification. Especially for studying more than two targets and additional reference transcripts, we suggest increasing the amount of cDNA through this preliminary cDNA amplification step, which does not significantly distort relative mRNA levels [14, 17]. It has been shown that the combination of random primer-based cDNA synthesis and primer-specific preamplification provide best results . The proposed primer-specific cDNA synthesis is not an option in our samples . Furthermore, the variability of the preamplification—introduced into the experimental workflow of reverse transcription quantitative qPCR– is lower than the variability caused by the reverse transcription step . Moreover, the amount of preamplified transcripts correlates with the initial cDNA target copy numbers in both good quality  and poor quality samples . In order to avoid a bias, appropriate controls have to be used, like defined amounts of cell culture RNA, no-reverse-transcription controls and no-template controls.
Another aspect in case of micro-dissected glomerular transections with small amounts of starting material is how to manage reactions that do not give rise to any Cq value (missing Cq values). In order to avoid missing any low-level expressed transcripts, we suggest the pragmatic approach of replacing all off-scale Cq values (Cq > 35) and missing Cq values by Cq = 36 .
Expression analysis in a single glomerulus
For certain projects, it could be relevant to compare gene expression levels between a normal (or less injured) glomerulus and a sclerotic glomerulus in the same biopsy sample. The diameter of a human glomerulus is on average about 150 µm. Therefore, in theory about 75 transections (cut at a thickness of 2 µm) could be cut from one single glomerulus. Indeed, this number of transections could provide sufficient RNA. However, in practice it will be difficult to identify a single glomerulus over multiple section levels. Therefore, it would be an approach to pool glomeruli with normal (or less injured) morphology as well as sclerotic glomeruli. For example, we would suggest micro-dissecting about 20 transections and pooling three glomeruli per biopsy.
Twelve different candidate reference transcripts were tested for expression stability. All candidate reference transcripts are located on different chromosomes and are involved in different basic cellular processes, as shown in Table 2. This is a necessary requirement for reference transcripts . We suggest using the geometric mean of PGK1 and PPIA as normalization factor in unstained glomerular transections. PPIA was also suggested as a reference transcript for the comparison of malignant and non-malignant kidney specimens [45, 46]. In contrast to our study Schmid et al. have shown that PPIA was expressed at low levels in glomeruli and thus it was claimed not to be suitable as a reference transcript for glomerular analysis . One has to keep in mind that biopsies stored in RNA later were used in their study for micro-dissection. Furthermore, they demonstrated that gene expression studies of human kidney biopsies published from 1999 to 2002 predominantly used GAPDH as a reference transcript for RT-PCR studies . In our study GAPDH does not belong to the most stable glomerular transcripts. GAPDH, PPIA and PGK1 were also recommended as reference transcripts in studies about mice with cystic kidney disease . ACTB—used in some glomerular studies [12, 13, 48]—is no longer suggested as a reference transcript for RT-PCR in general . This is confirmed by our results of less stably ACTB. Eukaryotic 18S rRNA is used as reference transcript in one other renal study , although this transcript is less stably expressed in our cohort. HPRT1, RPLPO  and B2M  have been described as renal reference transcripts and were detected in renal tumors. Here, we describe for the first time that HPRT1 and GUSB belong to the top four glomerular reference transcripts in human non-neoplastic kidney biopsies. However, RPLPO and B2M are less suitable. GUSB has also been described as a suitable reference gene in metastatic renal cell carcinoma .
Limitations of the study
As discussed in Cohen et al. , a considerable variability of relative mRNA expression can be observed between consecutive sections. This variation could be a reason (I) for the decreasing inter-assay variability of relative ADAMTS13 expression with an increasing number of isolated glomerular transections and (II) for significant different expression levels of ADAMTS13, POLR2A and GAPDH in cases with 300 compared to 150 (ADAMTS13, POLR2A) or 60 (GAPDH) micro-dissected glomerular transections. This bias is negligible when relative expression values of different mRNA transcripts are compared in one defined amount of transections.
Furthermore, the number of biopsies for analyzing the effect of input RNA (patient cohort II) is difficult to increase due to the size limit of renal biopsies in clinical practice. However, even in our eight cases POLR2A revealed a reliable reproducibility for different amounts of RNA recovered from 60 to 150 micro-dissected glomerular transections when normalized to PGK1 and PPIA. There may be two reasons for the missing correlation of GAPDH and ADAMTS13 relative expression: (I) the small number of biopsies analyzed; or (II) the lower expression level of GAPDH and ADAMTS13 compared to POLR2A. Consequently, the inverse conclusion might be that fewer values of lower expressed transcripts have to be approached with caution.
Another aspect is the use of markers for discrimination between Bowman’s capsules and the tubulointerstitium. Due to the low amount of starting material, we focused on two positive markers being highly expressed in podocytes but being also expressed at lower level in Bowman’s capsules (WT1 and GLEPP1). However, these markers are not expressed in the tubulointerstitium. To the best of our knowledge, it is difficult to find negative markers for Bowman’s capsules, which could further exclude the contamination of tubular cells and interstitial cells in RNA from Bowman’s capsules.
Furthermore, due to the low amount of starting material we were not able to perform a comprehensive reference gene analysis in hemalaun-stained micro-dissected Bowman’s capsule transections. Therefore, we have chosen the geometric mean of PGK1 and GAPDH as normalization factor for PEC, tubulointerstitium and glomeruli on the basis of the current literature [23, 24].
In summary, we show the feasibility for micro-dissecting glomerular and Bowman’s capsule transections followed by evidence-based recommendations for quantitative mRNA analysis. We recommend micro-dissecting at least 60 unstained glomerular or 300 hemalaun-stained Bowman’s capsule transections from FFPE renal biopsies to obtain reliable results equivalent to larger amounts of tissue. We provide evidence that RNA quantity is within an acceptable range even for low RNA purity. Therefore, we suggest the following pragmatic approach: in case of 60 available unstained glomerular transections, quantity and purity control is optional for documentation but not gainful regarding further analytical steps, because (I) the range of RNA concentrations in equal amounts of micro-dissected glomerular transections is low and (II) appropriate normalization of Cq values using our suggested reference transcripts (PGK1 and PPIA) can compensate for different amounts of RNA purity and quantity . However, with fewer than 60 available unstained glomerular transections either quality and quantity control should be performed (e.g. measuring the DV200 value and fluorometric analysis) or the sample has to be excluded.
a disintegrin and metalloproteinase with thrombospondin type 1 motives 13
acute tubular injury
- Cq :
- DV200 :
percentage of RNA fragments > 200 nucleotides in size
hypoxanthine phosphoribosyltransferase 1
laser capture micro-dissection
Minimum Information for Publication of Quantitative Real-Time PCR Experiments
- MMI® :
Molecular Machines and Industries
paired box 2
polymerase chain reaction
parietal epithelial cells
phosphoglycerate kinase 1
RNA polymerase II subunit A
peptidylprolyl isomerase A
- PTPRO (alias GLEPP1):
protein tyrosine phosphatase, receptor type O, alias GLEPP1
ribosomal protein lateral stalk subunit P0
reverse-transcription quantitative PCR
TATA-box binding protein
average pairwise variation
wilms tumor 1
Espina V, Wulfkuhle JD, Calvert VS, VanMeter A, Zhou W, Coukos G, et al. Laser-capture microdissection. Nat Protoc. 2006;1(2):586–603.
Hussein K. Gene expression profiling in laser-microdissected bone marrow megakaryocytes. Methods Mol Biol. 2011;755:429–39.
Nagasawa Y, Takenaka M, Matsuoka Y, Imai E, Hori M. Quantitation of mRNA expression in glomeruli using laser-manipulated microdissection and laser pressure catapulting. Kidney Int. 2000;57(2):717–23.
Hewitson TD, Christie M, Smith ER. Laser capture microdissection of archival kidney tissue for qRT-PCR. Methods Mol Biol. 2016;1397:251–65.
Jonigk D, Modde F, Bockmeyer CL, Becker JU, Lehmann U. Optimized RNA extraction from non-deparaffinized, laser-microdissected material. Methods Mol Biol. 2011;755:67–75.
Specht K, Richter T, Muller U, Walch A, Hofler MW. Quantitative gene expression analysis in microdissected archival tissue by real-time RT-PCR. J Mol Med (Berl). 2000;78(7):B27.
Blakey GL, Laszik ZG. Laser-assisted microdissection of the kidney: fundamentals and applications. J Mol Histol. 2004;35(6):581–7.
De Spiegelaere W, Cornillie P, Van Poucke M, Peelman L, Burvenich C, Van den Broeck W. Quantitative mRNA expression analysis in kidney glomeruli using microdissection techniques. Histol Histopathol. 2011;26(2):267–75.
Landolt L, Marti HP, Beisland C, Flatberg A, Eikrem OS. RNA extraction for RNA sequencing of archival renal tissues. Scand J Clin Lab Invest. 2016;76(5):426–34.
Cohen CD, Grone HJ, Grone EF, Nelson PJ, Schlondorff D, Kretzler M. Laser microdissection and gene expression analysis on formaldehyde-fixed archival tissue. Kidney Int. 2002;61(1):125–32.
Schmid H, Henger A, Cohen CD, Frach K, Grone HJ, Schlondorff D, et al. Gene expression profiles of podocyte-associated molecules as diagnostic markers in acquired proteinuric diseases. J Am Soc Nephrol. 2003;14(11):2958–66.
Fries JW, Roth T, Dienes HP, Weber M, Odenthal M. A novel evaluation method for paraffinized human renal biopsies using quantitative analysis of microdissected glomeruli and VCAM-1 as marker of inflammatory mesangial cell activation. Nephrol Dial Transplant. 2003;18(4):710–6.
Bao H, Hu S, Zhang C, Shi S, Qin W, Zeng C, et al. Inhibition of miRNA-21 prevents fibrogenic activation in podocytes and tubular cells in IgA nephropathy. Biochem Biophys Res Commun. 2014;444(4):455–60.
Agustian PA, Bockmeyer CL, Modde F, Wittig J, Heinemann FM, Brundiers S, et al. Glomerular mRNA expression of prothrombotic and antithrombotic factors in renal transplants with thrombotic microangiopathy. Transplantation. 2013;95(10):1242–8.
Schmid H, Cohen CD, Henger A, Irrgang S, Schlondorff D, Kretzler M. Validation of endogenous controls for gene expression analysis in microdissected human renal biopsies. Kidney Int. 2003;64(1):356–60.
Ceol M, Del Prete D, Tosetto E, Graziotto R, Gambaro G, D’Angelo A, et al. GAPDH as housekeeping gene at renal level. Kidney Int. 2004;65(5):1972–3.
Theophile K, Jonigk D, Kreipe H, Bock O. Amplification of mRNA from laser-microdissected single or clustered cells in formalin-fixed and paraffin-embedded tissues for application in quantitative real-time PCR. Diagn Mol Pathol. 2008;17(2):101–6.
Bockmeyer CL, Jonigk D, Kreipe H, Lehmann U. MicroRNA profiling using RNA from microdissected immunostained tissue. Methods Mol Biol. 2011;755:85–94.
Stahlberg A, Rusnakova V, Forootan A, Anderova M, Kubista M. RT-qPCR work-flow for single-cell data analysis. Methods. 2013;59(1):80–8.
Bockmeyer CL, Sauberlich K, Wittig J, Esser M, Roeder SS, Vester U, et al. Comparison of different normalization strategies for the analysis of glomerular microRNAs in IgA nephropathy. Sci Rep. 2016;6:31992.
Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007;8(2):R19.
Ratert N, Meyer HA, Jung M, Lioudmer P, Mollenkopf HJ, Wagner I, et al. miRNA profiling identifies candidate mirnas for bladder cancer diagnosis and clinical outcome. J Mol Diagn. 2013;15(5):695–705.
Singh S, Manson SR, Lee H, Kim Y, Liu T, Guo Q, et al. Tubular overexpression of angiopoietin-1 attenuates renal fibrosis. PLoS ONE. 2016;11(7):e0158908.
Arvaniti E, Moulos P, Vakrakou A, Chatziantoniou C, Chadjichristos C, Kavvadas P, et al. Whole-transcriptome analysis of UUO mouse model of renal fibrosis reveals new molecular players in kidney diseases. Sci Rep. 2016;6:26235.
Chung JY, Braunschweig T, Hewitt SM. Optimization of recovery of RNA from formalin-fixed, paraffin-embedded tissue. Diagn Mol Pathol. 2006;15(4):229–36.
Srinivasan M, Sedmak D, Jewell S. Effect of fixatives and tissue processing on the content and integrity of nucleic acids. Am J Pathol. 2002;161(6):1961–71.
von Ahlfen S, Missel A, Bendrat K, Schlumpberger M. Determinants of RNA quality from FFPE samples. PLoS ONE. 2007;2(12):e1261.
Jung M, Schaefer A, Steiner I, Kempkensteffen C, Stephan C, Erbersdobler A, et al. Robust microRNA stability in degraded RNA preparations from human tissue and cell samples. Clin Chem. 2010;56(6):998–1006.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22.
Auer H, Lyianarachchi S, Newsom D, Klisovic MI, Marcucci G, Kornacker K. Chipping away at the chip bias: RNA degradation in microarray analysis. Nat Genet. 2003;35(4):292–3.
Hlubek F, Jung A. Fast protocol for DNA extraction from formalin-fixed paraffin-embedded tissues. In: Stanta G, editor. Guidelines for molecular analysis in archive tissues. Heidelberg: Springer; 2011. p. 41–4.
Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M, et al. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006;7:3.
Chung JY, Cho H, Hewitt SM. The paraffin-embedded RNA metric (PERM) for RNA isolated from formalin-fixed, paraffin-embedded tissue. Biotechniques. 2016;60(5):239–44.
Wang L, Nie J, Sicotte H, Li Y, Eckel-Passow JE, Dasari S, Vedell PT, Barman P, Wang L, Weinshiboum R, et al. Measure transcript integrity using RNA-seq data. BMC Bioinform. 2016;17:58.
Eikrem O, Beisland C, Hjelle K, Flatberg A, Scherer A, Landolt L, et al. Transcriptome sequencing (RNAseq) enables utilization of formalin-fixed, paraffin-embedded biopsies with clear cell renal cell carcinoma for exploration of disease biology and biomarker development. PLoS ONE. 2016;11(2):e0149743.
Endrullat C, Glokler J, Franke P, Frohme M. Standardization and quality management in next-generation sequencing. Appl Transl Genom. 2016;10:2–9.
Specht K, Richter T, Muller U, Walch A, Werner M, Hofler H. Quantitative gene expression analysis in microdissected archival formalin-fixed and paraffin-embedded tumor tissue. Am J Pathol. 2001;158(2):419–29.
Deben C, Zwaenepoel K, Boeckx C, Wouters A, Pauwels P, Peeters M, et al. Expression analysis on archival material revisited: isolation and quantification of RNA extracted from FFPE samples. Diagn Mol Pathol. 2013;22(1):59–64.
Okello JB, Zurek J, Devault AM, Kuch M, Okwi AL, Sewankambo NK, et al. Comparison of methods in the recovery of nucleic acids from archival formalin-fixed paraffin-embedded autopsy tissues. Anal Biochem. 2010;400(1):110–7.
Votavova H, Forsterova K, Stritesky J, Velenska Z, Trneny M. Optimized protocol for gene expression analysis in formalin-fixed, paraffin-embedded tissue using real-time quantitative polymerase chain reaction. Diagn Mol Pathol. 2009;18(3):176–82.
Glenn ST, Head KL, Teh BT, Gross KW, Kim HL. Maximizing RNA yield from archival renal tumors and optimizing gene expression analysis. J Biomol Screen. 2010;15(1):80–5.
Korenkova V, Scott J, Novosadova V, Jindrichova M, Langerova L, Svec D, et al. Pre-amplification in the context of high-throughput qPCR gene expression experiment. BMC Mol Biol. 2015;16:5.
Fox BC, Devonshire AS, Baradez MO, Marshall D, Foy CA. Comparison of reverse transcription-quantitative polymerase chain reaction methods and platforms for single cell gene expression analysis. Anal Biochem. 2012;427(2):178–86.
Li J, Smyth P, Cahill S, Denning K, Flavin R, Aherne S, et al. Improved RNA quality and TaqMan Pre-amplification method (PreAmp) to enhance expression analysis from formalin fixed paraffin embedded (FFPE) materials. BMC Biotechnol. 2008;8:10.
Bjerregaard H, Pedersen S, Kristensen SR, Marcussen N. Reference genes for gene expression analysis by real-time reverse transcription polymerase chain reaction of renal cell carcinoma. Diagn Mol Pathol. 2011;20(4):212–7.
Jung M, Ramankulov A, Roigas J, Johannsen M, Ringsdorf M, Kristiansen G, et al. In search of suitable reference genes for gene expression studies of human renal cell carcinoma by real-time PCR. BMC Mol Biol. 2007;8:47.
Cui X, Zhou J, Qiu J, Johnson MR, Mrug M. Validation of endogenous internal real-time PCR controls in renal tissues. Am J Nephrol. 2009;30(5):413–7.
Kohda Y, Murakami H, Moe OW, Star RA. Analysis of segmental renal gene expression by laser capture microdissection. Kidney Int. 2000;57(1):321–31.
Selvey S, Thompson EW, Matthaei K, Lea RA, Irving MG, Griffiths LR. Beta-actin-an unsuitable internal control for RT-PCR. Mol Cell Probes. 2001;15(5):307–11.
Gholami K, Loh SY, Salleh N, Lam SK, Hoe SZ. Selection of suitable endogenous reference genes for qPCR in kidney and hypothalamus of rats under testosterone influence. PLoS ONE. 2017;12(6):e0176368.
Wierzbicki PM, Klacz J, Rybarczyk A, Slebioda T, Stanislawowski M, Wronska A, et al. Identification of a suitable qPCR reference gene in metastatic clear cell renal cell carcinoma. Tumour Biol. 2014;35(12):12473–87.
Antonov J, Goldstein DR, Oberli A, Baltzer A, Pirotta M, Fleischmann A, et al. Reliable gene expression measurements from degraded RNA by quantitative real-time PCR depend on short amplicons and a proper normalization. Lab Invest. 2005;85(8):1040–50.
CLB: research design, laboratory work, data collection, data analysis, discussion of the results, writing of the manuscript, new reagents and analytic tools. JW: discussion of the results, laboratory work, data collection and data analysis. KS, PS, FM: laboratory work, data collection and data analysis. ME, PZ and KA: discussion of the results. CD: research design, laboratory work, discussion of the results, new reagents and analytic tools. All authors read and approved the final manuscript.
The excellent technical assistance of Andrea Kosel and Stefan Söllner (both Nephropathology Erlangen) is gratefully acknowledged.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated and analysed during this study are included in this published article and its additional information file (see Additional file 1).
Consent for publication
Ethics approval and consent to participate
All studies were carried out according to the latest revision of the Declaration of Helsinki (https://www.wma.net/policies-post/wma-declaration-of-helsinki-ethical-principles-for-medical-research-involving-human-subjects) and were approved by the Ethics Committee of Hannover Medical School (Reference Number 2241–2014) and Friedrich-Alexander-Universität Erlangen-Nuremberg (Reference Number 4415) waiving the need for retrospective consent for the use of archived residual FFPE tissue.
ELAN start-up funding (Fond for research and teaching at University Hospital Erlangen-Nuremberg, 15-06-14-1-Bockmeyer) for financing of a technical staff position is gratefully acknowledged. Furthermore, the authors would like to thank for funding for a laboratory rotation (Interdisciplinary Center for Clinical Research, University Hospital Erlangen-Nuremberg) and for consumables from the Dr. Werner Jackstädt-Stiftung for C.L.B. The authors acknowledge the support by Deutsche Forschungsgemeinschaft and Friedrich-Alexander-Universität Erlangen-Nuremberg within the funding programme Open Access Publishing. All funding bodies had no influence on the design of the study as well as collection, analysis, and interpretation of data and in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Bockmeyer, C.L., Wittig, J., Säuberlich, K. et al. Recommendations for mRNA analysis of micro-dissected glomerular tufts from paraffin-embedded human kidney biopsy samples. BMC Molecular Biol 19, 2 (2018). https://doi.org/10.1186/s12867-018-0103-x
- Laser capture micro-dissection
- Laser microdissection
- Kidney biopsies
- Endogenous controls
- Reference transcript
- Housekeeping gene