|Home | About | Journals | Submit | Contact Us | Français|
miR-200b plays a role in epithelial-to-mesenchymal transition (EMT) in cancer. We recently reported abnormal expression of miR-200b in the context of human pulmonary hypoplasia in congenital diaphragmatic hernia (CDH). Smaller lung size, a lower number of airway generations, and a thicker mesenchyme characterize pulmonary hypoplasia in CDH. The aim of this study was to define the role of miR-200b during lung development. Here we show that miR-200b−/− mice have abnormal lung function due to dysfunctional surfactant, increased fibroblast-like cells and thicker mesenchyme in between the alveolar walls. We profiled the lung transcriptome in miR-200b−/− mice, and, using Gene Ontology analysis, we determined that the most affected biological processes include cell cycle, apoptosis and protein transport. Our results demonstrate that miR-200b regulates distal airway development through maintaining an epithelial cell phenotype. The lung abnormalities observed in miR-200b−/− mice recapitulate lung hypoplasia in CDH.
Every year, over 50,000 children are born with congenital diaphragmatic hernia (CDH) associated with abnormal lung development resulting in lung hypoplasia and persistent pulmonary hypertension1–3. CDH occurs as frequently as cystic fibrosis, but the pathogenesis is poorly understood. Smaller lung size, lower number of airway generations and a thicker mesenchyme characterize the abnormal lungs in CDH1. We have previously shown an inherent lung development defect in CDH4.
Lung development is a continuous process comprised of five developmental stages: embryonic, pseudo-glandular, canalicular, terminal saccular and alveolar5. Functionally, the early two stages are characterized by lung branching morphogenesis. Cell specification, vascularization and reduction of mesenchyme to form thin air-blood interfaces for gas exchange characterize the later stages6, 7.
MicroRNAs (miRNA) are small, non-coding RNAs that regulate gene expression through mRNA stability and translation8–10. They are essential for development and homeostasis of organs11–14. More than 1800 microRNAs have been identified in human15. Research focusing on the role of microRNAs in lung development and disease is limited. We recently discovered that miR-200b is elevated in abnormal lungs of human CDH babies. In the same study, we found that higher miR-200b expression in the fetal tracheal fluid of CDH fetus is associated with a better response to fetoscopic endoluminal tracheal occlusion (FETO, a prenatal therapy to promote lung growth)16.
MiR-200b belongs to the miR-200 family (miR-141, miR-429, miR-200a, miR-200b and miR-200c) and regulates epithelial-to-mesenchymal transition (EMT) in cancer and organ fibrosis17–20. Others have shown that miR-200 is down-regulated in a mouse model of fibrotic lung disease and human patients with idiopathic pulmonary fibrosis (IPF)21. Also, miR-200b can inhibit migration and invasion of non-small cell lung cancer cells22. The role of miR-200b during normal lung development has yet to be defined.
The goal of this study was to delineate the role of miR-200b during lung development using loss of function models in vivo. We generated a miR-200b−/−(KO) mouse to evaluate the functional impact of miR-200b absence on development in vivo. miR-200b deficient mice had stiffer lungs due to disturbed distal airway branching, thicker alveolar walls and downregulation of epithelial cell differentiation. Our data suggest that miR-200b is required to achieve the necessary balance in development of lung epithelial cells and fibroblasts to ensure development of a structurally and functionally effective respiratory organ.
We generated miR-200b−/− (KO) mice by replacing the complete miR-200b gene with a LacZ-reporter by targeted homologous recombination in C57Bl/6N mouse embryonic stem cells using the NorCOMM cassette (Fig. 1a,b, Supplementary Fig. 1). MiR-200b+/− mice were inter-crossed with the C57Bl/6N mice for at least 8 generations before use in experimental cohort studies. We confirmed complete knockout of miR-200b expression by RT-qPCR in fetal lungs (Fig. 1f) and lungs from 8-week old mice (Supplementary Fig. 5). We showed that mature microRNAs transcribed in the same cluster - miR-200a and miR-429 - were still expressed, albeit lower compared to miR-200b+/+ lungs (Fig. 1f, Supplementary Fig. 5). The expression of miR-200c and miR-141 did not change in the miR-200b−/− lungs compared to miR-200b+/+ lungs, suggesting that there were no compensatory effects on other family members (Fig. 1f and Supplementary Fig. 5).
Genotyping data from breeding miR-200b+/− x miR-200b+/− mice revealed an expected Mendelian distribution, indicating that these mice experience no embryonic lethality, and they were viable, fertile and appeared morphologically normal. In contrast to miR-200b/miR-429−/− mice reported by others23, our miR-200b−/− mice carried a lacZ-reporter gene that we exploited to localize miR-200b promoter activity during development (Fig. 1c). We observed a high expression in the endoderm of lung buds and developing airways during lung development. At E12.5, during branching morphogenesis, we observed lacZ staining mainly in the endodermal cells of the lungs and a faint expression in the mesenchymal cells. Notably, we observed high LacZ-miR-200b expression in other organs where development also hinges on epithelial-mesenchymal interactions and branching morphogenesis. These include the developing inner ear24, palate25 and mammary buds26 (Supplementary Fig. 3), consistent with prior reports suggesting that miR-200b expression is highly coordinated and associated with control of the development of these organs27–29.
Based on this high and dynamic expression of miR-200b, we focused first on the role of miR-200b during the early stages of lung development. We cultured mutant lungs from E11.5 mice as described before4. Branching patterns between the different lung explants were not affected in a significant manner (Fig. 1d,e). To evaluate if the absence of miR-200b resulted in compensatory upregulation of other family members, we assessed the abundance of all family members in fetal lung explants using RT-qPCR (Fig. 1f). Mature microRNAs transcribed in the same cluster–miR-200a and miR-429–were still expressed, albeit lower compared to miR-200b+/+ lungs (wt). The expression of miR-200c and miR-141 did not change in the miR-200b−/− fetal lungs compared to miR-200b+/+ lungs, suggesting that there were no compensatory effects on other family members. Newborn miR-200b+/− and miR-200b−/− mice did not display any breathing difficulties after birth, and together, these results indicate that proximal airway branching is not influenced by absence of miR-200b during development.
To determine whether lung function is affected in adult miR-200b−/− mice, we performed in vivo lung mechanics analyses using a flexiVENT small animal ventilator in 8-week-old miR-200b−/− mice and compared these with studies from miR-200b+/− and miR-200b+/+ mice. Mice were subjected to increasing doses of nebulized methacholine (MCh) to assess concentration dependent response characteristics of respiratory mechanics. Inhaled methacholine causes constriction of airway smooth muscle cells (bronchoconstriction)30. At rest, miR-200b−/− mice did not exhibit altered peripheral tissue/airway resistance (tissue damping), however after challenge with methacholine (MCh; 6mg/ml and higher) a substantive increase in tissue damping was revealed (Fig. 2a). Consistent with development of a low compliance or “fibrotic” lung, miR-200b−/− mice also showed significantly higher tissue elastance upon MCh-challenge (12mg/ml and higher) (Fig. 2b). At baseline, airflow resistance in conducting airways (Newtonian resistance) was not different between different groups of mice, but was significantly higher in miR-200b−/− mice challenged with high concentrations of MCh (Fig. 2c). Finally, miR-200b−/− mice demonstrated significantly lower hysteresivity (elastic hysteresis) in pressure-volume loops obtained before and after MCh challenge, consistent with the increased elastance we observed in these animals (Fig. 2d, Supplementary Fig. 4). To examine if the absence of miR-200b influenced the expression of other family members and therefore contributed to the observed phenotype, we measured the abundance of all miR-200b family members in 8-week-old wt and miR-200b−/− lungs. Like what we observed in fetal lungs, adult miR-200b−/− lungs had lower miR-200ba and miR-429 abundance, but no changes in miR-200c and miR-141 abundance (Supplementary Fig. 5). We also observed lower miR-200a and miR-429 in kidney tissues of miR-200b−/− mice (Supplementary Fig. 6).
Some of the lung function abnormalities in miR-200−/− could result directly from changes in surfactant function or distal airway branching. Therefore, we assessed the biophysical function of surfactant from miR-200b−/− and miR-200+/+ lungs using capillary surfactometry. We found that the surface tension-reducing capacity of surfactant in miR-200b−/− mice was markedly reduced (Fig. 2e), a finding that correlated with compromised labeling of Surfactant protein-B (SP-B) (Fig. 2f,g) and pro-Surfactant Protein-C (SP-C) in these lung tissues. Our data suggest a functional role for miR-200b in development of surfactant properties, resulting in reduced lung compliance in miR-200b−/− mice.
We evaluated lung morphometry in miR-200b−/− mice in a quantitative fashion, using two different, complementary approaches. We first performed in vivo high-resolution micro-CT scanning on live animals to eliminate the effects of tissue processing and inflation on the morphometry of the lungs. We made reconstructions of the lungs so that data between mice could be directly compared (i.e. if two images show the same level of grey scale, then they exhibit the same degree of x-ray attenuation). We then set a threshold to segment the data into solid tissue and air (using the same threshold value for all images) and calculated the volume of air in the lungs. To ensure that we used comparable areas of the lungs, we identified the carina in each lung and then moved 3.25mm above this. We then calculated the air volume from this point to the base of the lungs. Using in vivo high-resolution micro-CT scanning, we found that the density of lung parenchyma in miR-200b−/− mice is greater (gray area), and by measuring the airspace volume we demonstrated that distal alveoli are less air-filled compared to lungs from wildtype mice (Fig. 3a,b). This can result from increased alveolar space collapse or smaller alveolar air volume.
Hematoxylin and eosin (H&E) staining of miR-200b−/− lungs that were inflated and embedded in paraffin confirmed that these lungs have fewer septae and thicker alveolar walls when compared to miR-200+/+ lungs (Fig. 3c). We then scanned three complete sections per lung of three miR-200b+/+ and three miR-200−/− lungs using an Axio Scan.Z1 and used the ZEN Image Analysis software module to calculate the airspace area percentage of each lung. We observed lower percentages of airspace in miR-200−/− lungs (Fig. 3d).
Immunofluorescence studies for vimentin showed that miR-200b−/− lungs have more vimentin-positive cells, suggesting an increased presence of fibroblasts-like, mesenchymal cells (Fig. 3c,d). Also, miR-200b−/−lungs showed higher expression of Twist 1 protein – a transcription factor and marker for EMT - compared to miR-200b+/+ lungs (Fig. 3e,f). There were no differences in immunofluorescence patterns for other markers of lung fibroblast cell differentiation (Fibroblast growth factor-10) and epithelial cell differentiation (E-cadherin, cytokeratin, and CC-10) (data not shown). Taken together, these results indicate that even though miR-200b−/− mice do not experience obvious breathing difficulties, their lungs display a functional phenotype similar to lung fibrosis and lung hypoplasia observed in children with CDH.
To evaluate the effect of miR-200b on maintaining human bronchial epithelial cell phenotype, we transfected BEAS-2B cells with miR-200b inhibitors and performed double-immunofluorescence with cytokeratin and vimentin to mark epithelial cells and fibroblasts, respectively. We found that miR-200b inhibitors promoted accumulation of fibroblast cell markers following down regulation of miR-200b (Fig. 3g). To study if this also resulted in functional mesenchymal properties, we performed a scratch wound-healing assay (Fig. 3h). We transfected BEAS-2B cells with miR-200b inhibitors and 18h later, we observed 25% less scratch wound closure in the control group, whereas, in cultures transfected with miR-200b inhibitors, wound closure was enhanced by 75%.
In order to evaluate the effect on downstream targets and associated pathways owing to loss of miR-200b on the lung tissue transcriptome, we performed Next Generation Sequencing (NGS) on total RNA samples from lungs of three 8-week-old miR-200b−/− mice and three miR-200+/+ (wt) mice. Heat Map and unsupervised hierarchical clustering by sample and transcripts was performed on all samples passing QC using the top 500 genes that have the largest coefficient of variation based on FPKM counts (Fig. 4a). Table 1 shows the ten most differentially expressed mRNAs. The full list of differentially expressed transcripts is shown in Supplementary Table 3. We used Gene ontology (GO - Gene Ontology Consortium, 2000) enrichment analysis to identify GO terms that are significantly associated with differentially expressed protein coding genes. Using PANTHER Gene Ontology classification system31, we identified Notch and Wnt signalling among the most affected biological pathways and cytoskeletal as well as immunity proteins among the most affected protein class in miR-200b−/− lungs compare to miR-200b+/+(Fig. 4b,c). We used r package topGO32 in order to generate Go network for Biogolical process and evaluate Biological function. Our current findings revealed that the three most affected biological processes in miR-200b−/− mice lungs were related to cell cycle, apoptosis and protein transport (Fig. 4d). The most affected Biological Function is shown in Table 2. Among differentially expressed mRNAs, we confirmed the expression of Plunc, Cyp2a5 and Cdh26 (palate lung and nasal epithelial clone) by q-PCR (Fig. 4e).
We demonstrate for the first time that miR-200b plays a role in peripheral lung development by maintaining an epithelial cell phenotype. Our miR-200b deficient mice have lung function abnormalities, surfactant biophysical dysfunction with compromised pro-Surfactant Protein-C and surfactant protein-B expression. MiR-200b−/− lungs have decreased distal airway branching, a denser lung parenchyma with thicker alveolar walls, a higher number of fibroblast-like cells and over-expression of a marker for EMT: Twist.
We generated miR-200b−/− mice by targeted deletion of miR-200b in ES cells. We evaluated mature forms of all miR-200 family members to determine if there are any compensatory effects between the family members. Although miR-200a and miR-429 were expressed lower in miR-200b−/− lungs, their expression was not undetectable like miR-200b. Based on our studies, we cannot exclude that downregulation of miR-200a and miR-429 contributed to the observed lung phenotype in miR-200b deficient mice. Moreover, miR-200c and miR-141, which are transcribed independently, are expressed normally suggesting that there are no compensatory effects.
Our miR-200b−/− mice did not show any obvious breathing problems after birth, but lung function studies demonstrated severe peripheral airway obstruction in 8-week-old miR-200b−/− mice. These observed lung function abnormalities can be due to lung surfactant deficiency or lung fibrosis. Using human fetal lung cultures, Benlhabib, et al. showed previously that miR-200 family members regulate epithelial type II cell differentiation and function33. They found that miR-200 family inhibitors down-regulated pro-SP-C and SP-B expression. Here, we show that miR-200b−/− lungs have dysfunctional surfactant and compromised expression of these two proteins.
Using in vivo micro-CT scanning on alive animals, we found higher parenchymal density with significantly less air-filled distal alveoli in miR-200b−/− mice. Our lung morphometry studies corroborated these results by showing that the percentage of airspace in miR-200b−/− lungs was lower than in miR-200b+/+ lungs. Thus, both studies suggest that the observed peripheral airway obstruction is due to thicker alveolar walls reducing the airspace inside miR-200b−/− lungs. Our immunofluorescence studies indicate that Vimentin and Twist expression are upregulated in the peripheral lung tissues of miR-200b−/− mice. Interestingly, this was not confirmed in our NGS studies of miR-200b−/− lungs. It has been shown that regulation of vimentin and twist expression are mainly at the translational level34, 35. Our data suggests that miR-200b is involved in the translation regulation of these genes in the lung. Taken together, the observed lung function abnormalities in miR-200b−/− mice can be due to a mesenchymal-skewed, “fibrosis”-like lung phenotype. Others have suggested a role for miR-200b in lung fibrosis before21.
Using NGS analysis, we identified changes in the transcriptome in miR-200b−/− lungs. Of note, we did not observe upregulation of the primary direct targets of the miR-200 family in the transcriptome analysis of miR-200b−/− lungs. We hypothesize that this can be explained by the fact that most microRNAs directly regulate gene expression at the translational level. Therefore, changes in proteins of the direct gene targets of the miR-200 family might not be reflected in the transcriptome of miR-200b−/− lungs. Gene ontology analysis showed that different signaling pathways are affected by miR-200b absence.
Cyp2a5 (cytochrome P450, family 2, subfamily a, polypeptide 5) mRNA is one of the most down-regulated mRNAs in miR-200b−/− lungs. Cytochrome P450s comprise a superfamily of enzymes crucial for metabolism of a diverse group of compounds, drugs and environmental pollutants. These enzymes are mainly present in the liver, but lower levels have been detected in the lung36. Down-regulation of CYP2A5 suggests the involvement of miR-200b in lung metabolism.
Palate lung and nasal epithelial clone (Plunc) is one of the mRNAs that was down-regulated more than four times in miR-200b−/− lungs. PLUNC is the most abundant secretory protein in the lung and is expressed in nasal, oropharyngeal, and lung epithelial cells. Plunc acts as airway surfactant and plays a role in pulmonary host defence37, 38. It is essential for maintaining normal airway surface liquid homeostasis39, 40. Lower levels of Plunc can explain the increased elastance and airway resistance observed in the miR-200b−/− lungs.
Cadherin-26 (CDH26) is another significantly down-regulated mRNA in the miR-200b−/− lungs. Recently, others showed that CDH26 is involved in regulating lung epithelial cell polarity and differentiation. Knockdown of CDH26 results in a lack of epithelial cell polarity and differentiation41. Down-regulation of CDH26 in our knockout lungs can explain the fibroblast-like phenotype of the lungs in our knockout mice and the role of miR-200b in maintaining an epithelial cell phenotype.
Interestingly, these lungs still express the epithelial marker cytokeratin normally. These findings suggest a “partial EMT” in the lung parenchyma of miR-200b deficient mice. Others have shown that patients with idiopathic pulmonary fibrosis (IPF) have epithelial basal cells with a partial mesenchymal phenotype surrounding the fibroblastic foci42. Members of the miR-200 family are down-regulated in the lungs of patients with IPF and a mouse model of lung fibrosis21. We have previously discovered that lung hypoplasia in CDH lungs is characterized by reduced airway branching and a thickened interstitial mesenchymal cell layer recapitulating a fibrotic lung phenotype4, 43.
Knocking down miR-200b in our mice resulted in down-regulation of mature miR-200a and miR-429. These two miR-200 family members share the same transcript with miR-200b. This down-regulation can result from the removal of the miR-200b gene on the promotor activity or can be a post transcriptional effect of miR-200b on processing of the other microRNAs. Although miR-200c has the exact same seed sequence as miR-200b, the abundance of this microRNA along with miR141 was not changed in lungs and kidneys of 8-week old miR-200b−/− mice suggesting the absence of any compensatory effects of these two miR-200 family members in miR-200b deficient mice. However, it is unclear at this point what the influence of downregulation of miR-200a and miR-429 on the lung phenotype of miR-200b deficient mice is.
Our in vivo studies demonstrate a new role for miR-200b in distal lung airway development by regulating epithelial and fibroblast cell differentiation.
Mice were maintained in accordance with the guidelines of “Guide to the Care and Use of Experimental Animals”. The in vivo experiments were approved by the Animal Research Review Committee at the University of Manitoba.
We generated a miR-200b−/− (KO) mouse by replacing the complete miR-200b gene by targeted homologous recombination in C57Bl/6N mouse embryonic stem cells using the NorCOMM knockout cassette (online data supplement).
Whole mouse embryos are stained for β-galactosidase (lacZ) activity using X-gal as described before44. For E14.5 and older, we used a razor blade to section the embryos in half to facilitate penetration.
Eight-week-old male miR-200b+/+(wt) or miR-200b−/− (ko) mice (at least six mice for each group) were anesthetized with intra-peritoneal sodium pentobarbital (90mg/kg). The trachea was dissected using fine dissection scissors and a 20-gauge polyethylene catheter was inserted which was further connected to a flexiVent small animal ventilator (Scireq Inc. Montreal). Mice were ventilated with a tidal volume of 10ml/kg body weight, 150 times per minute. A positive end expiratory pressure (PEEP) of 3 cmH2O was used for all studies. Mice were subjected to an increased dose of nebulized methacholine (MCh) challenge protocol to assess concentration response characteristics of respiratory mechanics. For MCh challenge, ~30μL of saline containing from 0 to 50mg/ml MCh was delivered over 10seconds using an in-line ultrasonic nebulizer. To assess the effects of MCh challenge on respiratory mechanics we used low frequency forced oscillations (1–20Hz). Respiratory mechanical input impedance (Zrs) was derived from the displacement of the ventilator’s piston and the pressure in its cylinder. Correction for gas compressibility, and resistive and accelerative losses in ventilator, tubing and catheter were performed per the manufacturer instructions, using dynamic calibration data obtained from snap shot perturbation applied to the system in an open and closed configuration. By fitting Zrs to the constant phase model and Prime-3 perturbation flexiVent software calculated conducting airway resistance known as Newtonian resistance (Rn), peripheral tissue/airway resistance known as tissue damping (G) and tissue elastance or stiffness (H). Values for each parameter were calculated as the mean of all 12 perturbation cycles performed after each MCh challenge.
Bronchoalveolar lavage fluid (BALF) was collected with 4 repeated washes of excised lungs using 2ml saline in total. The fresh supernatant was used for assessing biophysical surfactant function using a capillary surfactometer according to the manufacture’s protocol (Calmia Medical, Inc.). We used 3 miR-200b−/− male mice and 3 miR-200b+/+ for this measurement and each sample was measured 5 times.
Three miR-200b−/− and three miR-200b+/+ male mice were anesthetized in an anesthetic chamber with 5% isoflurane and scanned for 32minutes, while the mice were breathing normally under an anesthetic mask, using the SkyScan 1176 x-ray microtomography system equipped with a large format 11 megapixel x-ray camera (Small animal model imaging core facility, University of Manitoba). Images were acquired at 18μm resolution with an exposure time of 310ms and 0.5° rotation step using a 0.5mm aluminum filter and source current and voltage of 500μA and 50kV respectively. Images were reconstructed using NRecon (Bruker MicroCT, Kontich, Belgium) with the dynamic range set to the same values for each mouse (0–0.08) and a beam hardening correction of 30%. Reconstructed images were processed using CTan (Bruker MicroCT, Kontich, Belgium) to visualize airways. To ensure that the same region of tissue was used for comparison between mice, the slice showing initial branching of the trachea into bronchi was found and the first slice for analysis was set to be 3.25mm (181 slices) above this slice. The last slice included in the analysis was set to the slice corresponding to the base of the lungs (i.e. the first slice showing only diaphragm and no lung tissue). Image size was reduced by drawing a region of interest around the lung tissue, excluding muscle and bone, and saving the reduced volume within the region of interest between the upper and lower slice limits. Slices were then loaded in CTVox (Bruker MicroCT, Kontich, Belgium) as minimum intensity projection images. The reduced data set was further processed in CTan to allow calculation of airspace volume and structure thickness. The 8-bit images were binarised using a lower threshold of 15 and an upper threshold of 255 to segment tissue containing air from surrounding tissues. An ROI shrink-wrap was performed in 2D followed by the bitwise operation image=(image) XOR (ROI) to produce an image of the air-containing tissue. This image was then de-speckled to remove white speckles smaller than 15 voxels. The volume of the segmented regions was calculated to produce an estimate of the volume of air in the lungs.
Three miR-200b−/− and three miR-200b+/+ 8-week-old lungs were inflation-fixated and embedded in paraffin. After performing H&E staining on three sections per lungs, we scanned the whole sections using an Axio Scan.Z1 microscope from Zeiss. We measured area percentage of the airspace using ZEN Image Analysis software and a module based on color coding the tissue and empty areas (airspace). Using the software, we eliminated the large airways from the calculation (Supplementary Fig. 7).
Three 8-week-old miR-200b−/− or miR-200b+/+ lungs were inflated intratracheally with 1ml of 10% formalin via the cannula by gentle infusion. All tissues were fixed in 10% formalin for overnight. Lungs were embedded in paraffin after dehydration in a graded ethanol series followed by xylene. Immunofluorescence was performed on 5μm sections using vimentin (Abcam; ab92547; 1:800), pro-Surfactant Protein-C (Abcam; ab40879; 1:800), pro+mature Surfactant Protein-B (Abcam; ab40876, 1:50) and twist (Abcam; ab50887; 1:200) antibodies as previously described45.
Human bronchial epithelial cells, BEAS-2B (ATCC, Manassas, VA, USA), cells were transfected with 0.01μg/ml of LNA-hsa-miR200b inhibitor or negative control oligonucleotides Exiqon, Denmark). After 48h cell were fixed with 4% PFA and immunofluorescence was performed using vimentin (Abcam; ab92547; 1:800) and cytokeratin (BioLegend, 628602, 1:200).
For mice lung explant culture, lungs were isolated from E11.5 embryos (offspring from a miR-200b +/− cross) and transferred to porous membranes (IsoporeTM) filters with dimensions of 1mm×1.5mm pore size (Millipore, USA) in a 12-well plate for a semidry floating explant culture and cultured for four days in a 1:1 mixture of DMEM and Ham’s F-12 Nutrient supplemented with 100μg/ml streptomycin, 100units/ml penicillin, 0.25mg/ml ascorbic acid. Branching morphogenesis and epithelial perimeter length were monitored daily in all groups by stereomicroscopy, photographs taken and measurements performed using ImageJ software. The difference between day 0 (D0: 0hours) and day 4 (D4: 96hours) of culture, were expressed as D4/D0 ratio.
The scratch (wound healing) or migration assay was performed as previously described by others46. Briefly, BEAS-2B cells were plated in 12 well plates and transfected with 0.01μg/ml of LNA-hsa-miR200b inhibitors or negative control oligonucleotides. After 12hours, we scraped the cell monolayer in a straight line with a 200-µl pipette tip. The cell cultures were washed with PBS to remove the debris followed by transfection with miR-200b inhibitors. Cell migration was monitored and photographed over time. Distance travelled from the initial scratch site was measured after 18h and the migration distance was quantified with ImageJ software.
The library preparation was done using the TruSeq® Stranded mRNA Sample preparation kit (Illumina inc). The starting material of total RNA (100ng) was mRNA enriched using the oligodT bead system. The isolated mRNA was subsequently fragmented using enzymatic fragmentation. Then first strand synthesis and second strand synthesis were performed and the double stranded cDNA was purified (AMPure XP, Beckman Coulter). The cDNA was end repaired, 3′ adenylated and Illumina sequencing adaptors were ligated onto the fragment ends, and the library was purified (AMPure XP). The mRNA stranded libraries were pre-amplified with PCR and purified (AMPure XP). The libraries size distribution was validated and quality inspected on a Bioanalyzer high sensitivity DNA chip (Agilent Technologies). High quality libraries were quantified using qPCR, the concentration normalized and the samples pooled according to the project specification (number of reads). The library pool(s) were re-quantified with qPCR and optimal concentration of the library pool used to generate the clusters on the surface of a flowcell before sequencing using Nextseq500/ High Output sequencing kit (51 cycles according to the manufacturer instructions (Illumina Inc.). All experiments were conducted at Exiqon Services, Denmark.
All data are presented as mean +/− standard error of mean, from a minimum of three independent experiments. Statistical significance was determined by two-way ANOVA or by student’s t-test as indicated in the figure legends. A p-value≤0.05 was considered significant.
We would like to thank Songyan Liu for bioinformatics analysis and Luke DeLange for assistance with vector construction. This work was supported by the Canadian Institutes of Health Research; the Manitoba Medical Services Foundation; the Children’s Hospital Research Institute of Manitoba; the Molly Towell Perinatal Research Foundation; Department of Surgery GFT Surgeons; the Thorlakson Foundation; RK is the recipient of a Career Development Award from the Canadian Child Health Clinician Scientist Program and a New Investigator Salary Award from the Canadian Institutes of Health Research, Manitoba Lung Association and the Children’s Hospital Research Institute of Manitoba. Canadian Institutes of Health Research and the Manitoba Research Innovation Fund supports GGH. AJH is supported by the Canada Research Chair Program.
N.K., R.V., L.F., C.D., D.A., B.I., F.Z., A.O., S.B., M.P., A.F., M.J., V.S., G.S., G.H., A.H. and R.K. acquired and analyzed the data for this study. N.K., R.V., A.O., G.H., A.H. and R.K. delineated the hypothesis and designed the studies. N.K. and R.K. wrote the manuscript. All authors reviewed and revised the manuscript.
R.K. holds a patent application (PCT/CA2015/051028) containing technology described in the manuscript.
Electronic supplementary material
Supplementary information accompanies this paper at doi:10.1038/s41598-017-05412-y
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.