A subscription to JoVE is required to view this content. Sign in or start your free trial.

In This Article

  • Summary
  • Abstract
  • Introduction
  • Protocol
  • Results
  • Discussion
  • Disclosures
  • Acknowledgements
  • Materials
  • References
  • Reprints and Permissions

Summary

Here we describe a method to assess lung expression of miRNAs that are predicted to regulate inflammatory genes using mice exposed to ozone or filtered air at different stages of the estrous cycle.

Abstract

MicroRNA (miRNA) profiling has become of interest to researchers working in various research areas of biology and medicine. Current studies show a promising future of using miRNAs in the diagnosis and care of lung diseases. Here, we define a protocol for miRNA profiling to measure the relative abundance of a group of miRNAs predicted to regulate inflammatory genes in the lung tissue from of an ozone-induced airway inflammation mouse model. Because it has been shown that circulating sex hormone levels can affect the regulation of lung innate immunity in females, the purpose of this method is to describe an inflammatory miRNA profiling protocol in female mice, taking into consideration the estrous cycle stage of each animal at the time of ozone exposure. We also address applicable bioinformatics approaches to miRNA discovery and target identification methods using limma, an R/Bioconductor software, and functional analysis software to understand the biological context and pathways associated with differential miRNA expression.

Introduction

microRNAs (miRNAs) are short (19 to 25 nucleotides), naturally occurring, non-coding RNA molecules.Sequences of miRNAs are evolutionary conserved across species, suggesting the importance of miRNAs in regulating physiological functions1. microRNA expression profiling has been proven to be helpful for identifying miRNAs that are important in the regulation of a variety of processes, including the immune response, cell differentiation, developmental processes, and apoptosis2. More recently, miRNAs have been recognized for their potential use in disease diagnostics and therapeutics. For researchers studying mechanisms of gene regulation, measuring miRNA expression can enlighten systems-level models of regulatory processes, especially when miRNA information is merged with mRNA profiling and other genome-scale data3. On the other hand, miRNAs have also been shown to be more stable than mRNAs in a range of specimen types and are also measurable with greater sensitivity than proteins4. This has led to considerable interest in the development of miRNAs as biomarkers for diverse molecular diagnostic applications, including lung diseases.

In the lung, miRNAs play important roles in developmental processes and the maintenance of homeostasis. Moreover, their abnormal expression has been associated with the development and progression of various pulmonary diseases5. Inflammatory lung disease induced by air pollution has demonstrated greater severity and poorer prognosis in females, indicating that hormones and the estrous cycle can regulate lung innate immunity and miRNA expression in response to environmental challenges6. In this protocol, we use ozone exposure, which is a major component of air pollution, to induce a form of lung inflammation in female mice that occurs in the absence of adaptive immunity. By using ozone, we are inducing the development of airway hyperresponsiveness that is associated with airway epithelial cell damage and an increase in neutrophils and inflammatory mediators in proximal airways7. Currently, there are not well-described protocols to characterize and analyze miRNAs across the estrous cycle in ozone-exposed mice.

Below, we describe a simple method to identify estrous cycle stages and miRNA expression in lung tissue of female mice exposed to ozone. We also address effective bioinformatics approaches to miRNA discovery and target identification, with an emphasis on computational biology. We analyze the microarray data using limma, an R/Bioconductor software that provides an integrated solution for analyzing data from gene expression experiments8. Analysis of PCR array data from limma has an advantage in terms of power over t-test based procedures when using small number of arrays/samples to compare expression. To comprehend the biological context of miRNA expression results, we then used the functional analysis software. In order to understand the mechanisms regulating transcriptional changes and to predict likely outcomes, the software combines miRNA-expression datasets and knowledge from the literature9. This is an advantage when compared with software that just look for statistical enrichment in overlapping to sets of miRNAs.

Protocol

All methods described here have been approved by the Institutional Animal Care and Use Committee (IACUC) of Penn State University.

1. Assessment of the Estrous Cycle Stage

  1. Properly restrain a female C57BL/6 mouse (8–9 weeks old) using the one-handed mouse restraint technique described in Machholz et al.10.
  2. Fill the sterile plastic pipette with 10 μL of ultra-pure water.
  3. Introduce the tip of plastic pipette into the vagina.
  4. Gently flush the liquid 4–5 times to collect the sample.
  5. Place the final flush containing vaginal fluid on a glass slide.
  6. Observe the unstained vaginal flush under a light microscope with a 20x objective.
    NOTE: Animals that do not show regular cycles due to pseudopregnancy or other causes need to be excluded from the experiment. It is recommended to perform daily vaginal secretions for at least three consecutive cycles to confirm cyclicity.

2. Exposure to Ozone

  1. Place a maximum of 4 mice in two 1.2 L glass containers with wire mesh lids, and water ad libitum.
  2. Put one glass container in the ozone chamber and the other in the filtered air exposure chamber.
  3. Adjust the ozone concentration to 2 ppm and monitor ozone levels regularly.
    NOTE: The ozone apparatus delivers a regulated airflow (>30 air changes/h) with controlled temperature (25 °C) and relative humidity (50%). The system generates ozone by an electrical discharge ozonizer, which is monitored and controlled by an ultraviolet ozone analyzer and mass flow controllers, as described previously11.
  4. Remove glass containers after 3 h of ozone/filtered air exposure. Return animals to cage with bedding, food, and water ad libitum.

3. Lung Collection

  1. 4 h after exposure, anesthetize animals with an intraperitoneal injection of a ketamine/xylazine cocktail (90 mg/kg ketamine, 10 mg/kg xylazine).
    NOTE: To confirm a proper level of anesthesia, check the mouse for a pedal reflex (firm toe pinch) and adjust the anesthetic as needed.
  2. Dampen the mouse skin with 70% ethanol.
  3. Make a 2 cm midline incision using an operating scissor and surgical tweezers to expose the vena cava.
  4. Sacrifice mice by transection of the vena cava and aorta. If needed, insert a 21 G gauge needle into the vena cava above the renal veins to collect blood prior to exsanguination. Alternatively, collect blood via heart puncture following standard protocols.
  5. Use a surgical scissor to cut open the abdominal cavity and remove skin/upper muscle, moving upwards toward the ribs.
  6. Use a surgical scissor to puncture the diaphragm.
    NOTE: Lungs will collapse away from the diaphragm.
  7. Cut away the ribcage using a surgical scissor to expose the heart and lungs.
  8. Using forceps, take a 1.5 mL RNase-free microcentrifuge tube and submerge it in liquid nitrogen to fill the tube.
    NOTE: Use goggles and protective gloves to handle liquid nitrogen.
  9. Remove the lungs, place them into the RNase-free microcentrifuge 1.5 mL tubes filled with liquid nitrogen to snap-freeze the tissue, and wait a few seconds until the liquid evaporates.
  10. Close the tube lid and store the tissue at -80 °C until use.

4. RNA Preparation

  1. Pulverize whole lungs using a stainless-steel tissue pulverizer.
    NOTE: The tissue pulverizer needs to be placed in liquid nitrogen prior to use. Clean the pulverizer after each use with RNAse solution.
  2. Split the pulverized lungs and place into two 1.5 mL tubes (half lung each).
  3. Add 500 µL of guanidinium thiocyanate per sample tube and mix.  Homogenize each sample using 18 G, 21 G, and 23 G needles, respectively.
    NOTE: Samples can be spiked with 5.6  x  108 copies of a small RNA spike-in control (from a different species) before proceeding with the extraction.
  4. Add 500 µL of ethanol to each sample and vortex for 15 s.
  5. Load the mixture in a spin column in a collection tube and centrifuge at 12,000 x g for 1 min. Discard the flow-through.
  6. For DNase I treatment (in-column);
    1. Add 400 µL of RNA Wash Buffer and centrifuge at 12,000 x g for 1 min.
    2. In an RNase-free tube, add 5 µL of DNase I and 75 µL of 1x DNA digestion buffer and mix. Add the mix directly to the column matrix.
    3. Incubate at room temperature for 15 min.
  7. Add 400 µL of RNA prewash solution to the column and centrifuge at 12,000 x g for 1 min. Discard the flow-through and repeat this step.
  8. Add 700 µL of RNA wash buffer to the column and centrifuge at 12,000 x g for 1 min. Discard the flow-through.
  9. Centrifuge at 12,000 x g for 2 min to remove remaining buffer. Transfer the column into an RNase-free tube.
  10. To elute RNA, add 35 µL of DNase/RNase-free water directly to the column matrix and centrifuge at 12,000 x g for 1.5 min.
  11. Measure total RNA concentration (260 nm) and purity using a spectrophotometer. Follow instructions to perform RNA quantification in a 1.5 µL sample aliquot. Blank the instrument with DNase/RNase-free water used for elution.
    NOTE: A 260/280 ratio of ~2.0 is generally accepted as “pure” for RNA. Typical RNA concentration usually ranges between 750 and 2,500 ng/µL.
  12. Store at -80 °C.

5. miRNA Profiling

  1. To retro-transcribe small RNAs, use 200 ng of total RNA.
    1. Prepare the reverse-transcription reaction mix on ice (total volume per reaction is 20 µL). For each reaction, add 4 µL of 5x buffer, 2 µL of 10x nucleotides mix, 2 µL of reverse transcriptase, and 2 µL of RNase-free water. Mix all the components and aliquot in 600 µL RNAse-free plastic tubes (10 µL of mix per reaction).
      NOTE: The reverse-transcription master mix contains all components required for first-strand cDNA synthesis except template RNA.
      NOTE: Calculate excess volume (10%) when preparing the master mix.
    2. Add the template RNA (200 ng in 10 µL) to each tube containing reverse-transcription master mix. Mix, centrifuge for 15 s at 1,000 x g, and store them on ice until placing in thermocycler or dry block.
    3. Incubate for 60 min at 37 °C.
    4. Incubate for 5 min at 95 °C and place the tubes on ice.
    5. Dilute the cDNA by adding 200 µL of RNase-free water to each 20 µL reverse-transcription reaction.
  2. Perform real-time PCR using the Mouse Inflammatory Response and Autoimmunity miRNA PCR Array.
    1. Prepare a reaction mix (total volume 1100 µL): For each reaction, add 550 µL of 2x PCR master mix, 110 µL of 10x universal primer mix, 340 µL of RNase-free water, and 340 µL of template cDNA (diluted reaction from step 5.1.5).
    2. Add 10 µL of reaction mix to each well of the pre-loaded miRNA PCR Array using a multichannel pipettor.
    3. Seal the miRNA PCR Array plate with optical adhesive film.
    4. Centrifuge the plate for 1 min at 1,000 x g at room temperature to remove bubbles.
    5. Program the real-time cycler, PCR initial activation step for 15 min at 95 °C, 3-step cycling containing denaturation for 15 s at 94 °C, annealing for 30 s at 55 °C, and extension for 30 s at 70 °C for 40 cycles number.
      NOTE: Follow manufacturer’s cycling conditions instructions to set up the real-time cycler. Perform the dissociation curve step built into the real-time cycler software.
    6. Perform data analysis.

6. Data Analysis

  1. Extract Ct values from the real time PCR software for each sample into an analysis software.
    NOTE: A Ct value of 34 is considered as cutoff. If samples contain spike-in control (e.g., cel-mir-39), normalize Ct values to the spike in control for each sample. Threshold values may need to be manually set. Baseline values are automatically set.
  2. Normalize Ct values to the average Ct of six miRNA housekeeping controls: SNORD61, SNORD68, SNORD72, SNORD95, SNORD96A, RNU6-2, using the following equation:
    ΔCt = (Ct_Target − Ct_housekeeping)
  3. For fold change calculations, calculate ΔΔCt values using a specific sample as the control, using the relative expression equation12:
    miRNA relative expression:
    2-∆∆Ct, where - ∆∆Ct = -[∆Ct test - ∆Ct control]
    NOTE: A fold change of 200 is considered as cutoff.
  4. Export fold change expression values to perform statistical analysis on R using the limma package in Bioconductor8.
  5. Correct for multiple comparisons using the Benjamini-Hochberg method13.
    NOTE:  A copy of the R script is available at: http://psilveyra.github.io/silveyralab/. Datasets and analyzed data are also available at the Gene Expression Omnibus under number GSE111667, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE111667.

7. Data Analysis: Functional Analysis Software

  1. Arrange the dataset. Include miRNAs with their respective expression log ratios and p-values. See Table 1 for specific dataset formatting.
  2. Open functional analysis software (version 01-10).
  3. Upload the dataset using the following format: File format: “Flexible Format”, Contains Column Header: “Yes”, Select Identifier type: “miRBase (mature)”, Array platform used for experiments: Choose the array platform.
    NOTE: Accepted files formats are .txt (tab delimited text files), .xls (excel files), and .diff (cuffdiff files).
  4. Select “Infer Observations” and verify that the experimental group labeling is correct.
  5. Go to “Dataset Summary” to revise the total amount of mapped and unmapped miRNAs.
  6. Click on the “new” button at the top left of program. Select “New MicroRNA Target Filter” and upload a microRNA dataset.
    1. Set Source to: TarBase, Ingenuity Experts Findings, miRecords.
    2. Set Confidence to: Experimentally observed or high (predicted).
    3. Select “Add Columns” to include a variety of biological information about the targets such as species, diseases, tissues, pathways and more.
    4. To focus on targets that have changed expression in the experiment, select “Add/ replace mRNA dataset”. Use “Expression Pairing” to locate microRNAs with the same or different expression levels.
      NOTE: The filter analysis will provide microRNA names and symbols, mRNA targets, the source that describes the target relationship, and the confidence level of the predicted relationship (Figure 2).
  7. Click “Add to My Pathway” to send the filtered dataset to a pathway canvas and explore more biological relationships.
  8. Use Path Designer to create a publication-quality model of microRNA effects.
  9. An alternate option is to create a core analysis.
    1. For core analysis type, select “Expression Analysis”.
    2. For measurement type, select “Expr Log Ratio”.
    3. Analysis Filter Summary: consider only molecules and/or relationships where: (species = mouse) AND (confidence = experimentally observed) AND (data sources = “Ingenuity Expert Findings”, “Ingenuity ExpertAssist Findings”, “miRecords”, “TarBase”, OR “TargetScan Human”).
    4. Select a cutoff of p-value = 0.05.
    5. Run the analysis.
      NOTE: The report will include: canonical pathways, upstream regulators analysis, diseases and functions, regulator effects, networks, molecules and more.

Results

The different cell types observed in smears are used to identify the mouse estrous cycle stage (Figure 1). These are identified by cell morphology. During proestrus, cells are almost exclusively clusters of round-shaped, well-formed nucleated epithelial cells (Figure 1A). When the mouse is in the estrus stage, cells are cornified squamous epithelial cells, present in densely packed clusters (Figure 1B

Discussion

MicroRNA profiling is an advantageous technique for both disease diagnosis and mechanistic research. In this manuscript, we defined a protocol to evaluate the expression of miRNAs that are predicted to regulate inflammatory genes in the lungs of female mice exposed to ozone in different estrous cycle stages. Methods for the determination of the estrous cycle, such as the visual detection method, have been described16. However, these rely on one-time measurements, and therefore are unreliable. To a...

Disclosures

The authors declare that they have no competing interests.

Acknowledgements

This research was supported by grants from NIH K01HL133520 (PS) and K12HD055882 (PS). The authors thank Dr. Joanna Floros for the assistance with ozone exposure experiments.

Materials

NameCompanyCatalog NumberComments
C57BL/6J miceThe Jackson Laboratory0006648 weeks old
UltraPure WaterThermo Fisher Scientific10813012
Sterile plastic pipetteFisher Scientific13-711-25Capacity: 1.7mL
Frosted Microscope SlidesThermo Fisher Scientific2951TS
Light microscopeMicroscope WorldMW3-H510X and 20X objective
Ketathesia- Ketamine HCl Injection USPHenry Schein Animal Health5585390 mg/kg. Controlled drug.
Xylazine Sterile SolutionLloyd Laboratories139-23610mg/kg. Controlled Drug.
EthanolFisher ScientificBP2818100Dilute to 70% ethanol with water.
21G gauge needleBD Biosciences305165
SyringeFisher Scientific3296541mL
Operating ScissorsWorld Precision Instruments501221, 50461314cm, Sharp/Blunt, Curved and 9 cm, Straight, Fine Sharp Tip
Tweezer KitWorld Precision Instruments504616
-80 ˚C freezerForma7240
Spectrum Bessman Tissue PulverizersFisher Scientific08-418-1Capacity: 10 to 50mg
RNase-free Microfuge TubesThermo Fisher ScientificAM124001.5 mL
TRIzol ReagentThermo Fisher Scientific15596026
Direct-zol RNA MiniPrep PlusZymo ResearchR2071
NanoDropThermo Fisher ScientificND-ONE-W
miScript II RT kitQiagen218161
Mouse Inflammatory Response & Autoimmunity miRNA PCR ArrayQiagenMIMM-105Z
Thin-walled, DNase-free, RNase-free PCR tubesThermo Fisher ScientificAM12225for 20 μl reactions
miRNeasy Serum/Plasma Spike-in ControlQiagen219610
Microsoft ExcelMicrosoft Corporationhttps://office.microsoft.com/excel/
Ingenuity Pathway AnalysisQiagenhttps://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/
R SoftwareThe R Foundationhttps://www.r-project.org/
Thermal cycler or chilling/heating blockGeneral Lab Supplier
MicrocentrifugeGeneral Lab Supplier
Real-time PCR cyclerGeneral Lab Supplier
Multichannel pipettorGeneral Lab Supplier
RNA wash bufferZymo ResearchR1003-3-4848 mL
DNA digestion bufferZymo ResearchE1010-1-44 mL
RNA pre-wash bufferZymo ResearchR1020-2-2525 mL
Ultraviolet ozone analyzerTeledyne APIModel T400http://www.teledyne-api.com/products/oxygen-compound-instruments/t400
Mass flow controllersSierra Instruments IncFlobox 951/954http://www.sierrainstruments.com/products/954p.html

References

  1. Rebane, A., Akdis, C. A. MicroRNAs: Essential players in the regulation of inflammation. Journal of Allergy and Clinical Immunology. 132 (1), 15-26 (2013).
  2. Cannell, I. G., Kong, Y. W., Bushell, M. How do microRNAs regulate gene expression?. Biochemical Society Transactions. 36 (Pt 6), 1224-1231 (2008).
  3. Pritchard, C. C., Cheng, H. H., Tewari, M. MicroRNA profiling: approaches and considerations. Nature Reviews Genetics. 13 (5), 358-369 (2012).
  4. Mi, S., Zhang, J., Zhang, W., Huang, R. S. Circulating microRNAs as biomarkers for inflammatory diseases. Microrna. 2 (1), 63-71 (2013).
  5. Sessa, R., Hata, A. Role of microRNAs in lung development and pulmonary diseases. Pulmonary Circulation. 3 (2), 315-328 (2013).
  6. Fuentes, N., Roy, A., Mishra, V., Cabello, N., Silveyra, P. Sex-specific microRNA expression networks in an acute mouse model of ozone-induced lung inflammation. Biology of Sex Differences. 9 (1), 18 (2018).
  7. Aris, R. M., et al. Ozone-induced airway inflammation in human subjects as determined by airway lavage and biopsy. American Review of Respiratory Disease. 148 (5), 1363-1372 (1993).
  8. Ritchie, M. E., et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Research. 43 (7), e47 (2015).
  9. Krämer, A., Green, J., Pollard, J., Tugendreich, S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics. 30 (4), 523-530 (2014).
  10. Machholz, E., Mulder, G., Ruiz, C., Corning, B. F., Pritchett-Corning, K. R. Manual restraint and common compound administration routes in mice and rats. Journal of Visualized Experiments. (67), (2012).
  11. Umstead, T. M., Phelps, D. S., Wang, G., Floros, J., Tarkington, B. K. In vitro exposure of proteins to ozone. Toxicology Mechanisms and Methods. 12 (1), 1-16 (2002).
  12. Livak, K. J., Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 25 (4), 402-408 (2001).
  13. Phipson, B., Lee, S., Majewski, I. J., Alexander, W. S., Smyth, G. K. Robust hyperparameter estimation protects against hypervariable genes and improves power to detect differential expression. Annals of Applied Statistics. 10 (2), 946-963 (2016).
  14. Smyth, G. K., et al. Linear Models for Microarray and RNA-Seq Data User's Guide. Bioconductor. , (2002).
  15. Benjamini, Y., Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society. Series B (Methodological). 57 (1), 289-300 (1995).
  16. Byers, S. L., Wiles, M. V., Dunn, S. L., Taft, R. A. Mouse estrous cycle identification tool and images. Public Library of Science ONE. 7 (4), e35538 (2012).
  17. Alves, M. G., et al. Comparison of RNA Extraction Methods for Molecular Analysis of Oral Cytology. Acta Stomatologica Croatica. 50 (2), 108-115 (2016).
  18. Wilfinger, W. W., Mackey, K., Chomczynski, P. Effect of pH and ionic strength on the spectrophotometric assessment of nucleic acid purity. Biotechniques. 22 (3), 478-481 (1997).
  19. Bustin, S. A., et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clinical Chemistry. 55 (4), 611-622 (2009).
  20. Walker, S. E., Lorsch, J. RNA purification- precipitation methods. Methods in Enzymology. 530, 337-343 (2013).
  21. Git, A., et al. Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression. RNA. 16 (5), 991-1006 (2010).
  22. Smyth, G. K. Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions Using R and Bioconductor. , 397-420 (2005).
  23. Griffiths-Jones, S. miRBase: the microRNA sequence database. Methods Mol Biol. 342, 129-138 (2006).
  24. Sethupathy, P., Corda, B., Hatzigeorgiou, A. TarBase: A comprehensive database of experimentally supported animal microRNA targets. RNA. 12 (2), 192-197 (2006).
  25. Agarwal, V., Bell, G. W., Nam, J., Bartel, D. P. Predicting effective microRNA target sites in mammalian mRNAs. eLife. 4, e05005 (2015).
  26. Xiao, F., Zuo, Z., Cai, G., Kang, S., Gao, X., Li, T. miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res. 37, D105-D110 (2009).
  27. Mullany, L. E., Wolff, R. K., Slattery, M. L. Effectiveness and Usability of Bioinformatics Tools to Analyze Pathways Associated with miRNA Expression. Cancer Informatics. 14, 121-130 (2015).

Reprints and Permissions

Request permission to reuse the text or figures of this JoVE article

Request Permission

Explore More Articles

Lung MicroRNAEstrous CycleOzone ExposureFemale MiceInflammatory GenesHormonal ContributionsTissue HarvestingRNA ExtractionDNase TreatmentCentrifugationHomogenizationExperimental ProtocolSample Collection Methods

This article has been published

Video Coming Soon

JoVE Logo

Privacy

Terms of Use

Policies

Research

Education

ABOUT JoVE

Copyright © 2025 MyJoVE Corporation. All rights reserved