|Home | About | Journals | Submit | Contact Us | Français|
Dicer is a ribonuclease whose major role is to generate mature microRNAs although additional functions have been proposed. Deletion of Dicer leads to embryonic lethality in mice. To study the role of Dicer in adults, we generated mice in which administration of tamoxifen induces deletion of Dicer. Surprisingly, disruption of Dicer in adult mice induced lipid accumulation in the small intestine. To dissect the underlying mechanisms, we carried out miRNA, mRNA and proteomic profiling of small intestine. The proteomic analysis was done using mice metabolically labeled with heavy lysine (SILAC mice) for an in vivo readout. We identified 646 proteins of which 80 were upregulated >2-fold and 75 were downregulated. Consistent with the accumulation of lipids, Dicer disruption caused a marked decrease of microsomal triglyceride transfer protein, long-chain fatty acyl-CoA ligase 5, fatty acid binding protein, and very-long-chain fatty acyl-CoA dehydrogenase, among others. We validated these results using multiple reaction monitoring (MRM) experiments by targeting proteotypic peptides. Our data reveal a previously unappreciated role of Dicer in lipid metabolism. These studies demonstrate a systems biology approach by integrating mouse models, metabolic labeling, gene expression profiling and quantitative proteomics can be a powerful tool for understanding complex biological systems.
Small non-coding RNA-mediated control of gene-expression has been demonstrated in a number of biological pathways. To date, three major classes of small non-coding RNAs—microRNAs (miRNAs), small interfering RNAs (siRNAs) and Piwi-interacting RNAs—have been widely studied. The biogenesis of miRNAs and siRNAs involves precursor molecule processing by the ribonuclease III enzyme Dicer.1, 2 Dicer, a ubiquitously expressed and evolutionarily conserved protein, consists of a number of domains including DRBM, RNase III, PAZ, helicase and a DECH-box motif. To generate mature miRNAs, the PAZ domain of Dicer binds the 2-nucleotide 3’ overhang of precursor miRNAs while the two RNase III domains excise the loop portion of precursor miRNAs.3 Dicer associates with TRBP, another RNA binding protein, to load mature miRNAs onto the RNA-induced silencing complex, which mediates base-pairing between miRNAs and their cognate mRNA targets.
A variety of biological systems have been used to investigate the many functions of Dicer. Overexpression of Dicer has been shown to promote cell growth in interferon-defective Li-Fraumeni syndrome fibroblasts.4 Conversely, global deletion of Dicer in mice leads to lethality at early stages of embryonic development.5 Bernstein et al. deleted exon 22 of Dicer1 in mice to demonstrate that homozygous Dicer1−/− mouse embryos died as early as embryonic day 7.5. Wienholds et al. corroborated these findings in zebrafish by introducing nonsense mutations in the RNase III domain of Dicer.6 Growth of homozygous dicer1−/− fish embryos arrests on day 10 post-fertilization and fish die three weeks post-fertilization. To bypass this embryonic lethality caused by global deletion of Dicer, tissue-specific promoters for Cre recombinase were used to disrupt Dicer expression in different organs.7–9 Conditional ablation of Dicer increases apoptosis of myoblasts with a subsequent decrease in muscle mass.10 Liver-specific Dicer deficiency results in lipid accumulation and glycogen depletion in hepatocytes as well as autochthonous development of liver cancers.11 Dicer deletion has also been found to affect the complete differentiation of T and B lymphocytes.12, 13 Additionally, mouse embryonic stem cells that lack Dicer fail to properly silence their centromeric chromosomal region.14 A recent study in Caenorhabditis elegans also described a previously unexplored role of Dicer in DNA fragmentation during apoptosis.15 Mechanisms underlying most of the Dicer knockout phenotypes remain to be fully elucidated.
In instances where biological effects or pathways may be difficult to dissect, novel and integrated approaches are often the best way to gain more insight. For example, mass spectrometry-based quantification strategies, including label-free or in vitro labeling methods and in vivo metabolic labeling methods,16 represent powerful tools to address major biological questions. Indeed, many studies have adopted a label-free quantification approach owing to its flexibility in terms of multiplexing and cost-effectiveness.17, 18 Spectral counting and peak area ratios are among the commonly used features to quantify peptides in the label-free strategy.19, 20 However, label-free quantification is dependent on highly reproducible sample handling and LC conditions between samples—which is generally not practical for analyzing a large number of samples—and small changes between samples cannot be reliably quantified.21 Another commonly used in vitro labeling quantification method is isobaric tagging for relative and absolute quantification (iTRAQ).22 Although iTRAQ offers many advantages over the label-free approach, efficiency of labeling, co-elution of physicochemically similar peptides, and maximum multiplexing of 8 samples still limit its application. Further, iTRAQ quantification is generally not compatible with ion-trap mass analyzers. Thus, in vivo labeling offers the advantage of measuring changes by minimizing sample processing errors. Since 15N labeling generates complex spectra, Kruger et al. circumvented this problem by using 13C6-lysine metabolic labeling (SILAC) of mice to reveal the importance of kindlin-3 in structural integrity of red blood cells.23
To investigate the effects of Dicer deletion, we generated conditional knockout mice with a tamoxifen-regulated Cre-ERT2 system. Our studies reveal that global disruption of the Dicer gene in adult mice induces multiple changes that include abnormalities of bone marrow and small intestine. Here, we systematically evaluated global proteomic changes in the small intestine in the absence of functional Dicer. For in vivo quantification, we employed 13C6-lysine SILAC labeling of mice and found that Dicer deletion causes dysregulation of lipid metabolism, among other abnormalities. Further validation using multiple reactions monitoring (MRM) corroborated our findings in SILAC mice. Overall, our study demonstrates the strength of quantitative proteomic approaches such as SILAC in animal systems to study gene functions. Our findings indicate a previously unappreciated role for Dicer in lipid metabolism that can now be studied in greater detail.
C57BL/6 transgenic mice were maintained in a pathogen-free environment. ROSA26-CreERT2 mice were crossed to Dicer1lox/lox mice to generate tamoxifen inducible deletion of exons 21 and 22. Briefly, 8-week-old female mice were administered tamoxifen (1 mg/day gavage) or vehicle control for five days. Dicer1lox/lox mice without ROSA26-CreERT2 were also administered tamoxifen and used as controls. Tissues were harvested after a 3-hour fasting period on day 8. Mice were perfused with 10% formalin prior to necropsy. For Oil Red O staining, the jejunum was harvested without formalin fixation and immediately immersed in O.C.T. Compound (Tissue-Tek) on dry ice. Tissue sections were rinsed with 70% ethanol, placed in Oil Red O stain for 10 minutes, washed twice in 70% ethanol and rinsed in water. Harris Hematoxylin was used for counterstaining. All animal studies were approved by the Institutional Animal Care and Use Committee of the Johns Hopkins University School of Medicine.
qRT-PCR was used to confirm deletion of Dicer1 exons 21 and 22 and to validate microarray results. Total RNA was purified using RNeasy Mini Kit (Qiagen, Valencia, CA) after the integrity of RNA was checked (2100 Bioanalyzer, Agilent Technology, Santa Clara, CA). Total RNA was reverse-transcribed using SuperScript III Reverse Transcriptase with oligo(dT)20 primers (Invitrogen, Carlsbad, CA) for confirming Dicer1 deletion while RNA to cDNA EcoDry Premix Oligo dT primers (Clontech, Mountain View, CA) was used for validating the gene expression microarray results. Biological triplicate samples were analyzed using iQ™SYBR Green Supermix Assay (Bio-Rad, Hercules, CA). All primer sequences are listed in Supplementary Table 1. To validate miRNA microarray results, we analyzed biological triplicate samples employing TaqMan MicroRNA Assay (Life Technologies, Grand Island, NY). The relative abundance of transcripts was calculated using Cq-Log10 standard curve method with β-actin or snoRNA135 as controls.
Total RNA was labeled and hybridized with Affymetrix GeneChip Mouse Exon 1.0 ST array for gene expression analysis as the manufacturer's instructions. Total RNA was labeled and hybridized with Affymetrix GeneChip miRNA array for miRNA analysis. All experiments were performed as biological triplicates. Analysis was carried out in GeneSpring GX 10 (Agilent, Santa Clara, CA). Benjamini-Hochberg method built in GeneSpring GX 10 was used for the multiple testing correction.
Stable isotope-labeled mouse food was obtained from Cambridge Isotope Laboratories (Mouse Feed Labeling Kit, catalog number: MLK-LYS-C). A twenty-eight-day-old female mouse (F0) was fed the diet containing 13C6-lysine and mated 14 days later. Heavy diet was continued until F1 pups were weaned. One F1 female mouse was maintained on heavy diet and propagated to F2. A littermate female mouse (F0) was fed light diet to generate unlabeled control mice. The labeling efficiency was monitored by mass spectrometric analysis of blood, liver and lung specimens collected at 10 weeks of age from F0 and 4 weeks of age from the F1 generation.
Small intestine was harvested and homogenized in 2% SDS buffer using POLYTRON PT 1200 E homogenizer (KINEMATICA, Bohemia, NY) and sonicated. Proteins from control (heavy) and Dicer-deleted (light) tissues were mixed in equal proportions. In-solution digestion was modified from the previously described.24 Briefly, disulfide bonds were reduced with 5mM dithiothreitol (60 °C, 45 min) and alkylated with 20 mM iodoacetamide. Amicon Ultra 30 kDa centrifugal filters (Milipore, Billerica, MA) were used for changing the lysate buffer to 8 M urea solution. Proteins were digested using Lys-C (1:100 (w/w), Wako Chemical USA, Richmond, VA) at 37 °C for 4 hours. Following a 4-fold dilution, the samples were incubated with trypsin (1:20 (w/w), modified sequencing grade, Promega, Madison, WI) at 37 °C for 16 hours.
SCX fractionation was modified from the previously described.24 Briefly, in-solution digests were dried and reconstituted in SCX solvent A (10 mM potassium phosphate buffer in 20% acetonitrile, pH 2.85). Peptides were loaded onto a POLYSULFOETHYL A column (100 × 4.6 mm, 300Å, 5µm, PolyLC, Columbia, MD) and resolved with an increasing gradient of SCX solvent B (solvent A containing 350 mM KCl, pH 2.85) over 70 minutes. The collected fractions were dried and desalted on a C18 column. For preparation of samples used in MRM analysis, basic RPLC fractionation was carried out on the XBridge C18 column (250 × 4.6 mm, 300Å, 5µm, Waters, Milford, MA) with mobile phase solvent A of 10 mM triethylammonium bicarbonate (TEABC, pH 9.5) and 10 mM TEABC with 90% acetonitrile (pH 9.5) as the mobile phase B increased from 5%, to 100% in 30 minutes, and persisted for 10 minutes.
LC-MS/MS analysis of day-8 samples was carried out on a quadrupole time-of-flight mass spectrometer (Agilent's 6538 Ultra High Definition Accurate-Mass Q-TOF) equipped with an HPLC-chip cube: separation-150 mm × 75 µm, enrichment-9 mm 160 nl, Zorbax 300SB-C18 5 µm. The Q-TOF was operated at capillary voltage of 1800 V, fragmenter voltage of 150 V, drying gas temperature of 300°C, drying gas flow rate of 6 l/min, medium isolation width of 4 m/z and a collision energy slope of 3 V plus offset of 2 V. In each duty cycle, 6 precursors were chosen, and spectra were acquired in the range of m/z 350–1700. Precursors with a single charge or unknown charge state were excluded. Database search was carried out using Spectrum Mill MS Proteomics Workbench (Agilent Technologies, Santa Clara, CA) against the mouse RefSeq database. The search parameters were as follows: trypsin and Lys-C as proteolytic enzymes, up to one missed cleavage, precursor peptide mass range 500 ~ 8000 Da, peptide mass error tolerance of 20 ppm, fragment mass error tolerance of 100 ppm, carbamidomethylation of cysteine and SILAC 2 (Lys 0-6Da) as a fixed modification, oxidation of methionine, acetylation of protein N-termini and deamidation of glutamine and asparagine as variable modifications. The forward and reverse scores provided in .spo files generated by Spectrum Mill were used to calculate the interpolated false discovery fate (FDR). FDR of 1% was used to select proteins for further analysis.
Small intestine tissue from four Dicer knockout mice and four controls were lysed in 8 M urea solution. Lysates were digested with Lys-C/trypsin and fractionated by basic RPLC. Fractions were analyzed separately on Agilent 6538 Ultra High Definition Accurate-Mass Q-TOF LC/MS. Proteins identified from the SILAC mice were selected as the candidates for MRM-based validation. Protein sequences of seven differentially regulated proteins in Dicer knockout mice were used as the input in Skyline (version 0.7.0.2556) to select the proteotypic peptides. Tryptic peptides with complete digestion and absence of cysteine and methionine residues were given preference for selection. For peptides containing histidine residues, +2 charged fragment ions were also included in the transition list. The majority of the selected precursors were with charge +2. At least two peptides with a minimum of 5 transitions were selected for each of the seven proteins.
For MS-based validation of differentially expressed proteins in the small intestine of Dicer knockout mice, selected peptides from seven proteins were monitored on Agilent's 6430 triple quadrupole mass spectrometer. MassHunter Optimizer software (version B.04.00, Agilent, Santa Clara, CA) was used for optimization of dwell time and collision energy to specify transitions of selected peptides. Titration curve of intensity was plotted for increasing amount of the standard heavy peptides. Sample was loaded and resolved onto ZORBAX 300SB-C18 chip using 1200 series HPLC. The triple quadrupole mass spectrometer was also interfaced with the HPLC-chip cube. A 35 min gradient with increasing solvent B (3%–100%) was applied to resolve peptides. Data were acquired using MassHunter, and the relative abundance was calculated by comparing the total intensity of 5 transitions.
Dicer is a ubiquitously expressed protein involved in the biogenesis of miRNAs. However, Dicer also plays an important role in transcriptional silencing and in the formation of heterochromatin.25, 26 Here, we employed metabolic labeling in the context of tamoxifen-inducible disruption of the Dicer gene to systematically investigate the function of Dicer in adult mammalian tissues.
To overcome embryonic lethality resulting from global disruption of the Dicer gene during development, we used the Cre-loxP system to disrupt the expression of floxed Dicer alleles in adult mice.27 For this purpose, tamoxifen responsive ROSA26-driven Cre recombinase (Cre-ERT2) was used to specifically delete floxed exons 21 and 22 of the Dicer gene (Figure 1). In the absence of tamoxifen, double-transgenic adult mice (CRE-ERT2:Dicer1lox/lox) were phenotypically normal and fertile. Tamoxifen administration facilitates the nuclear import of Cre recombinase, thereby triggering excision of floxed Dicer1 exons 21 and 22. qRT-PCR analysis confirmed the disruption of Dicer based on monitoring of the junctional region spanning exons 21 and 22. Excision of the Dicer floxed allele was ~91%. Thus, the inducible knockout of Dicer in adult animals enabled us to address its homeostatic roles in fully differentiated tissues.
Upon homozygous deletion of Dicer1 in Cre-ERT2:Dicer1lox/lox mice, we expected a poor survival rate and severe phenotypic changes. We monitored the mice daily and no gross physical abnormalities were observed in the first 6 days. Beginning day 7, half of the tamoxifen-treated mice developed diarrhea and 80% (4/5) of these mice died at day 10 post Dicer deletion. Histopathological examination of day-8 Dicer-deleted mice revealed no overt abnormalities in the lung, heart, liver, kidney or muscle. In contrast, the small intestinal epithelium showed villous distortion, abnormal cytoplasmic vacuolation, and enlarged nucleoplasm with prominent nucleoli in enterocytes (Figure 2, A, B). Bone marrow failure was also observed with a reduction in myeloid cells (Figure S1). Further, Dicer knockout mice showed evidence of splenic and thymic cortical atrophy. Given the striking pathology of the small intestine, we decided to examine in greater detail the role Dicer plays in maintaining intestinal tissue homeostasis.
Upon histological examination, extensive and abnormal vacuolation of enterocytes was evident in the small intestine of Dicer knockout but not control mice. Although vacuolation of enterocytes could result from a variety of reasons, we hypothesized that it likely results from aberrant lipid accumulation. To confirm this, we used Oil Red O staining to detect the presence of neutral lipids. As seen in panels C and D in Figure 2, the vacuolated areas were prominently stained by Oil Red O, indicating abnormal accumulation of lipid droplets in the enterocytes. A similar phenotype has been documented in humans with abetalipoproteinemia, where lipid droplets accumulate abnormally in the small intestinal epithelium. In contrast, the small intestine of control mice harbored little, if any, neutral lipid droplets. These droplets could contain any of the neutral lipids such as triglycerides, free fatty acids and cholesterol esters. Despite the aberrant lipid droplet accumulation observed in the enterocytes, the specific defect in lipid metabolism that was altered in the Dicer knockout mice was not apparent (i.e., whether the defect was in lipid absorption, catabolism, packaging and/or transport). To gain mechanistic insights into intestinal pathology of the Dicer knockout mouse model, we combined inducible Dicer deletion and SILAC mice to measure global proteomic changes resulting from disruption of the Dicer gene. This in vivo labeling strategy also enabled us to study the microenvironment of the small intestine, which cannot be recapitulated by in vitro studies.
Because Dicer is responsible for generating mature miRNAs, we predicted that Dicer knockout can deplete the whole repertoire of miRNAs. Thus, we carried out microRNA microarray profiling (as biological triplicates) and identified 150 miRNAs with the p-value cutoff of 0.05, expressed in the small intestine of Dicer knockout and control mice (GEO Accession #: GSE34909). Interestingly, only 24 miRNAs were downregulated ≥2-fold on day 8 after Dicer deletion. The most downregulated miRNAs were mmu-miR-215 (75.5-fold) and mmu-miR-194 (18.2-fold). Quite unexpectedly, we also observed 18 miRNAs to be upregulated ≥2-fold, including mmu-miR-195 (3.8-fold of control) and mmu-miR-199b (4.0-fold of control) (Table 1). qRT-PCR experiments carried out to validate the changes in a subset of miRNA microarray results revealed the same trend as microarray results (Figure S2).
Previous studies to examine differential regulation of protein abundance by miRNAs have shown that the magnitude of the effect observed is often subtle. Although both in vitro labeling as well as label-free approaches can provide quantitative data, they are not generally as reliable as in vivo labeling approaches such as SILAC if the changes are not dramatic. Therefore, we decided to use SILAC mice for our analysis. In this approach, proteins are labeled in vivo by feeding mice a 13C6-lysine diet. We observed incorporation of 13C6-lysine within two generations of mice fed continuously with the 13C6-lysine diet as assessed by mass spectrometric analysis of proteins extracted from their blood, liver, and lung. After 10 weeks on this diet, the labeling efficiency of proteins in F0 mice (as indicated by percent incorporation in proteins) ranged from ~88% to 99%. After a 4 week 13C6-lysine diet through the F1 generation, the labeling efficiency of proteins was >97%, indicating near complete incorporation of exogenously fed 13C6-lysine (Figure 3A and Supplementary Table 2).
SILAC mice were treated with vehicle control while unlabeled mice were administered tamoxifen to induce deletion of Dicer. Small intestine tissue from the two groups of mice was harvested on day 8 and equal amounts of tissue lysates were mixed from SILAC and light mice. The pooled lysates were digested with Lys-C followed by trypsin. The digests were subsequently separated by SCX to generate 24 fractions. These fractions were analyzed by LC-MS/MS (Figure 3B).
The data were searched against a mouse RefSeq database and a reversed protein sequence database using Spectrum Mill. Peptides scoring better than the 1% FDR cut-off were considered for further analysis. In all, we identified 3,087 peptides corresponding to 646 protein groups (Supplementary Table 3). Because only lysine containing peptides could be quantified, we obtained fold changes for ~54% of the identified peptides (derived from 313 proteins) (Figure 4A).
As Dicer generates mature miRNAs which suppress translation of mRNAs, we expected Dicer knockout to exhibit a widespread increase of protein abundance. To our surprise, we did not observed any major change in abundance in half of the proteins. We found 80 proteins to be upregulated (≥2 fold) and 75 proteins to be downregulated (≥2 fold) in the small intestine (Figure 4A). Representative lists of upregulated and downregulated proteins are provided in Table 2. Strikingly, 34 out of 80 upregulated proteins were ribosomal protein subunits. Additionally, proteins relevant to ribosomal biogenesis such as ribosomal RNA 2’-O-methyltransferase, nucleolin and nucleophosmin were also upregulated. Despite synchronized upregulation of proteins involved in ribosomal biogenesis and protein constituents of the translation machinery (Table 2), processes expected to increase overall protein synthesis, we still observed downregulation of 75 proteins in our study, suggesting production of these 75 proteins was significantly compromised. As expected from the abnormal lipid phenotype of Dicer knockout mice, enzymes that are involved in lipid metabolism appeared to be one of the prominent classes of proteins that were among these 75 proteins.
Considering deletion of Dicer causes a measurable change in the profile of miRNAs (see above) that interact with mRNAs and regulate protein abundance, we also performed transcriptomic analysis using DNA microarrays. Expression levels of 16,526 genes were quantifiable between Dicer deleted and control mice. After applying a p-value cutoff of 0.05 and correcting for multiple hypothesis testing, a total of 1,689 genes were found to exhibit >2-fold change in mRNA expression (GEO Accession #: GSE34905). We validated the changes at the mRNA level on a subset of transcripts using qRT-PCR, which confirmed the same trend as was observed in the microarray experiments (Figure S3). Overlaying the data from transcriptomic and proteomic analyses revealed that 219 genes have both mRNA and protein quantification values for direct comparison. Of these, ~76% exhibited a change in mRNA and protein levels in the same direction (Figure 4B, the correlation coefficient=0.68) while almost 15% genes showed a decrease in mRNA but an increase in protein abundance. Conversely, ~10% of the genes showed an increase in mRNA but a decrease in protein abundance. We did not observe any significant correlation between miRNA and their mRNA target levels and between miRNA and their protein target levels (the correlation coefficients are both <0.10). These data reveal a complex regulation of mRNA and protein levels by alterations in miRNA levels and presumable other functions of Dicer.
When we examined the repertoire of proteins whose abundance was altered in the small intestine of Dicer knockout mice, we used Ingenuity Pathway Analysis (IPA; http://www.ingenuity.com/products/pathways_analysis.html) and found that proteins involved directly or indirectly in lipid metabolism constituted a major fraction (16 out of 75 ≥2-fold, 14 out of 47 ≥4-fold) of the downregulated proteins observed in the small intestine of Dicer knockout mice. The list included proteins involved in lipid transport and catabolism (Table 3 and Supplementary Table 4). In abetalipoproteinemia patients, microsomal triglyceride transfer protein gene (MTTP) is mutated leading to abnormal lipid accumulation in enterocytes, a phenotype resembling that of the Dicer knockout mice. We found that microsomal triglyceride transfer protein (Mttp) was downregulated 11.4-fold in the small intestine of Dicer gene-deleted mice (Figure 5, A, B). The major function of Mttp is to package triglycerides into pre-chylomicron that is exported out of enterocytes into the lymphatic system and transported to peripheral tissues (e.g., adipose tissue, skeletal muscle, and liver) for storage or catabolism. Dihydroxyacetone kinase (Dak) was downregulated 32.3-fold in the small intestine of the Dicer knockout mice. Dak phosphorylates dihydroxyacetone, to generate dihydroxyacetone phosphate (DHAP), an intermediate metabolite in glycolysis. Downregulation of Dak can lead to the accumulation of dihydroxyacetone, which in turn may prevent glycerol (a precursor for triglyceride synthesis) from entering the glycolytic pathway, potentially channeling the glycerol substrate toward the triglyceride synthesis pathway.
Other downregulated enzymes included alcohol dehydrogenase 6A (Adh6a, 16.1-fold) and aldehyde dehydrogenase (Aldh2, 3.0-fold), fatty acid binding protein-1 (Fabp1, NP_059095.1, 31.3-fold), long-chain fatty-acyl-CoA ligase 5 (Acsl5, NP_082252.1, 4.6-fold), and very-long-chain fatty-acyl-CoA dehydrogenase (Acadvl, NP_059062.1, 1.7-fold). Fabp1 is a lipid chaperone, playing an important role in facilitating the entry of long-chain fatty acids into smooth endoplasmic reticulum. In Dicer knockout mice, Fabp1 protein levels reduced 31-fold compared to control animals (Figure 5, C, D). Acadvl plays an important role in mitochondria β-oxidation of long-chain fatty acids. Reducing Acadvl protein levels promotes accumulation of long-chain fatty acids. Overall, our quantitative proteomic analysis suggests that proteins involved in lipid metabolism are significantly altered in the small intestine of Dicer knockout mice, providing a possible mechanism accounting for abnormal lipid droplet accumulation in enterocytes of these animals (Figure 6). Our global transcriptomic analysis also revealed a general concordance between changes in mRNA level and protein levels for these molecules involved in metabolism (Figure 4B, the correlation coefficient=0.68). Reduction in the protein levels of Mttp, Acsl5, Fabp1 and Dak was mirrored by the downregulation of their corresponding mRNA levels (9.4-fold, 3.8-fold, 49.9-fold and 5.5-fold, respectively). In addition, we have also identified two downregulated proteins normally involved in glucose metabolism—fructose bisphosphate aldolase B (Aldob, NP_659152.1, 47.6-fold) and fructose bisphosphatase 1 (Fbp1, NP_062268.1, 18.2-fold).
Because changes in miRNAs can affect mRNA and protein abundance, we sought to correlate these parameters for the above-mentioned proteins involved in lipid metabolism. TargetScanMouse (Release 6.0) analysis of potential targets of upregulated miRNAs highlighted mmu-miR-23a and -23b, which were upregulated 1.8- and 1.2-fold, respectively, in Dicer knockout mice with their potential cognate targets being Aldob and Dak. mRNA levels of these two proteins were decreased 10- and 5.6-fold, respectively, in the transcriptomic analysis. Given their relatively long 3’ untranslated regions, Aldob (1086 nt) and Dak (1580 nt) can be targeted by multiple miRNAs including mmu-miR-23a and -23b, both of which were upregulated in the small intestine of Dicer knockout mice (Table 3).
To validate the altered levels of proteins identified from SILAC mice, we developed targeted MRM-based assays to quantify changes induced by Dicer gene deletion in a different cohort of unlabeled mice. We carried out biological experiments in quadruplicate for this purpose. Equal amounts of protein lysates were subjected to basic RPLC fractionation and analyzed on a triple quadrupole mass spectrometer. The list of all candidate peptides selected for MRM-based validation and the corresponding fold-change values obtained from MRM assays are provided in Supplementary Table 5. We monitored five transitions for each peptide. Of the downregulated proteins, we selected candidate peptides from Mttp and Dak for validation. We targeted three peptides from Mttp: SGFTTANQVLGVSSK (charge: +2, monoisotopic m/z: 748.3921) (Figure 7A), AFALNFQQTIAGK (charge: +2, monoisotopic m/z: 704.8828), and EFYSYENEPVGIENLK (charge: +2, monoisotopic m/z: 965.96118). The ratios of ion count between Dicer knockout and control for these three peptides were 0.09, 0.09, and 0.10; the ratio of this protein obtained from the SILAC mice was 0.09.
The proteotypic peptides targeted in Fabp1 were AIGLPEDLIQK (charge: +2, monoisotopic m/z: 598.8479) (Figure 7B) and NEFTLGEECELETMTGEK (charge: +2, monoisotopic m/z: 1058.9561). The ratios of ion counts between Dicer knockout and control were 0.04 and 0.03; the value obtained from SILAC mice result was 0.03. These two downregulated proteins in Dicer knockout mice are responsible for the packaging and transport for triglycerides. Among the upregulated proteins, we targeted two proteotypic peptides in nucleolin, TLVLSNLSYSATK (charge: +2, monoisotopic m/z: 698.8877) (Figure 7C) and GLSEDTTEETLK (charge: +2, monoisotopic m/z: 661.8197). Again, the results from MRM quantification corroborated results obtained from SILAC mice, which were 2.7 and 3.0 vs. 2.2. Overall, the quantification values from MRM validation and SILAC mice appeared robust and consistent, especially for those proteins with marked changes caused by Dicer gene deletion.
In this study, we have elucidated an important but previously unappreciated role of Dicer in maintaining small intestinal tissue homeostasis. Through inducible disruption of the Dicer gene in adult mice, we uncovered a Dicer-dependent mechanism regulating lipid transport and catabolism in enterocytes. In our mouse model, deletion of Dicer also reduced cellularity within the bone marrow niche. It is interesting that intestine and bone marrow exhibited pronounced phenotypes when Dicer was deleted because of high rates of proliferation and turnover that are common to both of these organs. Overall, our data clearly highlight the role of Dicer beyond embryonic development and reveal its critical role in maintaining tissue homeostasis in adult animals.
Our quantitative proteomic analysis provided mechanistic insights into metabolic processes regulated by Dicer, the absence of which resulted in abnormal lipid accumulation in enterocytes similar to abetalipoproteinemia in humans. Several enzymes involved in lipid transport (e.g. Mttp and Fabp1) and catabolism (e.g. Acadvl, Acsl5, Adh6a, Aldh2, and Dak) were strikingly downregulated. Downregulation of lipid-metabolizing enzymes in Dicer knockout mice occurred despite upregulation of ribosomal proteins and protein constituents of the translational machinery. In a simplified model, Dicer generates mature miRNAs, which suppress the translation of cognate mRNA targets or destabilize the transcripts. In this model, Dicer knockout is predicted to result in upregulation of proteins. However, we identified proteins that are downregulated as well as those that are upregulated. Since multiple functions of Dicer outside of microRNA biogenesis have been described, further studies are required to dissect additional roles of Dicer.
In addition to the discovery phase of our study, we developed MRM assays to target peptides identified in SILAC mice. Results from these MRM assays helped confirm that downregulation of lipid transport and metabolizing enzymes likely accounts for the abnormal lipid droplet accumulation in enterocytes of Dicer knockout mice. The use of these two quantitative methods in tandem provides a promising platform for discovering and confirming diverse roles of proteins in biological processes. Because of the regulatory role of miRNAs in protein abundance, we envisage that the approach presented here can be adopted for discovery and validation of miRNAs regulating proteins in metabolism research.
Figure S1. Bone marrow failure in mice lacking Dicer. The cellularity of bone marrow in Dicer knockout mice was drastically decreased, especially in the myeloid lineage, as compared to controls.
Figure S2. Validation of miRNA microarray experiments by qRT-PCR. TaqMan microRNA assays were used to quantitate miRNA expression levels for the indicated miRNAs. The red bars represent qRT-PCR results while the blue bars depict microarray results.
Figure S3. Validation of DNA microarray results by qRT-PCR. SYBR-Green assays were used to quantitate mRNA expression levels for the indicated transcripts. The green bars represent qRT-PCR results while the blue bars depict microarray results. A: genes observed to be upregulated in Dicer knockout mice; B: genes observed to be downregulated in Dicer knockout mice.
This study was supported in part by an NIH roadmap grant for Technology Centers of Networks and Pathways (U54RR020839) and a contract (HHSN268201000032C) from the National Heart, Lung and Blood Institute. G.W.W is supported by grants from the American Heart Association (SDG2260721) and the National Institute of Health (RO1 DK084171). J.M.P is supported by the NIH National Research Service Award (F32DK084607).