Genome-wide effect of pulmonary airway epithelial cell – specific Bmal1 deletion to

: Pulmonary airway epithelial cells (AECs) form a critical interface between host and environment. We investigated the role of the circadian clock using mice bearing targeted deletion of the circadian gene brain and muscle ARNT-like 1 ( Bmal1 ) in AECs. Pulmonary neutrophil infiltration, biomechanical function, and responses to influenza infection were all disrupted. A circadian time-series RNA sequencing study of laser-captured AECs revealedwidespreaddisruptioningenesofthecorecircadianclockandoutputpathwaysregulatingcellmetabolism(lipidsandxenobiotics),extracellularmatrix,andchemokinesignaling,butstrikinglyalsothegainofanovelrhythmictranscriptomein Bmal1 -targetedcells.ManyoftherhythmiccomponentswerereplicatedinprimaryAECs cultured in air-liquid interface, indicating significant cell autonomy for control of pulmonary circadian physiology. Finally, we found that metabolic cues dictate phasing of

For this reason, BMAL1 knockout models have limited utility in defining the impact on specific cell types and organs. To address this, many studies have used Cre-LoxP technology to achieve tissue-specific disruption of Bmal1, commonly leading to disrupted homeostatic physiology specific to that tissue (4).
There is mounting evidence for an important role of the circadian clock in pulmonary physiology and disease. Repeated phase shifting of environmental lighting schedules (i.e., jet lag or shift work protocols), which is known to disrupt circadian activity and organ-wide rhythmicity, leads to augmented lung tumor formation in p53 knockout mice (5), with a similar phenotype observed in a background of the Per2 knockout strain. We have previously shown that pulmonary responses to injury such as bleomycin-mediated oxidative damage are augmented in the circadian clock mutant ClockD19 mice, driven by dysregulation of the E-box-regulated nuclear factor erythroid-derived 2-like 2 and glutathione-mediated antioxidant defense pathway (6). Global loss of Bmal1 also leads to hyperinflammatory pulmonary responses to respiratory viral infections (7)(8)(9). Selective gene targeting of Bmal1 within bronchial airway epithelial cells (AECs) causes increased inflammatory responses to smoking (10), whereas our earlier studies have shown impaired pulmonary glucocorticoid signaling, leading to a selective augmentation of specific chemokines and a hyperinflammatory response to LPS (11). Thus, multiple pulmonary pathologies are linked to the circadian clock.
Here, we report the impact of circadian gene targeting using Cre-lox to delete Bmal1 in club cell secretory protein (CCSP)-expressing pulmonary AECs and reveal widespread perturbation of rhythmic transcriptional control within these cells. In vivo, we observed increased background pulmonary neutrophilia and impaired lung biomechanical function. We also observed impaired recovery to pulmonary influenza viral infection with activation of multiple cytokine and chemokine genes.
We next characterized the circadian transcriptome of AECs using laser-capture microdissection (LCM) from Bmal1 flox/flox and CCSP-Bmal1 2/2 mice over 2 circadian cycles. Rhythmic genes in AECs were enriched for metabolism, cell matrix interaction, and chemokine signaling, the majority of which were disrupted upon Bmal1 deletion. Strikingly, we also detected an emergent novel rhythmic transcriptome in Bmal1-targeted AECs, which lacked enrichment for canonical clock gene-regulating elements, implying that the intact clockwork in these cells may block the imposition of exogenous circadian signals to AECs. Using time-feeding schedules, we show that the pulmonary clockwork and responses to endotoxin challenge are phase reversed when mice are fed in the light phase, indicating a critical role for systemic metabolic cues in entrainment of the pulmonary clock. Finally, we cultured AECs using an air-liquid interface (ALI) culture model, and this revealed a remarkably similar transcriptional mechanism to those of intact tissue, suggesting an important autonomous role for the AEC in pulmonary physiology. Thus, the autonomous circadian clockwork in AECs is indispensable for normal pulmonary homeostasis and innate immunity and is entrained by systemic rhythmic feeding-associated metabolic cues.

Animals
All experiments on animals were carried out in accordance with the UK Home Office Animals (Scientific Procedures) Act of 1986 and European Directive 2010/63/EU. Mice were group housed under the following controlled conditions: 12-h light/dark cycle at 21°C with free access to standard rodent chow and water unless otherwise stated. CCSP-Bmal1 2/2 and control Bmal1 flox/flox mice were bred in a Per2-luc background. Global Bmal1 2/2 mice were bred by intercrossing of Bmal1 +/2 mice, a gift from Prof. Akhilesh Reddy (University of Cambridge, Cambridge, United Kingdom). C57BL/6 mice were commercially sourced from Envigo (Huntingdon, United Kingdom). All experiments used male mice aged between 10 and 20 wk unless otherwise stated. For collection of tissues in the light/dark cycle, mice were sampled at zeitgeber time (ZT)0-24, which by convention indicates lights on at ZT0. For circadian collections, mice were maintained in constant darkness and samples collected 1 cycle after transfer to darkness at circadian time (CT), which by convention anchors expected time of lights off and activity onset to CT12.

Lung function test
CCSP-Bmal1 2/2 and control Bmal1 flox/flox mice (males, n = 3-15, 4 or 12 mo old) were used for lung function tests in flexiVent Fx1 (Scireq, Montreal, Canada). Anesthetic reagents were hypnorm (0.315 mg/ml fentanyl; 10 mg/ml fluanisone) and midazolam (5 mg/ml). They were mixed with water for injection, with a ratio of 1:2:1 (hypnorm, water, and midazolam, respectively), and 0.1 ml/10 g mouse weight was used for intraperitoneal injection. After tracheal cannulation, mice underwent mechanical ventilation with default settings of the machine. The lungs were expanded with deep ventilation, and the pressure volume graph was examined for air flow leakage or voluntary respiration. Then, lung function was measured by the default mechanical scanning protocol and the average value of 3 repeated readings was calculated for comparisons.

Influenza infection
CCSP-Bmal1 2/2 and control Bmal1 flox/flox mice (males, n = 6-10, 10-20-wk-old) were intranasally infected with 5 plaqueforming units (pfu) of influenza A virus, Puerto Rico/8/34, and H1N1 in 50 ml PBS in the morning. Weight was checked every day to monitor disease development. In the first cohort, mice were culled at d 21 postinfection. In a second cohort, mice were culled at d 11 postinfection without BAL sample collection. Left lungs were fixed with 4% paraformaldehyde for paraffin embedding. Right lungs were frozen for RNA. Hematoxylin and eosin staining was performed in lung sections cut at 5-mm thickness, with images obtained from CaseView software (3DHistech, Budapest, Hungary) after scanning slides in Pannoramic 250 Flash III (3DHistech). Lung injury histology scoring was undertaken as previously described in Bayes et al. (12), with peribronchial infiltrate from 0 to 4 and alveolar involvement from 0 to 3. Six samples from each group were randomly selected and three 310 images were scored per mouse.  (13). After culling (pentobarbital, intraperitoneally), lungs were inflated with ;0.8 ml 50% optimal cutting temperature/PBS v/v solution via a tracheal cannula. Dissected lung samples were frozen on dry ice and stored at 280°C. Samples were cut into 10-mm sections in a cryostat machine onto polyethylene naphthalate membrane slides (11600288; Leica Microsystems, Buffalo Grove, IL, USA). After serial dehydration in ethanol solutions (100, 75, 50, 95, 100, and 100%, 30 s each), slides were air-dried and then dissected in LMD6500 (Leica Microsystems). Distal bronchiolar epithelium was collected ;200 nm from the alveolar opening. RNA was purified by the PicoPure RNA Isolation Kit (Thermo Fisher Scientific) with on-column DNase digestion. RNA integrity was checked by the Agilent RNA 6000 Pico Kit (Agilent Technologies, Santa Clara, CA, USA). Poly (A)-selected mRNA (100 ng total RNA/sample, RNA integrity number $6) was used for RNA sequencing (RNA-seq). The library was constructed with TruSeq SBS Kit v.3 (Illumina, San Diego, CA, USA) with low-input method. Samples were sequenced in HiSeq 2500 (Illumina) with 4 samples per lane, generating 100 3 100 base-pair end reads. After filtering and mapping, differentially expressed (DE) genes across all 12 time points were called using DE-Seq2 (14) with a false discovery rate of 0.05. For rhythmic gene detection, MetaCycle (https://cran.r-project.org/web/packages/MetaCycle/index.html) (15) was used with meta2d (P , 0.01). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed by the Database for Annotation, Visualization and Integrated Discovery (DAVID) v.6.8 (https://david.ncifcrf. gov/) (16,17), and gene ontology term analysis was carried out in Enrichr (https://cran.r-project.org/web/packages/enrichR/index. html) (18,19). Phase set enrichment analysis (20) was used to analyze rhythmic pathways at different times of day [q , 0.05 with rhythmic genes (P , 0.05 by MetaCycle)]. RNA-seq data were stored in Array Express (E-MTAB-6384; European Bioinformatics Institute, Cambridge, United Kingdom; https:// www.ebi.ac.uk/arrayexpress/). Normalized read counts, rhythmic gene detection, and DE genes are shown in the Supplemental Data.
ALI culture of primary tracheal epithelial cells ALI culture was performed as described in You et al. (21). After overnight digestion of tracheas in pronase solution [0.15% w/v in membrane-thermotolerant Escherichia coli (mTEC) basic medium (DMEM and Ham's F-12) supplemented with 15 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid, 4 mM L-glutamine, 3.6 mM NaHCO 3 , 100 U/ml penicillin, and 100 mg/ml streptomycin], cells were seeded in collagencoated cell inserts for 24-well plates. Cells were cultured for 7 d submerged in mTEC Plus buffer [mTEC basic medium supplemented with 10 mg/ml insulin, 5 mg/ml transferrin, 0.1 mg/ml cholera toxin, 25 ng/ml epidermal growth factor, 30 mg/ml bovine pituitary extract, and 5% fetal bovine serum (v/v)] to expend cell number. Then, cells were differentiated in air-liquid culture condition (no medium on top of the insert) with mTEC serum-free buffer (mTEC basic medium supplemented with 5 mg/ml insulin, 5 mg/ml transferrin, 0.025 mg/ml cholera toxin, 5 ng/ml epidermal growth factor, 30 mg/ml bovine pituitary extract, and 1 mg/ml bovine serum albumin) for 9 d to induce cilia genesis. Experiments were carried out in mixed male and female mice. At the end of the experiment, cells were fixed for immunofluorescence staining or frozen for RNA extraction.

Adenovirus Cre-mediated Bmal1 knockdown in primary tracheal cells and RNA-seq
The transfection was undertaken as described in Balasooriya et al. (23). High titter Ad5CMVCre virus (2 3 10 e11 pfu/ml) was purchased from the Viral Vector Core Facility (University of Iowa, Iowa City, IA). After overnight culture of freshly isolated tracheal epithelial cells from Bmal1 f/f mice, 0.5 ml virus solution was added into each well (;1000 pfu) with serum-free mTEC Plus solution. After 8 h, the virus medium was changed back to normal mTEC Plus medium and the cells were continually cultured as described above. For the Bmal1 knockdown efficiency study, cells differentiated in ALI condition were put into photomultiplier tubes to monitor Per2-luc activity after temperature synchronization [2 cycles of 12 h at 36.5°C and 12 h at 38.5°C (24)]. Two wells per unit of a 24-well plate were made in house to fit samples into the photomultiplier tubes. For time-serial RNA-seq, samples cultured at 37°C were collected after 2 d of temperature synchronization at 1 well per group per time point for 48 h. Poly(A) selected mRNA from 1 mg total RNA were used for RNA-seq in an Illumina HiSeq 4000 machine, with 12 samples per lane generating 75 3 75 base-pair end reads. The sequence data were analyzed as previously described except that rhythmic gene detection was carried out analyzing only the first cycle of RNA-seq data instead of 2 cycles of RNA-seq data due to weak gene oscillation in vitro. RNA-seq data were stored in Array Express (E-MTAB-6384).
Food entrainment and LPS challenge C57BL/6 mice (males, 8-9-wk-old) were bought from Envigo. They were housed in normal condition initially for 1-2 wk before being kept in light-controlled cabinets. Mice were entrained to the light/dark cycle for 2 wk and then fed either on their day or night time for 2 wk before tissue collection or LPS treatment. In the cabinet, half (n = 12 in the first cohort and n = 16 in the second cohort) were maintained in the normal light/dark cycle (lights on at 7 AM and lights off at 7 PM), and the other half in reverse light/dark cycle (lights on at 7 PM and lights off at 7 AM). After 2 wk in the light/dark condition described above, restricted feeding was introduced. Mice were kept in the same light/dark condition with free access to water, and standard chow food was provided from 10 AM to 4 PM every day for 2 wk. In this way, half of the mice were fed in their subjective night and half were fed in their subjective day. Mice in cohort 1 were culled at the end of restricted feeding for lung tissue harvest. For mice in cohort 2, the aerosolized LPS exposure experiment was carried out at ZT0 (7 AM) and ZT12 (7 PM) (n = 8, feeding group/time point) as previously described in Gibbs et al. (11). BAL was collected with lung samples 5 h later. Cell count and fluorescence-activated cellsorting analysis of BAL cells were performed as previously described. CXCL5 levels were measured with a CXCL5 ELISA Kit (DY254; Bio-Techne, Minneapolis, MN, USA), according to product protocol.

Statistics
Values are means 6 SD. Data were analyzed using Prism v.5.0 (GraphPad Software, La Jolla, CA, USA). Unpaired Student's t tests (with Welch's corrections when variance was unequal) and 1-and 2-way ANOVAs were used, as appropriate. Significance was set at P , 0.05.

RESULTS
Background BAL neutrophilia, altered lung functions, and impaired influenza defense in CCSP-Bmal1 2/2 mice Using flow cytometry, we analyzed BAL cells and identified significantly enhanced background neutrophilia in CCSP-Bmal1 2/2 mice (Fig. 1A, B). Next, we assessed lung mechanical function by flexiVent in 4-mo-old mice, revealing elevated lung compliance (the measure of the ability of the lung to stretch and expand), reduced elastance [the reciprocal of compliance (i.e., the pressure change that is required to elicit a unit volume change)] and a trend for reduced resistance of the respiratory tract to airflow. By 1 yr of age, these phenotypes were reversed in CCSP-Bmal1 2/2 mice (Fig. 1C-E). To test the potential involvement of changes to the extracellular matrix, we assessed collagen levels in the lung using picrosirius red staining. This revealed an increased collagen deposition in proximity to large airways of 1-yr-old mice (Fig. 1F). This indication of pulmonary scarring is the long-term consequence of increased neutrophil infiltration in CCSP-Bmal1 2/2 mouse lungs.
Our earlier studies have shown that CCSP-Bmal1 2/2 mice exhibit exaggerated acute inflammatory responses to LPS, mediated by abnormal regulation of selective chemokine signaling (11). To assess the generality of such responses, we adopted a chronic inflammatory model using an influenza challenge, which generates a strong local response in pulmonary tissues and damages airway epithelium, leading to a pattern of sustained inflammation with a well-defined time course over several days and weeks (25). Postinfection, we observed significantly greater weight loss in CCSP-Bmal1 2/2 mice over a 21-d course (Fig. 1G) but without an effect on mortality. This was confirmed in repeat experiments in which mice were culled on a single day (d 11 postinfection), and here histologic analysis revealed substantial regions of pulmonary consolidation in CCSP-Bmal1 2/2 mice compared with controls ( Fig. 1H and Supplemental Fig. S1C). This was characterized by extensive T-cell marker CD3 immunofluorescence staining (Supplemental Fig. S1A), and, using a NanoString inflammatory gene array, we defined augmented expression of a number of chemokines and cytokines, including the lymphocyte attractant chemokine Ccl20 (Supplemental Fig. S1B). Thus, Bmal1 deletion in AECs imposes detrimental effects on both pulmonary homoeostasis and viral defense.  intervals over 2 circadian cycles (48 h). RNA-quality RNA and RNA-seq read-distribution was similar among samples (Supplemental Fig. S2G). Additionally, we compared gene transcript levels from RNA-seq at a single time point (CT8, 2 replicates per group) with those obtained by qPCR, and this revealed a near 1:1 concordance in expression (Supplemental Fig. S2E).
First, we examined core canonical circadian clockregulating gene expression (Fig. 2). Strong circadian rhythms were observed in these genes from samples of Bmal1 flox/flox mice. Among them, E-box-and D-boxregulated genes (Tef, Chrono, Dbp, Nr1d1/Reverba, Nr1d2/ Reverbß, Per1, Per3, Dec1, and Dec2) were all suppressed in the CCSP-Bmal1 2/2 group, compatible with the established direct action of Bmal1 on these circadian elements. RAR-related orphan receptor element (RRE)-regulated transcripts (Npas2, Cry1, Rorc, Nifl3/E4bp4, and Clock) were all up-regulated in CCSP-Bmal1 2/2 group, compatible with the suppression of Bmal1-regulated Reverb expression, which normally represses RRE elements on these genes. A third group of genes showed little change, including Hlf, Cry2, Rora, and Per2, indicative of complex control via both E-box and RRE-regulating motifs and additional systemic cues. As an additional validation for potential cross-contamination with nontargeted cells, we examined RNA-seq read coverage over the floxed region (exon 8) of the Bmal1 gene and showed that this region was deleted in all samples obtained from AECs of CCSP-Bmal1 2/2 mice (Supplemental Fig. S2I).
We used the MetaCycle meta2d function to identify circadian rhythmic transcripts. A total of 916 transcripts were detected as rhythmic in Bmal1 flox/flox animals, which were nonrhythmic in CCSP-Bmal1 2/2 mice (Fig. 3A top  panels, group 1). An additional set of 126 transcripts retained rhythmicity in both genotypes (Fig. 3A, middle panels, group 2). Unexpectedly, we also identified 856 transcripts that gained rhythmicity in CCSP-Bmal1 2/2 mice (Fig. 3A, bottom panels, group 3). This latter cohort of genes also exhibited a markedly different peak phase of expression compared with groups 1 and 2, displaced by ;9 h (Fig. 3B). Examples of genes in groups 1, 2, and 3 were shown in Supplemental Fig. S4A, with KEGG pathway Figure 2. Disrupted core clock gene expression in CCSP-Bmal1 2/2 mice. Time-series (48 h) RNA-seq studies were carried out in LCM distal small AECs and RNA reads were normalized by DEseq2 method. Transcripts of core clock genes from Bmal1 flox/flox (black solid circle) and CCSP-Bmal1 2/2 (red solid square) mice are shown. Y-axis labels DESeq2-normalized reads count, and x-axis labels CT. enrichment (Fig. 3C). Relating to increased background neutrophilia phenotype, in addition to Cxcl5, both Cxcl3 and Cxcl15 were strongly rhythmic in showed no enrichment for these circadian regulatory elements (Fig. 3D). Rhythmic pathways were identified by different time of day in the lungs (Supplemental Fig. S3). The rhythmic pathways are largely different between genotypes (Supplemental Fig. S3A). In the common  (left). B) Phase distribution of rhythmic genes of groups 1-3, which centers at CT5 and CT15 in groups 1 and 2, and at CT9 and CT20 in group 3. C ) KEGG pathway enrichment analysis of genes in groups 1-3, ranked by significance, with most significant pathways as propionate metabolism in group 1, circadian rhythm in group 2, and antigen processing and presentation in group 3. D) Examples of enriched DNA motifs of rhythmic genes shown above in order of significance of fit. rhythmic pathway, some show changed phases (Supplemental Fig. S3B). For example, the phase of metabolism of lipids and lipid proteins was shifted from CT18 in Bmal flox/flox mice to ;CT10 in CCSP-Bmal1 2/2 mice. Genotype-specific rhythmic pathways are shown in Supplemental Figs. S3C, D, most of which are in the daytime. In Bmal1 flox/flox mice, pathways related to extracellular matrix are centered at early day, followed by pathways like pallet-activation signaling and aggregation, chemokine-signaling pathway in midday, and metabolism in the evening. In CCSP-Bmal1 2/2 mice, some of the rhythmic pathways are related to immune response, like antigen presentation and processing in early day.
Additionally, using DESeq2, we detected 814 transcripts with significantly increased (n = 507) or decreased (n = 307) expression across all time points in CCSP-Bmal1 2/2 mice, defined as DE. Of these, 642 were nonrhythmic in either genotype (Fig. 4A, B). KEGG pathway analysis revealed enrichment in metabolismrelated pathways and inflammatory pathways (Fig. 4C). Within the metabolic pathways, 20 of 82 genes were rhythmic in Bmal1 flox/flox samples, with a range of peak phases in expression (CT12-24, Fig. 4D). Increased Hmgcs2 expression was validated by Western blotting, showing strong induction of a gene in the lungs of both CCSP-Bmal1 2/2 mice and in mice bearing a global deletion of Bmal1 (Bmal1 2/2 ; Fig. 4E and Supplemental Fig.  S4B), the expression of which is normally restricted to the liver and gut (26), indicating an impact on differentiated cell function.
Thus, Bmal1 disruption leads to not only widespread impact on metabolic and inflammatory pathways in AECs but also the striking emergence of an additional novel rhythmic transcriptome in targeted cells.
Cell-autonomous effect of Bmal1 in ALI differentiated primary tracheal cells in vitro Given the importance of cellular communication in vivo, it is crucial to address the question of whether the effect of Bmal1 in AEC clockwork is cell-autonomous or is regulated by local pulmonary signals. For this purpose, we adopted an in vitro ALI model to culture primary pulmonary tracheal epithelial cells derived from Bmal1 2/2 knockout mice and wild-type (WT) littermate mice. The ALI model is widely regarded as a gold standard for the study of primary airway epithelia. In both genotypes, we observed clear development of characteristic-ciliated cells (Fig. 5A-D), and chemokines previously shown as Figure 5. Cell-autonomous Bmal1 function in ALI culture of primary mouse tracheal cells. ALI cultures were established using primary tracheal epithelial cells from WT and global Bmal1 2/2 mice, with 7 d under submerged condition and another 10 d exposed to air in the apical side of cells. A, B) Electron microscopy images for cilia in ALI primary tracheal cells from WT and global Bmal1 2/2 mice, respectively. C, D) Immunofluorescence staining of cilia markers (acetylated tubulin) and epithelial cell markers (E-cadherin) in ALI primary tracheal cells from WT and global Bmal1 2/2 mice, respectively. E, F) Cxcl5 and abnormally regulated in vivo in CCSP-Bmal1 2/2 mice were also significantly augmented by 20-40 folds in long-term cultured AECs (Fig. 5E, F). Thus, dysregulation of chemokine pathways operates in a cell-intrinsic manner upon Bmal1 deletion. In order to exclude potential contributions arising from tissue-wide reprogramming effects of cells derived from Bmal1-null mice (27), we next cultured AECs from Bmal1 flox/flox mice on a PER2::LUC reporter background and treated these cells with an adenovirusexpressing Cre driven by cytomegalovirus promoter, leading to loss of exon 8 in targeted cells (Fig. 5G). We confirmed the efficiency of targeting by tracking PER2:: LUC bioluminescence 2 wk after infection, revealing initial dampening and subsequent loss of overt circadian rhythmicity (Fig. 5H).
To define the impact on circadian-regulated pathways, AECs were seeded, Cre-transfected, and cultured in ALI for 7 d; then, cells were subjected to circadian temperature synchronization for 2 d, and subsequently RNA was collected every 4 h for 48 h (Fig. 5I). RNA-seq revealed 762 rhythmic transcripts in Bmal1 flox/flox cells, with 232 transcripts in Cre-transfected cells, only 6 of which were common to both data sets (Fig. 5J). Pathway enrichment of rhythmic genes in control cells revealed a strikingly consistent pattern with that observed in LCMderived Bmal1 2/2 AECs from whole lung (Supplemental Fig. S5 compared with Fig. 3C). Interestingly, despite the close similarity of cellular processes, the specific transcriptional composition of these 2 data sets showed relatively little overlap (Fig. 5J). Additionally, a total of 3041 genes showed significantly altered expression (Supplemental Fig. S6 and Supplemental Data) in which pathway analysis revealed common signatures in these data sets (Fig. 5K vs. Fig. 4C). These common elements included circadian clock-and chemokine-regulating genes (Supplemental Fig. S6).

Feeding is the zeitgeber of lung inflammatory response
The above studies indicate an important role for AEC circadian clock in regulating physiology, the emergence of a novel rhythmic transcriptome in Bmal1-deleted AECs in vivo, and the cell-autonomous nature of Bmal1 effect. But this leaves the unanswered issue of how the pulmonary clockwork may be entrained. Recent studies have highlighted the importance of systemic metabolic cues in entraining peripheral oscillators and, in particular, the role of timed feeding (28). We therefore examined the role of timed feeding cues on the circadian-gated pulmonary inflammatory response to LPS (Fig. 6A). Restriction of feeding to the light phase reversed the phase of expression of circadian clock genes like Nr1d1 and Dbp in the lung (Fig. 6B, C). We also assessed Bmal1 target genes such as Snx31, which also showed significant phase reversal in expression (Fig. 6D). We then tested the circadian gating mechanism driving pulmonary inflammatory responses to LPS (11). Night-fed mice treated with aerosolized LPS showed significantly higher inflammatory responses at ZT0 than ZT12, exhibited by higher numbers of BALtotal leukocyte, neutrophil infiltration, and CXCL5 levels. This is compatible with our earlier studies, with peak responses during the circadian daytime in mice (11). In contrast, mice confined to feeding during the day exhibited a clear phase reversal of inflammatory parameters (Fig. 6E-G). Thus, feeding cues act as a potent local synchronizing factor (ZT) in entraining Figure 6. Reverse feeding resets time of day variation in pulmonary LPS response. A) Schematic description of experiment design. The food-reversal experiment was undertaken using 2 separate cohorts. B-D) Food reversal without LPS treatment was performed in cohort 1. E-G) Aerosolized LPS exposure experiment at ZT0 vs. ZT12 was performed in cohort 2 and total BAL cells, neutrophils, and CXCL5 concentrations measured. Gene expression was determined in cohort 2. Data analyzed by 2-way ANOVA with post hoc test to examine time of day difference within genotypes (n = 6-8). *P , 0.05, **P , 0.01, ***P , 0.001, ****P , 0.0001 (significant time of day difference within genotypes).
circadian-regulated pathways in the lung and circadian gating of responses to inflammatory stimuli.

DISCUSSION
Here we show that Bmal1 plays an important role in a diverse range of cellular functions in AECs, regulating pulmonary neutrophil infiltration, biomechanical function, and influenza infection. Further, we show that the AEC clockwork retains a markedly autonomous phenotype in culture. Key findings are summarized in illustration Fig. 7. Our earlier studies have shown an important role for Bmal1 in AECs in the circadian gating of LPS-induced pulmonary inflammation (11), whereas other groups using the global Bmal1 knockout have described broadly similar pulmonary phenotypes, including an elevated pulmonary neutrophilia (7), aberrant responses to chronic inflammatory stimuli (8,9), and impaired mechanical function (29). We now show the pulmonary AEC is the critical cell type in mediating these effects.
Time-series analysis of laser-captured AECs revealed widespread disruption in the complex of canonical core clock genes in Bmal1-targeted cells, similar in directions to findings in liver and adipose tissue (30,31). We also identified a subset of rhythmic genes common to both WT and Bmal1-targeted AECs (group 2), with a similar phase angle for peak expression. These could be generated by systemic cues, 1 example of which is a glucocorticoidinduced gene Gilz (Supplemental Fig. S3). Interestingly, a substantial number of emergent rhythmic transcripts were detected in Bmal1-deleted mice (group 3), with a marked difference in phase angle. KEGG analysis of this latter data set revealed the emergence of pathways predominantly controlling immune responses, but this was not associated with migration of leukocytes to the pulmonary epithelium in Bmal1-deficient mice. Here, motif analysis revealed no significant enrichment of canonical circadian regulatory motifs. Remarkably, although disruption of Bmal1 in AECs impacted a significant proportion of the expressed transcriptome, overall there was relatively little change in the total proportion of rhythmic transcripts due to the emergence of novel rhythmic transcripts in targeted cells. This counterintuitive outcome is also consistent with  Figure 7. Illustration of key findings. Feeding sets the time of lung in vivo, driving rhythmic expression of core clock proteins like BMAL1 and more widespread rhythmic genes in club cells, involving mainly metabolic-and chemokine-signaling pathways. Deletion of Bmal1 in club cells abolishes these rhythmically expressed genes and uncovers a different set of newly rhythmic genes. This may be due to other rhythmic signals in vivo. CCSP-Bmal1 2/2 mice show altered lung mechanic functions and more severe influenza infection. In vitro, primary AECs were able to generate rhythmic gene expression under temperature entrainment, and deletion of Bmal1 disrupts genome-wide rhythmic gene expression without generating newly rhythmic gene expression. Flu stands for influenza in the picture. studies of the mouse liver in which Bmal1 is disrupted or during aging (32,33) and suggests that a hitherto unrecognized feature of the endogenous cellular clockwork may be the ability to mask exogenous rhythmic signals from adjacent cells and tissues.
The mechanisms involved in synchronization of peripheral tissue and rhythmic physiologic responses still remain very poorly understood (4). However, a number of studies have highlighted the importance of metabolicassociated cues and that restriction of feeding to abnormal circadian phases in mice is sufficient to phase-shift transcriptional rhythms of peripheral tissue (28). With a similar protocol, we showed clear phase reversal for expression of core clock-regulating genes, as well as matrix-associated and chemokine genes. This protocol reset the circadian gating mechanism driving endotoxin responses, with a 12 h phase-shift in feeding time-reversed mice. This implicates metabolic cues as important general regulators of rhythmic immunity. The clear translational implication is that timing of feeding should be considered in clinical studies of circadian pulmonary function (34).
A critical question is whether the transcriptional changes seen in targeted AECs arise as a consequence of cell-intrinsic Bmal1 function or because of the exposure of resident AECs to multiple extrinsic signals from neighboring cells and systemic cues. We tested this by using a cell culture model. This revealed strongly rhythmic expression for chemokine genes (e.g., Cxcl5), with loss of clock control and up-regulation in Bmal1-disrupted cells, phenocopying the situation in intact tissues (11). Further RNA-seq profiling and subsequent pathway analyses revealed a consistent and comparable enrichment of common pathways with those identified in laser-captured cells from whole lung and previous studies of whole-lung tissues (35). Processes involved in metabolism, xenobiotic detoxification, chemokine signaling, and extracellular matrix dominated these. Strikingly, despite the similarity of function, there was relatively little overlap of specific gene sets in our cell culture studies with those from lasercaptured cells, possibly reflecting differences in gene expression between small and large airway structures (36). Our data therefore show that many of the critical functions of Bmal1 in pulmonary epithelial cells are cell-intrinsic and not indirectly regulated by extrinsic factors.
Our approaches here employed LCM as a means to isolate pulmonary small airway cells. Such dissection does not allow fine-scale resolution to cellular subtypes, nor can we attribute high-resolution cell-specific phenotypes in ex vivo culture models. For this, single-cell imaging approaches could be used to allow dynamical assessment of individual gene products, currently an unrealized gold standard in the field. However, use of LCM does allow the instantaneous assessment of the state of specific cells from primary tissues, without concerns associated with cell sorting and subsequent phenotypical drift. Nonetheless, an impressive feature of our data is the persistence of primary phenotypes in cultured AECs.
In summary, we reveal that the circadian circuitry within the pulmonary AEC is essential for pulmonary physiology and immunity, the main functions of which were under circadian control (Fig. 7). An AEC-specific circadian transcriptome was revealed, as well as the effect of Bmal1 deletion, both in vivo and in vitro. Bmal1 regulation operates in a cell-intrinsic manner, and circadian clock in the lung is synchronized by metabolic feeding-associated cues. Collectively, this work has implications for the role of the clock in multiple pulmonary diseases such as cancer, asthma, and chronic obstructive pulmonary disease, in which the circadian clockwork is commonly disrupted.