We describe a protocol for laser microdissection of sub-segments of the human kidney, including the glomerulus, proximal tubule, thick ascending limb, collecting duct and interstitium. The RNA is then isolated from the obtained compartments and RNA sequencing is carried out to determine changes in the transcriptomic signature within each sub-segment.
Gene expression analysis of human kidney tissue is an important tool to understand homeostasis and disease pathophysiology. Increasing the resolution and depth of this technology and extending it to the level of cells within the tissue is needed. Although the use of single nuclear and single cell RNA sequencing has become widespread, the expression signatures of cells obtained from tissue dissociation do not maintain spatial context. Laser microdissection (LMD) based on specific fluorescent markers would allow the isolation of specific structures and cell groups of interest with known localization, thereby enabling the acquisition of spatially-anchored transcriptomic signatures in kidney tissue. We have optimized an LMD methodology, guided by a rapid fluorescence-based stain, to isolate five distinct compartments within the human kidney and conduct subsequent RNA sequencing from valuable human kidney tissue specimens. We also present quality control parameters to enable the assessment of adequacy of the collected specimens. The workflow outlined in this manuscript shows the feasibility of this approach to isolate sub-segmental transcriptomic signatures with high confidence. The methodological approach presented here may also be applied to other tissue types with substitution of relevant antibody markers.
Technological advances in studying tissue specimens have improved understanding of the state of health and disease in various organs. Such advances have underscored that pathology can start in limited regions or in specific cell types, yet have important implications on the entire organ. Therefore, in the current era of personalized medicine, it is important to understand the biology at both the cell and regional level and not only globally1. This is particularly true in the kidney, which is composed of various specialized cells and structures that differentially initiate and/or respond to pathological stress. The pathogenesis of various types of human kidney disease is still not well understood. Generating a methodology to study changes in gene expression in specific tubular segments, structures or areas of the interstitium in the human kidney will enhance the ability to uncover region specific changes that could inform on the pathogenesis of disease.
Human kidney biopsy specimens are a limited and precious resource. Therefore, technologies interrogating transcriptomics in kidney tissue should be optimized to economize tissue. The available methods to study transcriptomics at the cell and regional level include single cell RNA sequencing (scRNaseq), single nuclear RNaseq (snRNaseq), in situ spatial hybridization, and laser microdissection (LMD). The latter is well suited for precise isolation of regions or structures of interest within tissue sections, for downstream RNA sequencing and analysis2,3,4,5. LMD can be adopted to rely on identification of specific cell types or structures based on validated markers using fluorescence-based imaging during dissection.
The unique features of laser microdissection assisted regional transcriptomics include: 1) the preservation of the spatial context of cells and structures, which complements single cell technologies where cells are identified by expression rather than histologically; 2) the technology informs and is informed by other imaging technologies because an antibody marker defines expression signatures; 3) the ability to identify structures even when markers change in disease; 4) detection of lowly expressed transcripts in approximately 20,000 genes; and 5) remarkable tissue economy. The technology is scalable to a kidney biopsy with less than 100 µm thickness of a core necessary for sufficient RNA acquisition and enables the use of archived frozen tissue, which are commonly available in large repositories or academic centers6.
In the ensuing work, we describe the regional and bulk transcriptomics technology in detail, optimized with a novel rapid fluorescence staining protocol for use with human kidney tissue. This approach improves upon classic LMD explorations because it provides separate expression data for the interstitium and nephron sub-segments as opposed to aggregate tubulointerstitial expression. Included are the quality assurance and control measures implemented to ensure rigor and reproducibility. The protocol enables visualization of cells and regions of interest, resulting in satisfactory acquisition of RNA from these isolated areas to allow downstream RNA sequencing.
The study was approved for use by the Institutional Review Board (IRB) at Indiana University.
NOTE: Use this protocol with kidney nephrectomy tissue (up to 2 cm in both the X and Y dimensions) preserved in the Optimal Cutting Temperature (OCT) compound and stored at -80 °C. Perform all work in a manner that limits RNA contamination, use clean disposable gloves and a face mask. Ensure the cleanliness of all surfaces. The equipment for which this protocol was optimized is a laser microdissection system featuring pulsed UV laser.
1. Cryosectioning
2. Laser microdissection
3. RNA isolation
NOTE: For this RNA isolation protocol, we have adapted a protocol from a commercial RNA isolation kit. The manufacturer’s protocol has been modified to meet the quality control requirements set for the project.
4. RNA sequencing
Samples
We present data from nine reference nephrectomies (3 specimens obtained at Indiana University and 6 specimens obtained through Kidney Precision Medicine Project), utilizing the rapid fluorescence staining protocol to isolate kidney nephron segments and interstitial areas. The sections utilized in this process were obtained from deceased kidney donors or unaffected tumor nephrectomies. These samples did not have pathologic evidence of disease as visualized in the H&E stain taken from contiguous sections (Figure 1A).
Laser Microdissection Quality Control
Identification of tubular sub-segments in the kidney is accomplished through antibody staining of unique tubular markers, as well as morphology and structural landmarks. Figure 1B-C illustrates the staining and microdissection of tubular sub-segments from a representative nephrectomy. This is accomplished by staining with DAPI (nuclei), FITC-Phalloidin (F-actin), and an additional antibody as necessary. The fluorescence staining used made it possible to visualize renal sub-segments with a high degree of confidence, further guided by morphological features as well as spatial positioning of the imaged structures (Figure 2). For example, glomeruli are revealed by phalloidin and DAPI staining with very distinctive morphology. Similarly, collecting ducts are visualized using nuclear morphology and the absence of other stains. The megalin/LRP2 antibody (directly conjugated to Alexa Fluor 568) is used here to identify proximal tubules. The thick ascending limb is visualized using a Tamm-Horsfall protein (THP) antibody (directly conjugated to Alexa Fluor 546). In contrast, the interstitium is excised along the outer membrane of the tubules and includes stromal and immune cells as well as small capillaries.
To obtain sufficient RNA of 0.5 – 1 ng per sub-segment, we seek to dissect a minimum area of 500,000 µm2. This generally requires five to eight 12 µm thick sections (1 section per slide) to obtain sufficient material for all sub-segments. Dissection of any slide is completed in less than 2 hours to preserve RNA quality which allows one to cut ~4 segments per slide. Acquired tissue from the same sub-segment is pooled from multiple slides for downstream sequencing. A pre- and post-LMD immunofluorescence image is obtained to validate the collection of sub-segments of interest using an attached camera. Figure 3 delineates the success rates of dissection area and minimum RNA input. The success rate of meeting the minimum area input was >90% in this representative dataset. If the source tissue has a small amount of the CD or interstitium, there is a possibility the dissection area will be below 500,000 µm2 for these segments after utilizing 8 slides. Additional slides could be cut depending on the needs of the user; however, as seen in Figure 3, if the area is close to 500,000 µm2, the desired genes detected count has a high success rate (100%) and the total read count success rate is 98.6% (1 specimen fell below the QC metric). Thus, there was enough tissue to achieve sufficient gene and read counts without utilizing additional tissue.
Sequencing Input Quality Control
The selected commercial library preparation and sequencing platform allows reproducible measurements of transcript expression, even with low quantities of highly degraded RNA. The minimum RNA input is 500 pg and the minimum DV200 (proportion of reads longer than 200 nucleotides in length) is above 25%. There is no minimum RIN. Sequencing with 100 million reads per sample, we seek to saturate the available reads in order to detect lowly expressed transcripts. The total reads can be adjusted according to the needs of the user.
It is straightforward to obtain sufficient RNA from two 12 µm bulk cross-sectional sections, so we seek to obtain a higher minimum total mRNA quantity of 4 ng for bulk sequencing. As delineated in the protocol, each sub-segment requires 0.5-1 ng in total RNA. The optimum RNA concentration is above 50 pg/µL after isolation. RNA amounts lower than this may lead to reduced gene detection and read counts observed. The majority of the samples yield >20,000 genes detected and >1 million reads. A DV200 greater than 25% for all specimens is required by the manufacturer (>75% is considered optimal). Although RIN is measured, the process of laser microdissection reduces RIN and the RNA sequencing platforms have been optimized for fragmented RNA. RNA quantity and quality are assessed by a bioanalyzer prior to sequencing. The rapid stain decreases the amount of time the tissue is exposed to room temperature and aqueous conditions, thereby minimizing to the extent possible, the degradation of RNA. The DV200 quality control metric for sequencing input is also found in Figure 3.
Sequencing Output Quality Control
Downstream data processing uses quantile normalization to allow comparison between samples in different batches. A minimum gene detection count of 10,000 and a minimum read count of 1 million are employed as thresholds to include samples in quantile normalization. Samples with lower gene detection or read counts are excluded from quantile normalization and subsequent comparative analyses. Only 1 out of 98 dissected sub-segmental samples failed to meet these thresholds (Figure 3). The Q30 (proportion of reads mapped at a 99.9% confidence) has been greater than 93% for each sequencing experiment (Figure 4).
We sequence a human reference RNA (25 ng) with each sequencing run to allow correction for batch effect if such correction is desired or required. The measured batch effect should be within 1 standard deviation of mean expression with an R value > 0.9. If batch effect is higher, additional normalization is required for that sequencing run. Here, the batch effect has been below 3% in the four sequencing experiments depicted in Figure 4. Thus, batch effect is typically addressed with quantile normalization for all sequencing runs, agnostic to the reference RNA. No correction based on the reference RNA has been implemented yet, but this information is available and could be used depending on the goals of a particular analysis.
Rigor and Reproducibility
It is important to demonstrate rigor in identifying cell types and structures. After laser microdissection, we test for the enrichment of known markers in the glomeruli, PT, TAL, CD, and interstitium as compared to the other compartments (Table 1). An enrichment analysis is used to compare gene expression for pre-determined expected markers. Gene expression is conveyed as log2 ratios of the expression of the sub-segment of interest compared to the mean of the other sub-segments. This expression metric was chosen for human samples as it reduces inter-individual variability, facilitating comparisons between cell and compartment types across specimens. However, the raw reads and quantile normalized reads may also be compared in alternative analyses. Figure 5 illustrates examples of immunohistochemical staining of these 5 selected known markers, as presented in the Human Protein Atlas.
Thus, we present orthogonal data sources which lend confidence to the correct sub-segment being collected: 1) the imaging which include the antibody stain and morphology of the segment/compartment, and 2) the expression output showing known markers are expressed in the corresponding sub-segment.
The regional transcriptomic analysis of genes expressed in each sub-segment allows for novel and underappreciated marker identification. Figure 6 illustrates an example of one marker for each sub-segment identified through LMD regional transcriptomics that also yielded a specific immunohistochemical staining of the corresponding nephron sub-segment.
To demonstrate reproducibility and understand the effect of operator variability, we conducted an experiment on a tumor nephrectomy sample. This specimen had the glomeruli, PT, and TAL re-dissected to increase the confidence of the RNA sequencing results and now serves as a useful technical replicate. Months apart, this sample underwent a second instance of cryosectioning, antibody staining, LMD dissection, RNA extraction, library prep, and RNA sequencing (Figure 7) of the glomeruli, PT and TAL compartments. The two versions (v1 and v2) were compared. The glomerular, PT, and TAL compartments were highly correlated with r2 values >0.95, similar to the minimal batch effect of our reference RNA.
Figure 1: Representative images of the renal tissue.
(A) H&E stain of a human reference kidney tissue. (B) Immunofluorescent image of human reference kidney tissue with stained thick ascending limb segments prior to LMD and (C) immediately following dissection of single thick ascending limb structure. Scale bar = 100 µm. Please click here to view a larger version of this figure.
Figure 2: Representative images of reference renal sub-segments, visualized at 20x using a specific immunofluorescent staining approach.
(A)a Glomerulus, (B) Proximal Tubule, (C) Thick ascending limb, (D) Collecting duct, (E) Interstitium. (*= p < 0.05, **= p < 0.01, ***= p < 0.001 by ANOVA). Scale bar = 100 µm. Please click here to view a larger version of this figure.
Figure 3: Quality control metrics in segmental transcriptomics.
(A) Greater than 90% of pilot samples met the dissection area input threshold of 500,000 µm2. (B) All samples except one met the desired RNA input of at least 500 pg. (C) 100% of the pilot samples met the manufacturer’s minimum DV200 threshold of 25% and (D) 100% of samples met our minimum threshold of 10,000 genes detected. (E) All samples except one reached a minimum total read count of 1 million. As expected, lower dissection areas correlate with reduced RNA concentrations, lower gene counts, and lower read counts. Please click here to view a larger version of this figure.
Figure 4: Sequencing Quality Control.
(A) Sequencing runs employing a reference RNA (not quantile normalized) are compared to mean expression of all runs. There is strong correlation with minimum batch effect (all R >0.969). (B) Greater than 93% of our reads in all runs were mapped with high confidence (Q30 or 99.9%). Note that Q30 values are provided on a per lane basis with multiple lanes per run. Please click here to view a larger version of this figure.
Figure 5: Examples of immunohistochemical staining to identify the selected known markers.
Images are depicted for (A) NPHS1 (B) LRP2 (C) UMOD (D) SLC4A9 (E) COL6A2. Immunohistochemistry images are obtained from the Human Protein Atlas. Please click here to view a larger version of this figure.
Figure 6: Examples of immunohistochemical staining to illustrate expression of nontraditional marker genes in relevant sub-segments.
Depicted transcripts with corresponding immunohistochemistry include (A) SHANK3 in glomerulus, (B) ACSM2B in proximal tubule, (C) RAP1GAP in the thick ascending limb, and (D) L1CAM in the collecting duct. Immunohistochemistry images are obtained from the Human Protein Atlas. (*= p < 0.0001 by ANOVA) Please click here to view a larger version of this figure.
Figure 7: Correlation between two distinct dissections with separate sequencing of the same sample.
A high degree of correlation was observed between different dissections in the (A) glomerulus, (B) proximal tubule, and (C) thick ascending loop of Henle. Please click here to view a larger version of this figure.
Marker | Segment | Fold Change | p-Value |
NPHS1 | Glomerulus | 32.64 | 1.97E-08 |
LRP2 | Proximal Tubule | 4.69 | 1.21E-05 |
UMOD | Thick Ascending Limb | 24.92 | 2.00E-07 |
SLC4A9 | Collecting Duct | 4.09 | 0.000133 |
COL6A2 | Interstitium | 7.19 | 1.47E-06 |
Table 1: Representative average values for each sub-segment of interest among three nephrectomy samples, compared to the remaining sub-segments.
LMD based transcriptomics is a useful technology that anchors gene expression to specific areas within the tissue. The basis of this technology and its potential application in the kidney has been described previously8. However, optimization, modernization and streamlining of fluorescence-based dissection specifically aimed at high accuracy dissection for downstream RNA sequencing is less ubiquitous. Because this methodology is spatially grounded within the tissue, it has the potential to reveal novel or underappreciated markers in tissue structures, especially in disease states. In fact, many common markers of cells may change with disease, hence relying only on the cell RNA expression markers for identity adjudication without the spatial context may be challenging in disease states. Therefore, the spatially anchored LMD approach is likely to complement and provide a ground-truth platform for many other transcriptomic technologies like single cell and single nuclear RNA sequencing, which define cells by predominantly by their gene expression profile9. Furthermore, the depth of gene expression provided by LMD (up to 20,000 genes per sample) provides another advantage and important venue to cross-link data from multiple sources and account for changes in gene expression that may not be apparent without a certain depth. The consideration of tissue economy is another positive feature of this technology, whereby a thickness of less than 100 µm total from a frozen block could be sufficient for an entire LMD dataset. This allows leftover tissue to be used for other analytical purposes.
The LMD has limitations. One anticipated limitation of laser microdissection transcriptomic data acquisition is its lower throughput nature. Future automation may prove important in improving scalability10,11. Two additional limitations of LMD transcriptomics include the dependency on expertise and the impact of tissue quality on data results. Inadvertent collection of cells and material outside the segment of interest is a known limitation because enriched compartments of cells are collected, not a single cell. LMD relies on a user’s proficiency in understanding the tissue structure and therefore requires a certain domain expertise. Thus, individuals must be trained to identify relevant regions in the kidney with and without antibody staining. Finally, the effect of tissue quality may impact downstream data quality. The LMD protocol leads to tissue degradation, so the quality of starting material is important. However, this protocol was optimized on archived biopsies with several samples of only moderate quality. The data included shows that the selected transcriptomic platform is robust.
Of note, the isolation kit used and the downstream methodology of RNA sequencing must be tailored specifically for the expected degree of RNA degradation. The staining protocol described here was adopted because it had the most favorable effect on minimizing such effect. In this protocol, tissue was embedded in OCT in order to facilitate future comparisons across other transcriptomic technologies. However, formalin fixed tissue may be an alternative.
The data presented in this manuscript reveal the benefits of regional transcriptomics, including the optimization, quality control metrics and technology outcomes. The establishment of rigorous pre-analytical and analytical quality control criteria, and also the establishment of solid evidence of rigor and reproducibility is essential for any methodology. Future applications of regional transcriptomics can be broadly grouped into biologic applications, technical advancements, and analytic progress. Future biologic applications may be aimed at interrogation of tissue from uninjured, injured, and regenerating tubules, as well as tubules in close proximity to inflammation. These compartments can be identified in the same biopsy both by traditional morphological markers as well as antibody stains for “pathway” specific markers. The two antibody markers validated for use in this technology, megalin and uromodulin, are highly expressed in the proximal tubule and thick ascending limb respectively. However, it is possible that these markers may be reduced in severely diseased tissue. In these situations, additional antibody validation of pathway specific markers or novel markers that do not change their expression level may overcome this issue.
LMD is likely to be an appropriate method to use in transcriptomic interrogation of other organs. The selection and optimization of antibodies that work for rapid staining will be imperative. An antibody that works in a regular staining protocol may not necessarily work for the rapid staining protocol, which likely requires high affinity antibodies. A primary conjugated antibody will minimize the steps needed and may offer advantages. However, conjugation of an antibody to fluorophores may alter the binding, and pilot experiments need to be performed before incorporation of these reagents in the protocol.
In conclusion, we describe a pipeline for fluorescence-based laser microdissection to isolate specific areas and structures in the human kidney to enable a transcriptomic analysis. This approach offers important advantages and can be extended to other tissues to provide a tissue based molecular interrogation. Regional transcriptomics complements other transcriptomic technologies by providing a histopathologic anchor to understand expression, assisting in the interpretation of the biology and the signature of health and disease. This technology fits naturally within the vision for building a transcriptomic atlas of the kidney.
The authors have nothing to disclose.
General: The authors would like to thank the investigators of the Kidney Precision Medicine Project (www.kpmp.org) for their gracious support and advice.
Funding: Support for this work was provided by the NIH/NIDDK K08DK107864 (M.T.E.); NIH/NIDDK UG3DK114923 (T.M.E., P.C.D.); R01DK099345 (T.A.S.). Research reported in this manuscript was supported by the National Institute of Diabetes and Digestive and the Kidney Diseases (NIDDK) Kidney Precision Medicine Project (KPMP), (www.kpmp.org), under award number U2CDK114886.
Data and materials availability: Data is archived in the Gene Expression Omnibus (GEO # pending)
Acetone | Sigma-Aldrich | 270725-1L | |
AMPure Beads | Beckman Coulter | A63880 | |
Bioanalyzer | Agilent | 2100 | |
BSA | VWR | 0332-100G | |
DAPI | ThermoFisher | 62248 | |
Desiccant Cartridge | Bel-Art | F42046-0000 | |
DNAse | Qiagen | 79254 | RDD buffer is included in the pakage |
Laser Microdissection Microscope | Leica | LMD6500 | |
Megalin/LRP2 Antibody | Abcam | ab76969 | Directly conjugated to Alexa Fluor 568 |
Microcentrifuge tubes | ThermoFisher | AB-0350 | |
Microscope camera | Leica | DFC700T | |
PBS (RNAse Free) | VWR | K812-500ML | |
Phalloidin (Oregon Green 488) | ThermoFisher | O7466 | |
PicoPure RNA Isolation Kit | Applied Biosystems | KIT0204 | |
PPS-membrane slides | Leica | 11505268 | |
qPCR Human Reference Total RNA 25 µg | Takara Clontech | 636690 | |
RNA 6000 Eukaryote Total RNA Pico Chip | Agilent | 5067-1513 | |
RNAse Away | ThermoFisher | 7000 | |
RNAse Inhibitor | ThermoFisher | AM2696 | |
Sequencer (HiSeq or NovaSeq) | Illumina | NA | |
SMARTer Stranded Total RNAseq Pico Input v2 | Takara Clontech | 634411 | |
Tamm-Horsfall Protein Antibody | R&D Systems | AF5144 | Directly conjugated to Alexa Fluor 546 |
Tissue-Tek® O.C.T. Compound | Sakura | 4583 |