pathological tissue fibrosis is the abnormal accumulation of collagen - rich extracellular matrix ( ecm ) after a chronic or misregulated response to injury that progressively disrupts tissue architecture , leading to tissue stiffness , impaired organ function , and eventually organ failure . fibrosis is featured in diverse conditions , accounts for as much as 45% of all deaths worldwide , and appears to be increasing in prevalence . fibrosis complications arise in very different disease settings , from autoimmunity and environmentally induced inflammation to cancer , spanning multiple organs . to date , the treatment options are extremely limited for attenuating or reversing this process . thus , there is an urgent need to delineate underlying pathological mechanisms that may lead to new therapeutic approaches . both the initiation and persistence of pathological fibrosis involves the activation and differentiation of progenitors to myofibroblasts , the key effectors in fibrosis . myofibroblasts play an important role in executing physiologic tissue repair , leading to matrix deposition , wound contraction , and healing on one hand , and pathological fibrogenesis leading to chronic fibrosing conditions on the other . accordingly , prominent molecular pathways in acute tissue repair often recapitulate their fibrogenic function in chronic fibrotic conditions , essentially conditions in which wound healing has gone awry ; these include platelet derived growth factor receptor beta ( pdgfr ) and transforming growth factor beta ( tgf- ) . however , there are still many gaps in our understanding of the mechanisms underlying fibrogenesis . due to their morphology and function , these cells are incredibly difficult to study in  vitro . recently , elegant fate mapping studies pointed to a role for pericytes as a significant progenitor pool for myofibroblasts after injury in a variety of organ systems . liver pericytes , also known as hepatic stellate cells ( hscs ) , have been shown to be the major progenitor pool for myofibroblasts after carbon tetrachloride ( ccl4)-induced liver injury and fibrosis,5 , 6 and gene expression profiling of these cells isolated at various time points after ccl4 administration identified molecular alterations that might be functionally relevant to fibrogenesis . likewise , myofibroblasts and their precursors have been transcriptionally profiled during kidney fibrosis in mice using translational ribosome affinity purification technology . isolation of this cell type induces a lot of alteration in what these cells express ; therefore , it is important to study the function of the genes in their native disease environment without disrupting cell - cell and cell - matrix interactions . we aimed to identify novel modifiers of tissue fibrosis expressed in myofibroblasts or their progenitors through rna silencing in  vivo . although much more challenging than conducting a cell - based screen in  vitro , this in  vivo approach was pursued to enable the interrogation of gene function in the context of the complex tissue environment and thereby yield physiologically relevant fibrosis modifiers . in order to achieve this objective , we employed recently generated transcriptomes from myofibroblasts and their precursors in models of liver and kidney injury and fibrosis.6 , 7 to achieve effective gene silencing in  vivo , we used small interfering rna ( sirna ) delivery with lipid nanoparticles ( sirna - lnps ) , an emerging technology for gene silencing in  vivo , particularly in the liver , and hence deployed this technology in a model of liver fibrogenesis . the sirna - lnp technology has been previously used to show that silencing of the collagen type i alpha 1 gene ( col1a1 ) reduced its mrna levels and collagen accumulation in liver fibrosis models.9 , 10 however , the use of sirna - lnp technology in a targeted assay for identifying novel fibrosis modifiers has not yet been reported . herein , we report our use of a platform composed of novel sirna targets , sirna - lnp delivery technology , and ccl4-induced liver fibrosis , resulting in the identification of novel modifiers of fibrogenesis in  vivo . we aimed to silence genes in key fibrogenic cell lineages , myofibroblasts , and pericytes to identify novel mediators of tissue fibrosis . candidate genes were selected by the intersection of gene expression datasets for this cell lineage in two different organ systems . we identified genes that were commonly at least 2-fold upregulated in activated hscs isolated from livers of mice treated with ccl4 for 2  months and myofibroblasts and pericytes , their precursors , in the mouse kidney 2 and 5  days after unilateral ureteral obstruction ( uuo ) . we excluded transcripts that were not associated with a gene product , had no human homolog , or had well - known functions in fibrosis ( figure  1 ) . ultimately , we unbiasedly tested 24 genes by sirna - mediated gene knockdown ( kd ) in  vivo ( table s1 ) . ccl4-induced liver injury and fibrosis in mice is a well - characterized model consisting of repeated ccl4 administration that injures hepatocytes , followed by a fibrogenic reaction . hscs are activated to expand , differentiate to myofibroblasts , and produce increased levels of ecm and pro - inflammatory mediators , thereby also promoting macrophage accumulation ; all of these reactions constitute a coordinated response to tissue injury and the progressive accumulation of collagen . we found that animals dosed orally with 1  ml / kg of ccl4 in mineral oil on day 0 and day 7 and euthanized on day 10 exhibited significant liver fibrosis based on increased picrosirius red ( psr ) and collagen immunopositivity in liver tissue as compared to vehicle - treated mice . percent area of tissue immunoreactive with asma , the hallmark of myofibroblasts , and iba-1 , a macrophage - specific marker , were also increased in ccl4-treated animals , suggesting an increase in myofibroblasts and macrophage numbers in the tissue ( figures s1a  s1c ) . corresponding to increased collagen deposition in tissue , col1a1 mrna was markedly upregulated ( figure  s2a ) . the assessment of col1a1 mrna levels was used as an expedient approach to functionally assess a relatively large number of sirna targets and flag fibrogenic genes of interest . given the sirna - lnp delivery system targets multiple liver cell types ,  including hscs , hepatocytes , kupffer cells , but not endothelial cells,11 , 12 , 13 , 14 we validated the ability to kd genes expressed in hscs in mice using sirna - lnps ( figures s1d  s1h ) . mice were injected with a single dose of sirna - lnp ( 1  mg / kg ) against reelin ( reln - lnps ) , a gene prominently expressed in hscs . animals receiving reln - lnps , but not luc - lnps ( negative control ) , showed significant reduction in reln mrna in both oil- and ccl4-treated animals ( figure  s1e ) . we also demonstrate the ability to kd the hsc - specific gene col1a1 in liver . col1a1-lnp reduced col1a1 mrna by 50% , with no effect on reln mrna . because the tgf- pathway is a recognized fibrogenic mediator , we tested whether kd of transforming growth factor beta receptor 1 ( tgfbr1 ) decreased the transcription and accumulation of collagen indeed , administration of tgfbr1-lnps and col1a1-lnps , but not luc - lnps , reduced the transcription of col1a1 mrna ( figure  s1f ) as well as collagen accumulation ( figures s1 g and s1h ) . to test the role of genes identified in our transcriptomic analysis , we modified the sirna administration protocol to allow pre - existing protein to be turned over and therefore achieve pre - clearance of proteins encoded by the genes of interest and ensure continued kd throughout the experiment ( figure  2a ) . we first used luc - lnps to validate this protocol for detecting differences between oil- and ccl4-treated cohorts , with respect to the levels of col1a1 mrna , our primary parameter for defining target genes . we similarly evaluated differences in the mrna levels for a panel of other genes ( figure  s2a ) . the strictly standardized mean difference ( ssmd ) values for col1a1 , collagen type iii alpha 1 ( col3a1 ) , tissue inhibitor of metalloproteinases 1 ( timp1 ) , and platelet derived growth factor receptor beta ( pdgfrb ) , but not tgfbr1 , integrin subunit alpha m ( itgam ) , or alpha - actin-2 ( acta2 ) , suggested that these genes were suitable to use as readout genes ( figure  s2b ) . having established the suitability of the ccl4 treatment with the sirna - lnp kd platform , we proceeded to test our candidate genes . we tested the effect of kd of 24 individual genes in  vivo ( figure  2a ) . we identified seven genes that significantly reduced the amount of col1a1 mrna ( table 1 ) . five of these seven genes , namely , early growth response 2 ( egr2 ) , fk506-binding protein ( fkbp10 ) , atpase na / k transporting subunit alpha 2 ( atp1a2 ) , follistatin like 1 ( fstl1 ) , and hyaluronan synthase 2 ( has2 ) , were silenced by 75%100% in whole liver tissue and significantly reduced col1a1 mrna levels . silencing of these genes did not elicit any overt toxicity , as measured by body weight loss ( figure  s3 ) . because the other 19 genes did not meet these criteria ( tables 1 and s1 ) , we focused on further analyzing the five modifiers of fibrosis . the kd of the transcription factor ( tf ) egr2 had the broadest effect by reducing the levels of col1a1 and col3a1 mrna as well as the percent area immunopositive for asma , iba-1 , and collagen proteins , although the reduction in collagen by half did not achieve significance ( figure  2 ; table 1 ) . egr2 silencing also reduced pdgfrfb mrna levels ( figure  3a ) , suggesting a mechanism of decreased fibrosis and macrophage accumulation due to reduced expansion of activated hscs . given that egr2 has structural similarities to its family members , egr1 and egr3 , we confirmed the specificity of the egr2 sirna by transfecting human 293 t cells with plasmids encoding mouse egr1 , egr2 , or egr3 under the transcriptional control of the cmv promoter . in the tested construct , the expression of mouse egr1 , egr2 , and egr3 is not under the control of the endogenous promoter . thus , reduction in gene expression is due to the binding of the sirna egr2 to the mrna . using the same egr2 sirna that was used for the in  vivo experiments , we found that this sirna reduced egr2 mrna levels , but not egr1 or egr3 mrna , unambiguously demonstrating that egr1 and egr3 were not targeted by our egr2 sirna ( figures s4a  s4c ) . interestingly , even though the sirna against egr2 was specific , egr2-lnp treatment in  vivo also reduced egr1 and egr3 mrna levels ( figures s4d  s4f ) , suggesting transcriptional co - regulation of egr family members . similar to egr2 , kd of fkbp10 or atp1a2 reduced the collagen and iba-1 positive areas as did that of egr2 ( figure  2 ; table 1 ) ; however , kd of fkbp10 or atp1a2 did not alter the percent immunopositivity for asma . fkbp10 kd also decreased pdgfrb mrna levels , but this effect was somewhat weaker than that of egr2 ( 28% versus 45% reduction ; figures 3a and 3b ) , possibly accounting for its lack of effect on asma immunopositivity . these results suggest that egr2 , fkbp10 , and atp1a2 regulate fibrogenesis through different mechanisms . kd of either of the two other fibrosis modifiers , fstl1 and has2 , decreased col1a1 mrna levels as well as macrophage and collagen accumulation in the tissue , although has2 kd did not achieve significance for the latter ( figure  2 ) . notably , kd of these genes showed increased amounts of asma immunopositivity , which was highly significant for has2 kd , suggesting increased numbers of myofibroblasts . corresponding with the increased myofibroblasts , we detected higher levels of the timp1 mrna in the livers of both fstl1 and has2 kd animals ( figures 3d and 3e ) . despite the increased myofibroblast and timp1 mrna levels , , we also conducted in parallel in  vitro studies and showed that three independent fstl1 sirnas reduced col1a1 mrna levels . thus , the direct reduction of collagen transcription may be sufficient to reduce overall collagen accumulation , despite the elevated number of myofibroblasts and reduced collagen degradation ( figure  s5 ) . egr2 emerged as the strongest modifier of fibrosis in the ccl4 liver fibrosis / sirna kd platform . we aimed to directly demonstrate that egr2 kd in an hsc would reduce col1a1 and pdgfrb mrna levels . we tested the effect of egr2 kd in the human hsc line lx-2 using three independent sirnas ( figure  4 ) . this experiment required a different sirna to be used because no species cross - reactive sirna was designed ( see materials and methods ) . all sirnas against erg2 reduced its mrna levels as well as col1a1 and pdgfrb mrna levels . these results are consistent with the effect of egr2 kd in our in  vivo ccl4 assay platform . as expected , kd of col1a1 reduced its own expression , but not that of egr2 . interestingly , col1a1 silencing also reduced pdgfrb mrna levels , suggesting a regulatory feedback loop . taken together , these data support that egr2 is fibrogenic in both rodents and humans . myofibroblasts are key effectors of fibrosis and therefore the elucidation of pathways that promote this cell lineage would be a significant advance in the field . here , we report the use of a targeted in  vivo sirna kd platform to identify novel mediators of fibrogenesis . we leveraged novel biology - targeting genes upregulated during liver and kidney fibrosis in myofibroblasts and their progenitors and sirna - lnp technology to silence these genes in ccl4-induced liver injury and fibrosis in mice . we identified five genes , namely egr2 , atp1a2 , fkbp10 , fstl1 , and has2 , which modified fibrogenesis in this system because their kd reduced col1a1 mrna levels and collagen accumulation in the liver . prior to our study , egr2 , atp1a2 , and fkbp10 had not been functionally validated in fibrosis in  vivo , and our study is the first to validate fstl1 and has2 in liver fibrosis . we further directly demonstrated that egr2 kd in human hscs reduced the expression of fibrotic genes . we conclude that in  vivo sirna kd using sirna - lnp is a powerful tool to identify disease - relevant modifiers in  vivo . studying myofibroblasts is extremely difficult because they are in direct contact with multiple cell types , including other myofibroblasts , myeloid cells , mesothelial cells , endothelial cells , epithelial cells , and circulating fibrocytes in  situ . myofibroblasts respond to chemo - mechanical stimuli , and therefore disruption of their native environment during cell isolations alters their activation state . to circumvent these shortcomings this approach allowed us to study functional consequences of silencing of specific genes of interest in  vivo in their native environment . it is of note that our sirna - lnp delivery approach allows us to silence gene expression in hscs ; however , this lnp formulation also targets other cells in the liver , including kupffer cells and hepatocytes , but not endothelial cells.11 , 12 , 13 , 14 therefore , we can not rule out the contribution of gene silencing in kupffer cells and hepatocytes on the anti - fibrotic outcome . our strategy to identify novel players in tissue fibrosis took advantage of recent transcriptomic datasets for myofibroblasts and their precursors isolated from fibrotic liver or kidney in mouse models.6 , 7 using sirna - lnp delivery technology , we tested the effect of gene silencing in a physiologically relevant context in  vivo . we employed a ccl4 administration model that featured significant collagen accumulation and established several robust rna and histological measures of fibrogenic activity to identify fibrosis modifiers . using the criteria of significant reduction of col1a1 mrna levels of greatest interest are those that also reduced collagen accumulation in the tissue , namely , egr2 , atp1a2 , fkbp10 , fstl1 , and has2 . notably , these five genes were the most strongly silenced in  vivo ( by 75% or greater ) that significantly reduced col1a1 mrna levels . it is possible that we might have achieved a greater degree of in  vivo gene silencing for other target genes with a higher dose or a different dosing paradigm . thus , although we can not draw conclusions about the functional role of genes whose silencing was incomplete , our results clearly identify five genes as novel fibrogenic mediators in liver fibrosis in  vivo . these genes fell into different groups based on the effects of their silencing on a mini mrna array and histological measures . we propose grouping them as follows : egr2 , a gene with potential function in hsc expansion ; atp1a2 and fkbp10 , genes that regulate mechanosensing ; and fstl1 and has2 , genes involved in matrix - cell signaling . egr2 , a tf of the egr family , emerged as the strongest effector in the sirna screen . egr tfs are important for the induction of cellular programs , including cell proliferation , differentiation , and death , in response to stimuli such as cytokines , growth factors , and toxic substances.17 , 18 , 19 
 egr2 has not been previously reported to play a role in tissue fibrosis in  vivo . rather , egr2 is prominently known for its role in early myelination of the peripheral nervous system in humans and mice.21 , 22 , 23 egr2 kd silenced its expression by 94% , reduced col1a1 and col3a1 mrna levels , and diminished accumulation of myofibroblasts and macrophages , each to 15% of control . thus , our demonstration of reduced fibrogenic activity by egr2 silencing reveals a novel egr2 function . the downregulation of pdgfrb mrna levels by 45% suggests that egr2 may promote hsc expansion . this mechanism would explain the reduced asma - expressing hsc and thereby reduced signals for macrophage accumulation , all of which may have contributed to the reduction in fibrosis . importantly , we also directly demonstrate that egr2 kd in a human hsc line reduces col1a1 and pdgfrb transcripts , supporting the regulation of a fibrogenic program in hscs by this tf . the comparatively weaker effect of col1a1 silencing in  vitro ( 50% ) versus in  vivo sirna - lnp- mediated kd ( 80% ) of egr2 is due to the use of different delivery methods as well as the use of different sirna sequences to achieve gene silencing in human cells versus in  vivo in mice . interestingly , col1a1 kd also reduced pdgfrb mrna levels . because ecm accumulation and expression of pdgfr are linked , it is feasible that decreased collagen production would decrease pdgfr expression in a feedback loop . the relevance of our finding to fibrogenesis in humans is also supported by a prior in  vitro study , in which overexpression of egr2 in human skin fibroblasts upregulated collagen gene expression and myofibroblast differentiation , and these profibrotic tgf--induced responses were attenuated by egr2-depletion . we unambiguously confirmed that our sirna is specific for egr2 ( figures s4a  s4c ) and also found that egr2 kd in human hscs reduces col1a1 and egr2 , but not egr1 or egr3 mrna ( data not shown ) , confirming the anti - fibrotic effect . however , in  vivo egr2 kd also decreased the mrna levels of egr1 and egr3 , supporting the co - regulation of these transcription factors . it has been shown previously in  vitro that overexpression of egr1 induced egr2 and egr3 expression and likewise egr3 overexpression induced egr1 and egr2 expression , whereas overexpression of egr2 did not have such an effect.25 , 26 , 27 here , we report that egr2 kd also reduces egr1 and egr3 mrna levels in  vivo , a novel aspect of egr tf co - regulation . although egr2 has not been previously reported to play a role in tissue fibrosis in  vivo , 
 egr1 and egr3 have been implicated in the fibrotic response in  vitro and in  vivo previously.27 , 28 , 29 therefore , the extent of the anti - fibrotic effect may be due the co - regulation of egr tf . atp1a2 and fkbp10 were identified as fibrogenic modulators , with mechanisms of action distinct from that of egr2 . atp1a2 is the large catalytic subunit of an enzyme that catalyzes the hydrolysis of atp , coupled with the exchange of na and k ions across the plasma membrane . atp1a2 has a known function in the contractility of skeletal and cardiac muscle . on the basis of our results , it is intriguing to speculate that atp1a2 silencing may have reduced the contractility of myofibroblasts , effectively mimicking their behavior in a more compliant matrix environment and hence resulting in a less fibrogenic program.31 , 32 supporting this theory , we found that atp1a2 kd reduced fibrogenesis in  vivo by significantly limiting the transcription of the col1a1 and timp1 genes , macrophage accumulation , and collagen deposition . however , in contrast to egr2 , atp1a2 silencing had no effect on asma expression . this theme of altered mechanosensing may also explain the alleviated fibrosis upon silencing of fkbp10 . fkbp10 is a molecular chaperone in the endoplasmatic reticulum that directly interacts with collagen i , contributing to the maturation and molecular stability of type i procollagen , promoting normal secretion , stable inter - molecular crosslinking , and collagen deposition in the ecm . fkbp10 has a well - established role in bone mass in humans and mice , but its role in tissue fibrosis is largely unexplored . fkbp10 mouse embryos are post - natally lethal and display reduced collagen crosslinking in calvarial bones . in this study , we report our novel finding that fkbp10 silencing reduces collagen deposition to 64% ; however , we found no reduction in asma expression . we speculate that fkbp10 silencing alters the quantity or crosslinking of the ecm and thereby alters mechanosensing . consistent with our in  vivo results , a prior in  vitro study showed that inhibition of fkbp10 in primary human fibroblasts from pulmonary fibrosis patients attenuates expression of pro - fibrotic genes and decreases collagen secretion . thus , our findings expand the scope of atp1a2 and fkbp10 functions by demonstrating their novel roles in fibrogenesis in  vivo . fstl1 and has2 silencing also affected multiple fibrogenic readouts in our assay platform , but resulted in a pattern distinct from that of egr2 , atp1a2 , and fkbp10 . both fstl1 and has2 silencing reduced col1a1 mrna , macrophage accumulation , and collagen deposition but counterintuitively increased myofibroblast density as well as the level of timp1 mrna , a profibrotic pathway . fstl1 is a tgf--inducible , secreted , matricellular glycoprotein belonging to the secreted protein acidic rich in cysteines ( sparc ) family , with pleiotropic functions.37 , 38 , 39 during the course of our current study , it was reported that haplodepletion of fstl1 or a neutralizing anti - fstl1 antibody reduced bleomycin - induced pulmonary fibrosis . we found that fstl1 gene silencing reduced collagen transcription and accumulation in the context of liver fibrosis but surprisingly increased myofibroblast density , possibly reflecting a compensatory feedback loop . we hypothesize that the overall reduction in collagen accumulation is a direct result of decreased collagen production by fstl1 silencing , so that the increased number of myofibroblasts and reduced activity of collagen degradation would still result in decreased overall collagen accumulation . indeed , kd of fstl1 in the human hsc line lx-2 resulted in decreased levels of col1a1 mrna . like fstl1 , silencing of has2 resulted in unexpected effects , also likely through a matrix - modifying mechanism . has2 is one of three isoenzymes responsible for cellular hyaluronan ( ha ) synthesis , a matrix component normally highly abundant in joint synovium . ha is well known for its role in wound healing.40 , 41 
 has2 has been previously reported to be pathological in the bleomycin - induced lung fibrosis model , in which mesenchymal cell - specific deletion of has2 abrogated myofibroblast accumulation and inhibited the development of lung fibrosis . we show that has2 silencing also reduces liver fibrosis . has2 silencing also unexpectedly augmented myofibroblast density . taken together , our studies have identified fstl1 and has2 as fibrosis modifiers in our in  vivo liver fibrosis platform and reveal surprising effects of silencing these genes , supporting their complex role in the regulation of fibrogenesis . to the best of our knowledge , our study is the first example of a targeted sirna assay for the discovery of novel fibrosis modifiers in  vivo . the results implicate five novel modulators of the fibrogenic process in the liver and demonstrate the feasibility of sirna - lnp for interrogating gene function in the liver . because these genes are commonly upregulated in pericytes and myofibroblasts in both liver and kidney fibrosis models , these genes warrant follow - up to further explore their function and mechanism of action in liver and other organs and to  inform their role in both physiological and pathophysiological fibrosis . all animal procedures were conducted in accordance with cambridge laws and the institutional animal care and use committee approved protocol no . 657 . 8- to 10-week - old male balb / c mice ( taconic bioscience ) were given ccl4 by oral gavage ( sigma - aldrich ; 1  ml / kg dissolved in mineral oil ) on days 0 and 7 . sirna - lnps were either given on day 6 at a concentration of 1  mg / kg or on days 7 , 3 , 2 , and 6 at 0.5  mg / kg . these experiments were controlled with cohorts of mice that received sirna targeting luciferase ( luc - lnps ) with oil or ccl4 , respectively . animals were sacrificed on day 10 by co2 asphyxiation , followed by cardiac perfusion with 10  ml of pbs , and liver lobes were harvested for qpcr and histopathological analysis . target fibrosis modifier genes were selected by intersecting transcriptomic data generated from isolated hscs of mice treated with ccl4 for 2  months , and genes specifically expressed in pericytes / myofibroblasts in the mouse kidney 2 and 5  days after uuo.6 , 7 169 genes were common to both datasets and significantly upregulated at least 2-fold . among these , we excluded transcripts that were not associated with  a  gene product , genes with no known human homolog , and well - known modifiers in fibrosis ( validated in  vivo by independent investigators ) , resulting in 24 genes for further evaluation . a set of ten sirnas with canonical sirna structures was defined for each target gene . all sirnas were directed against the coding sequence of their respective target mrna and perfectly matched all known mrna transcript variants of the target gene available in the ncbi reference sequence database ( refseq db , release 65 , may 2014 ) . the sirnas used in this study were designed to be a perfect match only to their target mrna and to have  2 mismatches within positions 218 of the 19-mer sequence to any other genes . the design as described above is sufficient to ensure specificity.44 , 45 specifically , at least four mismatches between the sirna sequences of each of the genes of interest ( atp1a2 , egr2 , fkbp10 , fstl1 , and has2 ) and the readout genes col1a1 , col3a1 , tgfbr1 , pdgfrb , and timp1 . further , there were three mismatches between egr2 sirna and egr1 sequence and more than four mismatches between egr2 sirna and egr3 sequence . the sirna antisense strands lacked a seed region ( nucleotides 27 ) identical to a seed region ( nucleotides 27 ) of known mouse mirnas ( mirbase , release 20 , june 2013 ) . additional similarity analysis showed that there were at least four mismatches between the sirna sequences of each of the genes of interest ( atp1a2 , egr2 , fkbp10 , fstl1 , and has2 ) and the readout genes col1a1 , col3a1 , tgfbr1 , pdgfrb , and timp1 . further , there were three mismatches between egr2 sirna and egr1 sequence and more than four mismatches between egr2 sirna and egr3 sequence . from sirnas fulfilling those criteria , the sirnas for each final screening were selected for predicted activity based on analysis with proprietary algorithms ( axolabs , gmbh ) . in  vivo grade sirnas each of the ten aforementioned sirna molecules / target gene was screened for its potency to kd gene expression in a cell - based reporter assay . for each target gene , the regions of their coding sequences targeted by the sirnas were synthesized and inserted to the 3 utr of the firefly luciferase gene in pcmvluc(v ) validation vector ( origene ) , generating an mrna transcript that encodes both the luciferase ( luc ) and the target gene when expressed in cells . 293 t cells ( atcc ) were co - transfected with the validation vector and the sirnas of interest using lipofectamine 2000 ( thermo fisher scientific ) at 0.1- , 1- , and 10-nm concentrations to assess kd potency . sirna binding to the target gene leads to degradation of the mrna that encode both luc and the target gene , resulting in reduced luciferase expression . luciferase expression was measured 24  hr post - transfection in a functional assay using the promega luciferase assay system ( promega ) , and luciferase signal was detected using a perkinelmer envision plate reader ( perkinelmer ) . for each sirna tested , the efficiency of kd was calculated as percent of luciferase signal compared to cells that were transfected with the validation vector for the target gene alone . the most potent sirna for each target gene ( > 90% luciferase reduction at 1  nm sirna in  vitro ) was then formulated as lnps for in  vivo application . the sequences of the sirnas that resulted in efficient kd ( 75% ) of their target gene in  vivo are provided in table s2 . sirna - lnps were synthesized by mixing together an sirna - containing aqueous phase with a lipid - containing ethanol phase . the  aqueous phase contained sirna in 10  mm citrate buffer ( ph  3 ) . the ethanol phase contained ionizable lipid c12 - 200 ( wuxi apptec ) , 1,2-distearoyl - sn - glycero-3-phosphocholine ( dspc , avanti polar lipids ) , cholesterol ( sigma ) , and 1,2-dimyristoyl - sn - glycero-3-phosphoethanolamine - n-[methoxy(polyethylene glycol)-2000 ] ( c14 peg 2000 , avanti ) at a 50:10:38.5:1.5 molar ratio and 5:1 c12 - 200:sirna weight ratio . the aqueous and ethanol phases were mixed together at a 3:1 volume ratio in a microfluidic chip device using syringe pumps as previously described to a final sirna concentration of 1  mg / ml . the resultant sirna - lnps were dialyzed overnight in a 20,000 molecular weight cut - off ( mwco ) cassette against 1x pbs at 4c . on average , sirna - lnps had a mean diameter of approximately 100150  nm , with a polydispersity index between 0.1 and 0.2 as measured by dynamic light scattering ( zetasizer , malvern instruments ) . the sirna encapsulation efficiency ( approximately 60% ) was determined using a modified quant - it ribogreen assay ( invitrogen ) as previously described . this sirna delivery system can target multiple liver cell types , including hepatocytes , monocytes , tissue macrophages , dendritic cells ( dcs ) , and myofibroblasts / pericytes . rna was extracted from liver lobes with trizol ( thermo fisher scientific ) , followed by on - column dnase digestion using rneasy mini kits and rnase free dnase sets ( both from qiagen ) . cdna was generated using the high - capacity cdna reverse transcription kit , and qpcr analysis was performed using pre - designed taqman primer probes ( table s3 ) on the quantstudio 12k flex real - time pcr system ( all from thermo fisher scientific ) to determine the expression levels of sirna target genes and a panel of readout genes . the relative expression of each gene was normalized to mouse or human glyceraldehyde 3-phosphate dehydrogenase ( gapdh ) using the 2 method , and the mean expression of each gene in the control group was set to 100% . 5-m sections were stained with psr using an automated leica system according to standard protocols . anti - alpha - smooth muscle actin ( sma ) ihc was performed with fluorescein isothiocyanate ( fitc)-labeled mouse monoclonal antibody 1a4 , followed by an anti - fitc antibody and periodic acid - schiff ( pas ) counterstain ( abcam ) . anti - iba-1 and anti - col1a1 ihc was performed using rabbit polyclonal antibodies ( cat . # 019 - 19741 , wako , and genetex , cat . the iba-1 positive area was limited to infiltrating macrophages by differentiating isolated single cells ( largely kupffer cells lining sinusoids ) from larger round cells in clusters ( periportal , largely infiltrating macrophages , sometimes filled with mineral ) . de novo interstitial asma was differentiated from endogenous asma in muscular arteries by performing a pas counterstain and then programming the algorithm to distinguish vascular morphology for exclusion from quantification . the human hsc line lx-2 ( emd millipore ) was transfected with  sirnas against egr2 , col1a1 , or luc with lipofectamine rnaimax reagent following the manufacturer s protocol . we used the following predesigned sirnas ( listed as sirna i d ) : egr2 : s4540 , s4541 , and s4542 ; col1a1 : s3276 ; luc : am4629 ; and fstl1 : s22032 , s22033 , and s22034 ( thermo fisher scientific ) . rna was isolated on day 2 after transfection using rneasy kits with on - column dnase digestion ( qiagen ) . 293 t cells ( atcc ) were co - transfected with pcmv - egr1 , pcmv - egr2 , or pcmv - egr3 ( origene ) and either silencer select negative control no . 1 sirna ( thermo fisher scientific ) or sirna against egr2 ( axolabs ) . rna was isolated 24  hr after co - transfection using rneasy plus kits ( qiagen ) . statistical significance was assessed by two - tailed , unpaired student s t test for comparison of two groups or by one - way anova , followed by the dunnett s multiple comparisons test for more than two groups . we used an ssmd method to analyze the suitability of the mrna readout genes in our sirna - lnp / ccl4 liver fibrosis platform ( figure  s2 ) . the ssmd method was used to calculate the median of differences divided by the standard deviation of the differences between the group treated with oil and the control sirna targeting luciferase ( oil / luc - lnp ) and versus ccl4/luc - lnp . it represents the average fold change penalized by the variability of the fold change across animals . ssmd is a suitable method because it takes the variability between animals into account . for our study , a  threshold of   < 1.5 was taken as the cut - off value for suitable readout genes .
fibrotic diseases contribute to 45% of deaths in the industrialized world , and therefore a better understanding of the pathophysiological mechanisms underlying tissue fibrosis is sorely needed . 
 we aimed to identify novel modifiers of tissue fibrosis expressed by myofibroblasts and their progenitors in their disease microenvironment through rna silencing in  vivo . 
 we leveraged novel biology , targeting genes upregulated during liver and kidney fibrosis in this cell lineage , and employed small interfering rna ( sirna)-formulated lipid nanoparticles technology to silence these genes in carbon - tetrachloride - induced liver fibrosis in mice . 
 we identified five genes , egr2 , atp1a2 , fkbp10 , fstl1 , and has2 , which modified fibrogenesis based on their silencing , resulting in reduced col1a1 mrna levels and collagen accumulation in the liver . 
 these genes fell into  different groups based on the effects of their silencing on a transcriptional mini - array and histological outcomes . 
 silencing of egr2 had the broadest effects in  vivo and also reduced fibrogenic gene expression in a human fibroblast cell line . prior to our study , egr2 , atp1a2 , and 
 fkbp10 had not been functionally validated in fibrosis in  vivo . 
 thus , our results provide a major advance over the existing knowledge of fibrogenic pathways . 
 our study is the first example of a targeted sirna assay to identify novel fibrosis modifiers in  vivo .