Overexpression of Dimethylarginine Dimethylaminohydrolase 1 Attenuates Airway Inflammation in a Mouse Model of Asthma

Levels of asymmetric dimethylarginine (ADMA), an endogenous inhibitor of nitric oxide synthase, are increased in lung, sputum, exhaled breath condensate and plasma samples from asthma patients. ADMA is metabolized primarily by dimethylarginine dimethylaminohydrolase 1 (DDAH1) and DDAH2. We determined the effect of DDAH1 overexpression on development of allergic inflammation in a mouse model of asthma. The expression of DDAH1 and DDAH2 in mouse lungs was determined by RT-quantitative PCR (qPCR). ADMA levels in bronchoalveolar lavage fluid (BALF) and serum samples were determined by mass spectrometry. Wild type and DDAH1-transgenic mice were intratracheally challenged with PBS or house dust mite (HDM). Airway inflammation was assessed by bronchoalveolar lavage (BAL) total and differential cell counts. The levels of IgE and IgG1 in BALF and serum samples were determined by ELISA. Gene expression in lungs was determined by RNA-Seq and RT-qPCR. Our data showed that the expression of DDAH1 and DDAH2 was decreased in the lungs of mice following HDM exposure, which correlated with increased ADMA levels in BALF and serum. Transgenic overexpression of DDAH1 resulted in decreased BAL total cell and eosinophil numbers following HDM exposure. Total IgE levels in BALF and serum were decreased in HDM-exposed DDAH1-transgenic mice compared to HDM-exposed wild type mice. RNA-Seq results showed downregulation of genes in the inducible nitric oxide synthase (iNOS) signaling pathway in PBS-treated DDAH1-transgenic mice versus PBS-treated wild type mice and downregulation of genes in IL-13/FOXA2 signaling pathway in HDM-treated DDAH1-transgenic mice versus HDM-treated wild type mice. Our findings suggest that decreased expression of DDAH1 and DDAH2 in the lungs may contribute to allergic asthma and overexpression of DDAH1 attenuates allergen-induced airway inflammation through modulation of Th2 responses.


Introduction
Asthma is a chronic inflammatory lung disease characterized by airway hyperresponsiveness (AHR), airway inflammation, excess mucus production and pulmonary remodeling [1]. Arginine metabolism has been found to play a critical role in the pathogenesis of allergic asthma [1,2] and mounting evidence suggests that nitric oxide (NO) bioavailability plays an important role in the development of allergic inflammation [3,4]. Arginine is metabolized by arginase yielding urea and L-ornithine that is further metabolized to polyamines and prolines, which can modulate cell proliferation and collagen production [4]. Arginine is also metabolized by nitric oxide synthases (NOS) including neuronal NOS (nNOS), inducible NOS (iNOS) and endothelial NOS (eNOS) [5]. eNOS is expressed in bronchial epithelium and type II alveolar epithelium, nNOS is expressed in airway nervous tissue, and iNOS is expressed in type II alveolar epithelium, lung fibroblasts, airway and vascular smooth muscle cells. The constitutive NO produced by nNOS and eNOS is important for smooth muscle relaxation, bronchodilation and determination of vascular tone and blood pressure while the inducible NO produced by iNOS has pro-inflammatory effects [5].
Asymmetric dimethylarginine (ADMA) is a by-product released from proteolysis of methylated proteins and ADMA competitively inhibits all three NOS by displacing L-arginine from NOS [6]. ADMA also competes with arginine for cellular uptake by cationic amino-acid transporters, affecting the cellular ADMA/arginine ratio [7,8]. Analysis of methylarginine metabolism in the cardiovascular system showed that the lung is a major source of ADMA [9]. ADMA has been shown to have profound effects on multiple tissues. Microarray studies showed that pathophysiolog-ical concentrations of ADMA elicit significant changes in the gene expression in coronary artery endothelial cells [10]. Treatment of primary mouse lung fibroblasts with ADMA induced arginase activity and collagen production [11]. ADMA infusion resulted in increased lung resistance and decreased compliance in response to methacholine in mice, which was associated with significantly increased pulmonary collagen deposition [11]. ADMA potentiates ovalbumin-induced airway inflammation in a mouse model of asthma [12]. In humans, plasma ADMA levels are increased in severe asthma patients compared to nonsevere asthma patients and control subjects [3]. Other studies showed that ADMA levels are increased in sputum and exhaled breath condensate from asthma patients [13][14][15]. A more recent study showed that lower L-arginine/ADMA ratios are associated with reduced lung function and increased respiratory symptom frequency in subjects with late-onset asthma [16].
DDAH activity is a key determinant of intracellular ADMA concentration [17]. Ninety percent of ADMA is metabolized by DDAH and the rest is excreted through the kidneys. DDAH metabolizes ADMA to generate citrulline and dimethylamine. There are two isoforms of DDAH in human, mouse, rat and other species, DDAH1 and DDAH2 [18][19][20]. Immunostaining of human lung tissues showed expression of DDAH1 and DDAH2 in both alveolar and bronchiolar epithelium [21]. One study showed that homozygous DDAH1-deficient mice die before birth while the heterozygous DDAH1-deficient mice have a 20% increased level of ADMA and develop severe endothelial dysfunction [22]. Anther study showed that DDAH1-deficient mice are viable and have significantly increased ADMA levels [23]. ADMA levels are decreased in DDAH1-transgenic mice by 50% compared to wild type mice. Further, DDAH1-transgenic mice display a significant increase in NOS activity and a decreased risk of endothelial dysfunction compared to the wild type mice [24][25][26]. DDAH1 is the major enzyme responsible for metabolizing ADMA whereas DDAH2 has no detectable role in degrading ADMA in vivo [23]. In this paper, we determined the role of DDAH1 in the development of AHR and airway inflammation.

Mice
Mice were maintained and handled under Institutional Animal Care and Use Committee-approved procedures (Cincinnati Children's Hospital Medical Center, Protocol Number: 2D10082) and the Guide for the Care and Use of Laboratory Animals (Institute of Laboratory Animal Resources, National Research Council). C57BL/6 wild type and DDAH1-transgenic mice were purchased from Jackson Laboratory (Bar Harbor, ME). Immunization C57BL/6 wild type and DDAH1-transgenic mice were immunized with intratracheal instillation of 100 mg of house dust mite (HDM) (Dermatophagoides pteronyssinus) extract (Greer Laboratories, Lenoir, NC) in 50 ml of PBS or 50 ml of PBS alone 3 times per week for 3 weeks as previously described [27,28].

Analysis of AHR and airway inflammation
Twenty four hours after the last treatment, AHR in response to methacholine (0, 50, 100 and 200 mg/ml) was measured by flexiVent (SCIREQ, Montreal, Canada) as previously described  [29]. After AHR measurement, blood, bronchoalveolar lavage (BAL) samples and lung tissues were harvested as previously described [30]. The BAL total and differential cell counting as well as histological staining of lung sections by hematoxylin and eosin (H&E) or periodic acid Schiff (PAS) (Thermo Fisher Scientific, Waltham, MA) were performed as described by the manufacturer or described previously [30]. ELISAs for total IgE and IgG1 or HDM-specific IgE and IgG1 in BAL fluid (BALF) and serum were done as previously described [31].

RNA-Seq and Bioinformatics Analysis
Total RNA was isolated from mouse lungs using Trizol reagent (Invitrogen, Carlsbad, CA), digested with RNase-free DNase and purified using an RNeasy MinElute kit (QIAGEN, Valencia, CA). Equal amounts of RNA were pooled from each mouse lung in an experimental group (n = 4 per group) and analyzed in duplicate. RNA-Seq was performed by the Genomics Sequencing Core in the University of Cincinnati. The RNA-Seq library was constructed using a PrepX SPIA RNA-Seq kit (IntegenX, Pleasanton, CA) and Apollo 324 NGS Library Prep System (IntegenX). 10 ng of total RNA was converted into cDNA suitable for mRNA sequencing. The cDNA was then sheared by Covaris S2 (Covaris, Woburn, MA) under the conditions recommended by IntegenX, followed by Bioanalyzer assay of the size distribution with Agilent High Sensitivity DNA kit (Agilent, Santa Clara, CA). The properly sheared cDNA fragments were purified by Agencourt AMPure XP magnetic beads (Beckman Coulter, Brea CA). Using the IntegenX PrepX ILM DNA library kit for Illumina and Apollo 324 NGS Library Prep System, 500 ng of purified cDNA fragments were then put through end repair, addition of a single 'A' base and ligation of adapters, and indexed individually. The products were purified and enriched by PCR to create the final cDNA library targeting mRNAs. The size of the generated library was validated by Bioanalyzer and the library was quantified using the Kapa Library Quantification kit (Kapa Biosystems, Woburn, MA). Six individually indexed cDNA libraries were equal amount pooled for clustering in cBot system (Illumina, San Diego, CA). Libraries at the concentration of 6.5 pM were clustered onto a flow cell using Illumina's TruSeq SR Cluster Kit v3, and sequenced for 50 cycles using the TruSeq SBS kit on Illumina HiSeq system.
Sequence reads were aligned to the reference genome (mm10) using TopHat aligner [32]. The counts of reads aligning to each gene's coding region were summarized using ShortRead and associated Bioconductor packages for manipulating and analysis of next-generation sequencing data and custom-written R programs [33,34]. Statistical analysis to identify differentially expressed genes for each comparison was performed using the negative-binomial model of read counts as implemented in the DESeq Biocondoctor package [35]. P-values were adjusted for multiple comparisons based on false discovery rates (FDR) [36]. Differential expressions with adjusted p-values of ,0.05 were considered statistically significant.
Functional enrichment analysis was performed using the logistic regression based LRpath methodology [37]. The gene lists used in the functional enrichment analysis were from genes associated with Gene Ontology terms and KEGG pathways. The statistical significance of gene list enrichment was determined by the False Discovery Rate cut-off of 0.1. Genes that were both members of at least one statistically significant gene list and had differential expression p-values of ,0.01 were considered to be differentially expressed for the purpose of network analysis. Ingenuity Pathways Analysis (IPA) (Ingenuity Systems, Mountain View, CA) was used to identify integrated and interconnected biological networks and upstream targets for the differentially expressed genes between two groups (with adjusted p-values of ,0.05). RNA-Seq data have been deposited with the NCBI Gene Expression Omnibus (http:// www.ncbi.nlm.nih.gov/geo) under accession number GSE49047.

RT-PCR
Total RNA was isolated from lung tissues using Trizol reagent (Invitrogen), digested with RNase-free DNase and purified using an RNeasy MinElute kit (QIAGEN). Reverse transcription was done using Oligo-dT First-Strand cDNA Synthesis Kit (GE Healthcare, Piscataway, NJ). Quantitative PCR (qPCR) was done using the SYBR Green Master Kit and LightCyclerH 480 instrument (Roche Diagnostics, Indianapolis, IN). All primers used are listed in Table S1.

Mass Spectrometry (MS)
MS was performed by the Analytical & Mass Spectrometry of Small Molecules Core at the University of Cincinnati. The procedure to quantify ADMA was adapted from Schwedhelm et al

Statistical analysis
All values are expressed as mean 6 SD. The data were analyzed with a 2-tailed unpaired student's t-test with Welch's correction or 1-way ANOVA with Newman-Keuls' post test using Prism 5.0c for Mac OS X from GraphPad Software (San Diego, CA). A pvalue of ,0.05 was considered statistically significant.

Results
Expression of DDAH1 and DDAH2 is decreased in the lungs and ADMA levels are increased in BALF and serum from HDM-treated mice Since ADMA levels are increased in allergic asthma and DDAH is responsible for the majority of ADMA metabolism, we hypothesized that the increased levels of ADMA are due to decreased expression of DDAH following allergen exposure. We determined the expression of DDAH1 and DDAH2 in the lungs in a mouse model of asthma. C57BL/6 mice were challenged by intratracheal instillation of 100 mg HDM, 3 times per week for 3 weeks. This is a more clinically relevant protocol of chronic allergen exposure with an increased dosage of HDM for the C57BL/6 mouse strain. As shown in Fig. 1A and 1B, HDM exposure resulted in decreased expression of DDAH1 and DDAH2 in mouse lungs, which correlated with increased levels of ADMA in BALF and serum ( Fig. 1C and 1D).

Overexpression of DDAH1 attenuates HDM-induced airway inflammation
We next determined the effect of DDAH overexpression on airway inflammation and AHR in a HDM-induced asthma model using DDAH1-transgenic mice since DDAH1 is the major enzyme responsible for ADMA metabolism [23]. The mice were challenged by intratracheal instillation of 100 mg HDM 3 times per week for 3 weeks. Twenty-four hours after the last treatment AHR was measured by flexiVent and airway inflammation was assessed by BAL cell counts. BAL total cell and eosinophil numbers were decreased in HDM-treated DDAH1-transgenic mice compared to HDM-treated wild type mice ( Fig. 2A-2B). No significant differences in inflammatory cell infiltration in airways determined by H&E staining or in mucus production determined by PAS staining of lung sections were observed between HDMtreated DDAH1-transgenic mice and HDM-treated wild type mice (data not shown). No significant difference in AHR was observed in HDM-treated DDAH1-transgenic mice compared to HDM-treated wild type mice (data not shown).
IgE levels in BALF and serum are decreased in DDAH1transgenic mice compared to wild type mice after HDM treatment We next tested whether allergen sensitization was affected by the overexpression of DDAH1. As shown in Fig. 3A and 3B, the total IgE levels were decreased in BALF and serum in HDMtreated DDAH1-transgenic mice compared to HDM-treated wild    type mice but there was no difference in the levels of total IgG1 in BALF or serum ( Fig. 3C and 3D). The level of HDM-specific IgE was decreased in serum from HDM-treated DDAH1-transgenic mice compared to HDM-treated wild type mice but the level of HDM-specific IgG1 in serum was unaffected ( Fig. 3E and 3F).

Gene expression profiles are significantly different in the lungs of wild type and DDAH1-transgenic mice
To elucidate the potential mechanisms underlying attenuated HDM-induced airway inflammation in DDAH1-transgenic mice, we determined gene expression profiles in the lungs from C57BL/ 6 wild type and DDAH1-transgenic mice treated with PBS or HDM (intratracheal instillation of 100 mg of HDM in 50 ml of PBS or 50 ml of PBS alone 3 times per week for 3 weeks) by RNA-Seq.
There were 707 genes differentially expressed in the lungs of HDM-treated wild type and HDM-treated DDAH1-transgenic mice (Table S3). The top downregulated genes (with an adjusted p-value of ,0.05) are shown in Table 3 (grouped and sorted by fold change), including metabolism/transport genes (Fbp1, Dhrs9, Arg1, etc.), immune/defense response genes (Ear11, Chia, Sod3, etc.), cytokine/chemokine genes (Ccl22, Tnfsf12, Ccl9 and Il33), extracellular matrix genes (Col6a2, Muc5ac, Muc4 and Col1a1) and microRNAs (Mir5109 and Mir5107). The top upregulated genes (with an adjusted p-value of ,0.05) are shown in Table 4 (grouped and sorted by fold change), including ribosomal complex genes (Rps21, Rps12, Rpl30, etc.) and ribosomal pseudogenes (Gm11968, Gm12191, Gm13253, etc.). The top network for the downregulated immune/defense genes and cytokine/chemokine genes was associated with inflammatory response, cellular movement, and hematological system development and function (Fig. 5A). The top network for the upregulated ribosomal complex genes and ribosomal pseudogenes is associated with cellular development, cancer, and cellular assembly and organization (Fig. 5B). The top upstream targets for both downregulated and upregulated genes were IL-13 and FOXA2 (Fig. 5C).
The expression of IL-13, IL-4, IL-33, CCL11, ARG1, MMP-12, CHIA, CHI3L3 and CHI3L4 is decreased in the lungs from DDAH1-transgenic mice compared to wild type mice following allergen challenge Based on RNA-Seq results, we verified the expression of some key cytokines, chemokines, matrix metallopeptidases and chitinases including IL-13, IL-4, IL-33, CCL11, MMP-12 and CHIA.  As shown in Fig. 6A and 6B, induction of IL-13 and IL-4 expression in the lungs was attenuated in HDM-treated DDAH1transgenic mice compared to wild type mice. Expression of IL-33 was decreased in lungs from DDAH1-transgenic mice after HDM treatment compared to wild type mice (Fig. 6C). We found that the expression of CCL11 and MMP-12 was decreased in HDMtreated DDAH1-transgenic mice compared to HDM-treated wild type mice (Fig. 6D, 6E). The expression of ARG1, a key enzyme in arginine/NO metabolism, was significantly decreased in HDMtreated DDAH1-transgenic mice compared to HDM-treated wild type mice (Fig. 6F). The difference in expression of iNOS (NOS2) or FOXA2 was not significant between wild type and DDAH1transgenic mice (Fig. 6G, 6H). The expression of acidic chitinase (CHIA) and chitinase like proteins (CHI3L3 and CHI3L4) was decreased in the lungs of HDM-treated DDAH1-transgenic mice compared to HDM-treated wild type mice ( Fig. 7A-7C).

Discussion
Our data reveal that expression of DDAH1 and DDAH2 is decreased in the lungs in a mouse model of asthma, and overexpression of DDAH1 attenuates allergen-induced airway inflammation. Asthma is a condition of decreased NO bioavailability. ADMA is an endogenous inhibitor of NOS, which is a major source of NO. DDAH1 and DDAH2 are responsible for metabolism of over 90% of ADMA in vivo, and our data support a role for DDAH downregulation in asthma pathogenesis.
The mechanism of the observed downregulation of DDAH may be epigenetic modulation. Studies using mouse trophoblast stem cells and trophoblastic tissues of postimplantation mouse embryos showed DNA methylation-dependent epigenetic regulation of DDAH2 gene expression. The CpG island in the DDAH2 promoter was hypermethylated in trophoblast stem cells but hypomethylated in differentiated cells [39]. We found that the mouse DDAH1 promoter also contains a CpG island [40], but their methylation status and relationship to gene expression are unknown. Previous studies have shown that allergen exposure results in altered methylation status of IL-4 and IFNc promoter CpG islands [41]. It is likely that allergen exposure results in hypermethylation of mouse DDAH1 and DDAH2 promoter CpG islands, which results in decreased expression of DDAH1 and DDAH2 in mouse lungs.
Our data showed that overexpression of DDAH1 attenuated allergen-induced airway inflammation although it had no significant effect on AHR. It is not surprising as studies showed uncoupled airway inflammation and AHR in allergen challenged C57BL/6 mice due to strain-dependent genomic factors [42]. The infiltration of eosinophils into the lungs was decreased in HDMexposed DDAH1-transgenic mice, which is consistent with the result that expression of CCL11 was decreased in lungs from HDM-treated DDAH1-transgenic mice. The total IgE and HDMspecific IgE levels in BALF or serum were decreased in DDAH1transgenic mice after HDM exposure, suggesting that overexpression of DDAH1 affected Ig class switch.
The RNA-Seq data showed that overexpression of DDAH1 results in decreased lung expression of multiple immune/defense response genes that are associated with a network of inflammatory responses. The top upstream target is iNOS. A previous study showed that increased ADMA levels increase the expression of iNOS in mouse lungs [12], suggesting the effect of DDAH1 overexpression on iNOS expression/activity may be mediated by altered ADMA levels. Following HDM exposure, the expression of ARG1 is significantly decreased in the lungs from DDAH1 transgenic mice, suggesting arginine/NO pathways play important roles in the effect of DDAH1 overexpression on allergic airway inflammation. We found the expression of genes involved in mucus production (Clca3, Muc5ac and Muc4) and collagen synthesis (Col6a2 and Col1a1) is also decreased. Another previous  study showed that increased levels of ADMA resulted in increased pulmonary collagen deposition [11], suggesting DDAH1 may regulate collagen synthesis through modulation of ADMA levels. Network analysis suggests that DDAH1 regulates mucus production gene expression through IL-13/FOXA2. RT-qPCR data further showed that the expression of IL-13, IL-4 and CCL11 is decreased in HDM-treated DDAH1-transgenic mice, which is consistent with attenuated eosinophil infiltration in airways and decreased serum and BALF IgE levels. Although the difference in expression of iNOS and FOXA2 is not significant between the wild type and DDAH1-transgenic mice, overexpression of DDAH1 may directly or indirectly affect the activity of iNOS and FOXA2. Interestingly, the expression of acidic chitinase (CHIA) and chitinase like proteins (CHI3L3 and CHI3L4) was decreased in the lungs of HDM-treated DDAH1-transgenic mice. Acidic chitinase has been shown to paly important roles in asthma [43]. CHI3L3 and CHI3L4 are rodent specific chitinase like proteins that are induced by Th2 cytokines or allergen challenge and have chemotactic activity [44][45][46]. Whether acidic chitinase and chitinase like proteins can modulate the expression of DDAH1 requires further studies. Overexpression of DDAH1 may have non-specific effects on gene expression. Generation of cell type specific transgenic mice with different expression levels of the DDAH1-transgene would help minimize the non-specific effects.
In summary, our data suggest that decreased expression of DDAH1 and DDAH2 in lungs may contribute to allergic asthma and overexpression of DDAH1 attenuates allergen-induced airway inflammation through modulation of Th2 responses.