Skip to main content
Download PDF

Lung branching morphogenesis is accompanied by temporal metabolic changes towards a glycolytic preference

Cell & Bioscience volume 11, Article number: 134 (2021) Cite this article

Abstract

Background

Lung branching morphogenesis is characterized by epithelial-mesenchymal interactions that ultimately define the airway conducting system. Throughout this process, energy and structural macromolecules are necessary to sustain the high proliferative rates. The extensive knowledge of the molecular mechanisms underlying pulmonary development contrasts with the lack of data regarding the embryonic lung metabolic requirements. Here, we studied the metabolic profile associated with the early stages of chicken pulmonary branching.

Methods

In this study, we used an ex vivo lung explant culture system and analyzed the consumption/production of extracellular metabolic intermediates associated with glucose catabolism (alanine, lactate, and acetate) by 1H-NMR spectroscopy in the culture medium. Then, we characterized the transcript levels of metabolite membrane transporters (glut1, glut3, glut8, mct1, mct3, mct4, and mct8) and glycolytic enzymes (hk1, hk2, pfk1, ldha, ldhb, pdha, and pdhb) by qPCR. ldha and ldhb mRNA spatial localization was determined by in situ hybridization. Proliferation was analyzed by directly assessing DNA synthesis using an EdU-based assay. Additionally, we performed western blot to analyze LDHA and LDHT protein levels. Finally, we used a Clark-Type Electrode to assess the lung explant's respiratory capacity.

Results

Glucose consumption decreases, whereas alanine, lactate, and acetate production progressively increase as branching morphogenesis proceeds. mRNA analysis revealed variations in the expression levels of key enzymes and transporters from the glycolytic pathway. ldha and ldhb displayed a compartment-specific expression pattern that resembles proximal–distal markers. In addition, high proliferation levels were detected at active branching sites. LDH protein expression levels suggest that LDHB may account for the progressive rise in lactate. Concurrently, there is a stable oxygen consumption rate throughout branching morphogenesis.

Conclusions

This report describes the temporal metabolic changes that accompany the early stages of chicken lung branching morphogenesis. Overall, the embryonic chicken lung seems to shift to a glycolytic lactate-based metabolism as pulmonary branching occurs. Moreover, this metabolic rewiring might play a crucial role during lung development.

Background

Pulmonary branching morphogenesis is an intricate process governed by epithelial-mesenchymal interactions and is dependent on complex signaling events. This process occurs throughout the early stages of embryonic lung development and defines the respiratory airway structure [1]. In the chicken, Gallus gallus, the primordial lung appears around day 3 of embryogenesis as a protuberance from the primitive foregut [2]. During this process, the mesobronchus grows distally, and the new secondary bronchi sprout laterally into the surrounding mesenchymal compartment [1, 3]. This lateral or monopodial branching is exceptionally similar to the domain branching subroutine characteristic of the mammalian lung system [4, 5]. Moreover, the molecular events underlying the development of the avian respiratory system are highly conserved with the mammalian and point to similar functions [6]. For instance, FGF (Fibroblast Growth Factor), WNT (Wingless-related Integration Site), SHH (Sonic Hedgehog), and Retinoic Acid signaling pathways were described as playing critical roles in chicken pulmonary branching morphogenesis [7,8,9,10,11].

The signaling mechanisms involved in early lung development are quite well studied in several animal models [10, 12, 13]; however, little is known concerning the embryonic lung metabolic needs [14, 15]. As for the mammalian adult lung metabolic requirements, it has been shown that they are achieved through the uptake and catabolism of glucose, which represents the primary fuel to the adult lung tissue [16,17,18]. The coordination between signaling and metabolism is now emerging as a key concept for understanding developmental processes. Therefore, it is of major importance to investigate how metabolism contributes and is dynamically regulated during animal development [19, 20].

Recent studies have demonstrated that glycolysis can serve additional roles beyond the classical bioenergetics purpose and contribute to shaping embryonic development both in time and space. For instance, during the murine chorioallantoic branching stage, the developing embryo redirects glucose carbon flow into the pentose phosphate pathway by suppressing phosphofructokinase 1 (PFK1) and aldolase; concomitantly, there is an increase in the glycolytic fraction to serve lactate biosynthesis. This study suggests a rewiring of glycolytic metabolism in the whole embryo, and over time, to promote accurate chorioallantoic branching [21]. Moreover, Slaninova et al. showed that loss of NOTCH signaling in Drosophila and human microvascular cells resulted in the downregulation of pivotal glycolytic genes; in opposition, a short pulse of NOTCH signaling stimulates glycolysis and increases lactate production, in detriment of Tricarboxylic Acid Cycle (TCA) activity [22]. Overall, these results revealed that NOTCH signaling promotes a glycolytic shift that resembles the Warburg effect [22]. Furthermore, Drosophila larval neuroblast cells, while proliferating, depend on aerobic glycolysis. However, during cell cycle exit and terminal differentiation, cell metabolism switches from aerobic glycolysis to an OXPHOS-based metabolism [23]. More recently, Bulusu et al. described a glycolytic activity gradient that contributes to Presomitic Mesoderm (PSM) differentiation and development in the mouse embryo [24]. Additionally, it has been shown that there is an FGF/Wnt coordinated glycolytic gradient that regulates cell motility and controls specification, thus contributing to PSM development in the chicken embryo [25].

This study aimed to investigate the metabolic profile underpinning the early stages of pulmonary branching morphogenesis. Chicken embryonic lungs were used to perform ex vivo lung explant culture. Extracellular metabolites associated with glucose catabolism (alanine, lactate, and acetate) were evaluated by 1H Nuclear Magnetic Resonance (1H-NMR) spectroscopy in the culture medium. The expression patterns/levels of key enzymes and transporters from the glycolytic pathway were assessed by in situ hybridization and qPCR, namely: glucose transporters (glut1, glut3, and glut8), monocarboxylate transporters (mct1, mct3, mct4, and mct8), and metabolic-related enzymes such as hexokinase (hk1 and hk2), phosphofructokinase (pfk1), lactate dehydrogenase (ldha and ldhb) and pyruvate dehydrogenase (pdha and pdhb). Proliferation status was determined by directly assessing DNA synthesis using an EdU-based assay. Additionally, the protein expression levels of LDH were determined by Western blot. Lastly, the respiratory capacity of explants was evaluated by measuring the basal oxygen consumption rate (OCR).

This report describes, for the first time, the temporal metabolic changes that accompany the early stages of chicken lung branching morphogenesis. In our experimental setting, the embryonic chicken lung shows a glycolytic preference with a shift to lactate production as pulmonary branching proceeds.

Results

Lung branching morphogenesis is accompanied by temporal metabolite changes

To describe the metabolic alterations that occur during the early stages of chicken pulmonary branching, in vitro lung explant culture was performed using stages b1, b2, and b3 (1, 2, or 3 secondary buds formed per bronchus, respectively) that correspond to the first three branching stages (Fig. 1a); new branches are clearly seen with the epithelial marker l-cam (Fig. 1a). Explant culture allows assessing extracellular metabolite fluctuations while in a controlled environment. Briefly, the culture system was performed for 48 h and refreshed at D1 (24 h); medium was collected at D0 (0 h; to be used as reference/control), D1, and D2 (48 h), and then analyzed by 1H-NMR spectroscopy. Metabolite production/consumption was calculated following the mathematical formula | (D1-D0) + (D2-D0) |, expressed in pmol, and normalized to the total amount of protein. From the 1H-NMR spectra analysis, it was possible to detect the following metabolites: glucose, alanine, lactate, and acetate (Fig. 1b). Additionally, morphometric analysis was performed to assess branching, revealing a gradual increment of the epithelial compartment when comparing between stages (B1 vs B2 vs B3), which implies an increase in lung branching morphogenesis (Additional file 1: Figure S1).

Fig. 1
figure 1

Extracellular metabolite fluctuations throughout early stages of lung branching. a Representative examples of embryonic chicken lungs at DO (0 h) (b1, b2, b3) and D2 (48 h) (B1, B2, B3) of explant culture, and corresponding in situ hybridization for l-cam, (n = 5/stage). l-cam is an epithelial marker, showing active branching formation. Scale bar: 500 µm. b Schematic representation of the glucose catabolism pathway and the potential pyruvate destinations. Blue labeling indicates the metabolites that were detected and quantified in the 1H-NMR spectroscopy analysis. c Glucose consumption, d Alanine production, e Lactate production, and f Acetate production, during 48 h of explant culture. The medium was refreshed at D1 (24 h). Medium samples were collected at D1 (24 h) and D2 (48 h). D0 media samples were used as control. Metabolite consumption or production was calculated following the mathematical formula | (D1-D0) + (D2-D0) |. Metabolite data was normalized to the total amount of protein. Results are expressed as | mean |± SEM (n ≥ 5/stage). One-Way ANOVA and Fisher’s LSD test were performed. Significantly different results are indicated as: *p < 0.05; **p < 0.01

Full size image

The 1H-NMR results revealed that glucose consumption, by the embryonic lung, decreases from B1 (3.8 × 107 ± 1.7 × 106 pmol/mg protein) to B2 (3.2 × 107 ± 2.0 × 106 pmol/mg protein) (p < 0.05) and from B1 to B3 (3.3 × 107 ± 4.7 × 105 pmol/mg protein) (p < 0.05) (Fig. 1c). Alanine is produced (Fig. 1d), revealing a statistically significant increase from B1 (1.3 × 106 ± 4.6 × 105 pmol/mg protein) to B3 (2.9 × 106 ± 2.2 × 105 pmol/mg protein) (p < 0.05). Likewise, a progressive upsurge in lactate production is detected (Fig. 1e). Indeed, the production of lactate during branching increases from stage B1 (1.6 × 107 ± 2.3 × 106 pmol/mg protein) and from stage B2 (1.9 × 107 ± 1.8 × 106 pmol/mg protein) to B3 (2.5 × 107 ± 9.7 × 105 pmol/mg protein) (p < 0.01 and p < 0.05, respectively). Similarly, acetate production (Fig. 1f) increases through lung branching stages, namely from B1 (2.1 × 106 ± 3.0 × 105 pmol/mg protein) to B3 (3.2 × 106 ± 1.6 × 105 pmol/mg protein) (p < 0.05).

Key glycolytic enzymes and transporters are present in the embryonic chicken lung

To characterize the molecular machinery underlying the metabolic variations associated with early branching stages, key enzymes and transporters involved in glucose catabolism were evaluated by qPCR (Fig. 2a). For this purpose, b1 to b3 lung explants were collected at 0 h and after 48 h of culture (Fig. 1a) and assessed for the expression of glut1, glut3, glut8, mct1, mct3, mct4, mct8, hk1, hk2, pfk1, ldha, ldhb, pdha, and pdhb.

Fig. 2
figure 2

mRNA expression levels of glucose catabolism transporters in early stages of lung branching. a Schematic representation of the glucose catabolism pathway. Blue labeling indicates the enzymes and transporters that were evaluated by qPCR. Relative expression levels of b glut1, c glut3, d glut8, e mct1, and f mct8 in embryonic lungs at D0 (0 h) (b1, b2, b3) and D2 (48 h) (B1, B2, B3) of explant culture. g Expression fold variation, from 0 to 48 h of explant culture, for glut1, glut3, glut8, mct1, and mct8. mRNA expression levels were normalized for both 18 s and actin-β housekeeping genes. Results are expressed in arbitrary units, as mean ± SEM (n ≥ 5/stage/condition). One-Way ANOVA and Fisher’s LSD test were performed. Significantly different results are indicated as: *p < 0.05; **p < 0.01

Full size image

glut1, glut3, and glut8 transcripts (Fig. 2b–d) are present in all three pulmonary stages. glut1 displays an increase in the expression levels from b1 (0 h) to B1 (48 h) with p < 0.01 (Fig. 2b). Similarly, glut3 expression levels (Fig. 2c) increase between b1 and B1 (p < 0.05), and also between b2 and B2 (p < 0.05). After 48 h, both glut1 and glut3 expression levels decrease from B1 to B3 (p < 0.01; p < 0.05). Likewise, glut8 transcript (Fig. 2d) decreases from B2 to B3 (p < 0.05).

mct1 and mct8 transcripts (Fig. 2e, f) display a similar tendency in the embryonic chicken lung. Both mct1 and mct8 exhibit an increase from b1 to B1 (p < 0.05; p < 0.01) and then a decrease from B1 to B3 (p < 0.05; p < 0.01). After 48 h, MCT transcripts behave comparably to GLUTs. mct3 and mct4 expression levels were virtually undetectable in the embryonic chicken lung (data not shown).

hk1 transcript (Fig. 3a) is present in the three pulmonary stages studied and displays statistically significant differences between b2 and B2 (p < 0.05). Also, hk1 expression decreases from B2 to B3 (p < 0.05). hk2 mRNA levels are practically null in the embryonic chick lung (data not shown).

Fig. 3
figure 3

mRNA expression levels of glucose catabolism enzymes in early stages of lung branching. Relative expression levels of a hk1, b pfk1, c ldha, d ldhb, e pdha, and f pdhb, in embryonic lungs at D0 (0 h) (b1, b2, b3) and D2 (48 h) (B1, B2, B3) of explant culture. g Expression fold variation, from 0 to 48 h of explant culture, for hk1, pfk1, ldha, ldhb, pdha, and pdhb. mRNA expression levels were normalized for both 18 s and actin-β housekeeping genes. Results are expressed in arbitrary units, as mean ± SEM (n ≥ 5/stage/condition). One-Way ANOVA and Fisher’s LSD test were performed. Significantly different results are indicated as: *p < 0.05; **p < 0.01; ***p < 0.001

Full size image

pfk1 expression levels are maintained between groups (Fig. 3b), except for a statistically significant increase between B2 and B3 (p < 0.05).

Concerning ldha transcript (Fig. 3c), an increase in the expression levels from b1 to B1 (p < 0.001) and from b3 to B3 is detected (p < 0.05). After 48 h, both B2 and B3 lungs display a decrease in ldha levels compared to B1 lungs (p < 0.05). ldhb (Fig. 3d) is present in all the stages and time points analyzed, and statistically significant differences were observed only between b1 and B1 lungs (p < 0.05).

pdha expression levels (Fig. 3e) are moderately maintained among the groups studied. On the opposite, pdhb transcript levels (Fig. 3f) increase from b1 to B1 (p < 0.05) but decrease between B1 and B3 lungs (p < 0.01).

For all transcripts, fold variation was expressed as the change that occurs from 0 to 48 h of explant culture (Figs. 2g and  3g). Major variations were observed in stages b1 and b2 for glucose transporters (glut1 and glut3) and monocarboxylate transporters (mct1 and mct8) (Fig. 2g). Still, the b3/B3 stages exhibited only minor variations in the same groups of transcripts (Fig. 2g). hk1 and ldha showed a potential synergism with symmetric behavior regarding the fold variation (Fig. 3g). ldhb revealed a fold variation decrease for all stages, from 0 to 48 h. Still, no major expression variations were observed in pfk1, pdha, and pdhb (Fig. 3g).

ldha and lhdb exhibit region-specific expression patterns in the embryonic lung

To study the spatial distribution of lactate dehydrogenase, the expression pattern of ldha and ldhb was characterized by in situ hybridization. Afterward, representative lungs from the three branching stages were processed for histological sectioning.

ldha mRNA is present in the proximal epithelium of the lung, namely in the trachea region (Fig. 4a, dark arrowhead; 4b, black rectangle). ldha transcript is also expressed in the more distal region of the main bronchus (Fig. 4a, dagger). However, it is not expressed in the main bronchus epithelium and secondary bronchi (Fig. 4c, black arrow and asterisk, respectively). This expression pattern is conserved among the three branching stages but decreases throughout branching morphogenesis, namely for stage b3. Slide sectioning of the hybridized lungs confirmed the presence of ldha in both mesenchymal and epithelial compartments of the embryonic trachea (Fig. 4d).

Fig. 4
figure 4

ldha and ldhb mRNA expression pattern at early stages of chick lung branching. Representative examples of in situ hybridization of stage b1, b2, and b3 lungs for a-d ldha and e–h ldhb, n ≥ 9 per stage. Scale bar: whole mount, 500 μm; slide sections, 100 μm. The black rectangle in images b and f indicates the region shown in the corresponding slide section. Asterisk: secondary bronchi. Black arrow: main bronchus epithelium. Dagger: distal region. Dark arrowhead: trachea region

Full size image

ldhb transcript is completely absent from the proximal region of the lung (Fig. 4e, dark arrowhead) and the epithelium of the primary bronchus (Fig. 4f, black arrow). ldhb is strongly expressed in the distal-most region of the lung (Fig. 4e, dagger) and the secondary bronchi (Fig. 4g, asterisk). This expression pattern is maintained in the three stages studied. Lung sectioning confirmed that ldhb is not expressed in the distal epithelial compartment but exclusively present in the distal mesenchyme of the growing tips (Fig. 4h, dagger).

Lung active branching sites are associated with high proliferation

To assess the proliferation status of lung branching morphogenesis, we performed an EdU-based proliferation assay using B1 to B3 lung explants (Fig. 1a). Proliferation was determined by directly assessing EdU incorporation into new DNA strands using Alexa Fluor 488 (Green); nuclei were counterstained with Hoechst 33342 (Red) (Fig. 5).

Fig. 5
figure 5

Proliferation analysis of lung branching morphogenesis. Representative confocal microscopy fluorescence images of lungs explants after 48 h in culture (B1, B2, B3). a, d, g Nuclei were stained with Hoechst 33342 (Red). b, e, h Proliferation was assessed using an EdU-based assay; EdU incorporation in DNA was detected using Alexa Fluor 488 (Green). c, f, i Overlay represents the merged images of Hoechst (Red) and EdU (Green). All images are represented as maximum intensity projection of z-stacks. n ≥ 4 per stage. Scale bar: 500 μm. Asterisk: secondary bronchi/active branching sites. White arrow: distal tip. White arrowhead: trachea region

Full size image

After 48 h of culture, high proliferation levels were detected in the trachea region (white arrowhead) (Fig. 5b), in the distal tip of the lung (white arrow) (Fig. 5e), and in the secondary bronchi/active branching sites (asterisk) (Fig. 5h). The proliferation patterns are maintained in the three stages (B1, B2, and B3).

LDH protein levels support the increase of lactate production throughout branching

To understand how LDH protein levels may contribute to the observed increase in lactate production, the protein expression levels of LDHA and total LDH (that recognizes both LDHA and LDHB proteins) were assessed.

As seen in Fig. 6a, at 0 h, LDHA protein levels fluctuate between stage/condition. Semi-quantitative analysis (Fig. 6b) revealed that LDHA expression levels increase at 0 h, namely from b1 to b3 (p < 0.01). On the other hand, after 48 h of culture, this variation dissipates, and there are no differences between stages. Conversely, when comparing both time points for each stage individually, a statistically significant decrease in LDHA expression levels is observed (Fig. 6b), namely from b1 to B1 (p < 0.01), from b2 to B2 (p < 0.0001), and from b3 to B3 (p < 0.0001).

Fig. 6
figure 6

Western blot analysis of LDHA and total LDH throughout chick lung branching. Representative immunoblots for a LDHA and c LDHT (LDHA and LDHB contribution) of pooled-tissue samples of embryonic lungs at D0 (0 h) (b1, b2, b3) and D2 (48 h) (B1, B2, B3) of explant culture. Loading control was performed using β-tubulin (55 KDa). LDHA and LDHT correspond to 37 KDa. Semiquantitative analysis, of three independent experiments, for b LDHA and for d LDHT protein expression levels. Results are represented in arbitrary units, normalized for β-tubulin. Results are expressed as mean ± SEM (n = 3/stage/condition). One-Way ANOVA and Fisher’s LSD test were performed. Significantly different results are indicated as: *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001. Full-length blots are presented in Additional file 1: Figure S2

Full size image

In contrast, total LDH (LDHA plus LDHB) displays similar levels between stages at 0 h (Fig. 6c). Semi-quantitative analysis (Fig. 6d) exposed that there are no statistically significant differences between the three stages analyzed at 0 h. After 48 h of culture, total LDH protein levels vary between stages, namely from B1 to B2 (p < 0.05) and from B1 to B3 (p < 0.01). When comparing both time points for each stage individually, there is a statistically significant increase in total LDH expression levels from b2 to B2 (p < 0.001) and from b3 to B3 (p < 0.0001); in contrast, from b1 to B1 the protein levels remained unaltered.

Basal oxygen consumption rate is maintained during branching morphogenesis

To corroborate that the embryonic lung displays a glycolytic preference as branching morphogenesis develops, the respiratory capacity of explant tissue was evaluated. For this purpose, B1 to B3 lung explant tissue was collected after 48 h of culture (Fig. 1a), and the basal oxygen consumption rate was measured (Fig. 7).

Fig. 7
figure 7

Basal oxygen consumption rate (OCR) measurements throughout early stages of lung branching. Oxygen consumption rate of embryonic lungs after 48 h of explant culture (B1, B2, B3). Results are represented in nmol/min of oxygen consumption normalized for the total amount of protein. Results are expressed as mean ± SEM (n ≥ 5/stage/condition). One-Way ANOVA and Fisher’s LSD test were performed. p < 0.05 was considered

Full size image

After 48 h of culture, basal OCR was maintained between B1, B2, and B3 (Fig. 7). However, a slight decrease tendency in explants’ respiratory capacity throughout pulmonary branching is detected.

Discussion

Lung development is a complex process characterized by epithelial-mesenchymal interactions mediated by signaling cascades, leading to the formation of a fully functional organ [12, 26]. Altogether, these processes contribute to the formation of a fully functional adult organ responsible for gas exchange.

The metabolic profile of the mammalian adult pulmonary system was studied in the ‘70 and ‘80 s, revealing a metabolism that is primarily dependent on glucose utilization [16,17,18]. On the other hand, little is known about embryonic lung metabolic needs. The importance of metabolism during development, beyond its canonical role, has been recently reported; several studies have shown that metabolism is highly regulated both in time and space, affecting multiple aspects of animal development [19]. Moreover, glucose uptake and catabolism are essential for all cellular processes and are believed to be enhanced in high proliferative systems such as cancer and embryonic development [20, 27].

Branching morphogenesis is a crucial step in lung organogenesis since it establishes the future airway conducting system [28]. In the chicken lung, b1, b2, and b3 stages correspond to the sequential appearance of the first three secondary buds, with approximately an 8 h gap between them. These morphologically similar stages revealed to be a useful system to study the molecular mechanisms underlying avian branching morphogenesis. Overall, signaling pathways display similar patterns between stages, and the interaction between different signaling events ultimately shapes airway branching [7,8,9, 11]. In this sense, we wondered if b1, b2, and b3 stages were also similar at the metabolic level or if, despite the molecular similarities, the embryonic lung gradually adapts its metabolic needs to cope with the progressive tissue growth. Accordingly, in this work, we focused on describing the metabolic profile of the embryonic chicken lung during the abovementioned stages of branching morphogenesis. For this purpose, an ex vivo explant culture system was used to precisely measure the metabolic and associated molecular alterations during 48 h of culture [29,30,31]. This method is particularly relevant to study the whole organ while under very controlled conditions since it preserves the epithelial-mesenchymal interactions that contribute to branching morphogenesis [30, 32]. Lung explants are described to maintain the native physiological interactions between cells and tissues, thus mimicking the in vivo structure and function [29]. The culture medium from ex vivo lung explant culture was collected to study the production/consumption of extracellular metabolites related to glucose catabolism (alanine, lactate, and acetate) by 1H-NMR spectroscopy. Afterward, the molecular machinery associated with glucose catabolism (glut1, glut3, glut8, mct1, mct3, mct4, mct8, hk1, hk2, pfk1, ldha, ldhb, pdha, and pdhb) was characterized by qPCR. ldha and ldhb mRNA spatial localization was assessed by in situ hybridization. Proliferation status was determined by directly assessing DNA synthesis using an EdU-based assay. LDH protein expression levels were evaluated by western blot. The respiratory capacity of lung explants was evaluated by measuring the basal oxygen consumption rate.

In this study, 1H-NMR results revealed that glucose consumption by chicken lung explants decreases from B1 to B2 and from B1 to B3 stage, and remains stable between B2 and B3 stage (Fig. 1c). In B1, glucose consumption is higher, and we believe this is due to a higher nutrient requirement to cope with the very beginning of branching morphogenesis. As lung branching morphogenesis proceeds, glucose consumption profile variation may reflect the need to cope with specific energy and nutrients requirements needed to maintain high proliferative rates. Similarly, the adult mammalian lung metabolic requirements are mainly achieved through glucose catabolism, although other metabolic substrates might support lung energetic demands [16,17,18, 33].

Glucose cellular uptake relies on the presence of integral membrane glucose transporters (GLUTs) that move glucose from the extracellular space to the cytosol by facilitated diffusion; furthermore, this process depends on the number of glucose transporters expressed in the cell surface. In this sense, the expression levels of selected transporters, glut1, glut3, and glut8, in the chicken developing lung were characterized. glut1, glut3, and glut8 are expressed in the three pulmonary stages studied (Fig. 2b–d). After 48 h in culture, major variations of glut1 and glut3 are detected in B1 and B2 stages (Fig. 2b, c) but not in the B3 stage (Fig. 2b, c). glut1 and glut3 expression level variations may support the observed glucose consumption variations. glut1 and glut3 transcripts have already been identified during chicken embryonic development in other tissues than the lung, namely, in the brain, muscle, and heart; moreover, glut1 and glut3 expression levels are developmentally regulated in these structures [34]. In the adult, glut1 and glut8 are ubiquitously expressed, whereas glut3 is highly expressed in the brain, similarly to their mammalian counterparts [35,36,37]. In human fetal lung, GLUT1 is present in bronchial and primitive alveolar epithelial cells during the branching phase but becomes progressively less expressed from week 19 onwards [38]. GLUT1 is widely expressed in adult human tissues, and it is commonly overexpressed in several tumor types [39]. On its turn, GLUT3 is highly expressed in lung cancer and actively contributes to cellular glucose uptake and consequent high proliferation of pulmonary tumor cells [40]. Altogether glut1, glut3, and glut8 may likely contribute to sustaining glucose consumption throughout early pulmonary branching morphogenesis stages.

Once inside the cell, glucose is converted into glucose-6-Pi in a reaction catalyzed by hexokinase [41]. hk1 is expressed in the embryonic chicken lung during branching morphogenesis (Fig. 3a), whereas hk2 transcript is barely detected (data not shown).

PFK1 catalyzes a key step of glycolysis, the unidirectional conversion of fructose-6-Pi and ATP into fructose-1,6-Pi and ADP [41]. In the embryonic chicken lung, pfk1 transcript is present in all stages studied without major variations between groups (Fig. 3b). pfk1 expression levels may contribute to control the relative glucose flux into the glycolytic pathway throughout early branching.

From the glycolytic pathway, each molecule of glucose is converted into two molecules of pyruvate. In our explant culture medium, several pyruvate-derived metabolites (alanine, lactate, and acetate from acetyl-CoA) were detected (Fig. 1d–f). During branching morphogenesis, the embryonic chicken lung progressively increases alanine production that is finally exported to the extracellular medium (Fig. 1d). Previous studies revealed that, in the adult mammalian lung, part of the glucose carbons serves for alanine production through pyruvate transamination [17, 18, 33]. Similarly, in prostate cancer, the production of significantly high levels of alanine supports the need for proliferating cells for protein synthesis and membranogenesis [42]. Moreover, in pancreatic ductal adenocarcinoma, alanine produced from the surrounding environment is used by tumor cells for biosynthetic purposes [43]. Thus, alanine is an important end product of glucose catabolism and, in the developing lung, may contribute to protein biosynthesis to support active pulmonary growth.

Lactate was detected at high concentrations in the explant culture medium. Indeed, in our culture conditions, there is a sharp increase in lactate production from stage B1 to B3, and from stage B2 to B3 (Fig. 1e), concurrently with the increase in branching (Additional file 1: Figure S1). The lactate/glucose ratio revealed that in B1 and B2 explants, 41% and 58% of glucose is directed to lactate production, respectively; in B3 explants, around 76% of glucose is converted into lactate. Likewise, in the mammalian adult lung, around 50% of the glucose carbons are metabolized into lactate [16,17,18, 33]. Therefore, lactate production by the developing lung may be facilitating the uptake and incorporation of nutrients by promoting the activity of biosynthetic pathways to, for instance, form new biomass and regenerate NAD+ required for maximal glycolytic flux [44]. The increase in alanine and lactate production is characteristic of high proliferative systems in which fast energy and macromolecules are necessary to sustain the growth.

Lactate dehydrogenase (LDH) is the enzyme responsible for the interconversion of pyruvate into lactate and NADH into NAD+. This enzyme is a tetramer of two types of subunits, LDHA (formerly LDH-M) and LDHB (formerly LDH-H), encoded by ldha and ldhb, respectively [45]. In this work, we characterized the expression levels of both ldha and ldhb in the embryonic chicken lung. ldha expression increases from b1 to B1 and from b3 to B3 (Fig. 3c). Regarding 48 h’ time point, the expression levels are higher at B1 and then decrease substantially to B2 and B3 (Fig. 3c). On the other hand, ldhb expression levels decrease from 0 to 48 h for b1 stage (Fig. 3d). Although the relative expression levels of ldhb are greater than ldha, higher variations accompanying branching were observed in ldha (Fig. 3g). Afterward, we performed in situ hybridization to assess ldha and ldhb spatial localization. ldha and ldhb displayed region-specific expression patterns in the chicken developing lung. ldha is highly expressed in the proximal region of the lung (Fig. 4a and d) and non-existent in the main bronchus epithelium and secondary bronchi (Fig. 4c). In opposition, ldhb transcript is absent from the proximal region of the lung (Fig. 4e) but expressed in the distal region. ldhb is specifically expressed in the growing tips and secondary buds (Fig. 4e and g), suggesting an association with highly proliferative regions. Interestingly, in the chicken and mouse embryo, ldhb is found mostly in the posterior (distal-most) region of the embryo, the tailbud region [25]. ldha and ldhb specific patterns remind a spatial distribution typical of proximal–distal markers.

To determine the proliferation status of lung branching morphogenesis, we performed an EdU proliferation assay, using B1 to B3 lung explants (Fig. 1a). High levels of proliferation were detected in the trachea region (white arrowhead) (Fig. 5b), matching ldha mRNA expression pattern (Fig. 4d). In addition, high levels of proliferation were also present in the secondary bronchi/active branching sites (asterisk) (Fig. 5h), and in the distal tips of the developing lung (white arrow) (Fig. 5e), which coincides with ldhb transcripts spatial localization (Fig. 4e–h). These results point towards an association between ldha and ldhb expression localization and proliferation events occurring during branching morphogenesis.

After the mRNA expression studies, we performed western blot for LDHA and total LDH (LDHA and LDHB contribution). LDHA displayed a decrease between 0 and 48 h (Fig. 6a, b). In contrast, total LDH expression levels increased from 0 to 48 h of ex vivo culture (Fig. 6c, d). From these results, we speculate that LDHB may be contributing to the increase of total LDH protein levels. When we compare mRNA and protein expression data, an inverse pattern between the same ldh/LDH isoform is observed (Figs. 3c vs 6b; Figs. 3d vs 6d). It is generally assumed that mRNA expression levels directly correlate with protein levels. However, several studies have demonstrated that this not always occurs [46, 47]. Distinct factors may account for these differences, such as regulatory mechanisms at transcriptional and/or translational level, or mRNA and protein stability, that can collectively impact total protein levels. In addition, given western blot results and mRNA spatial localization, it is possible that during branching, distal regions produce more lactate to cope with high proliferation rates. The increase in LDH total levels may account for the observed rise in extracellular lactate levels, which seem to be a characteristic of the developing lung.

Since lactate exchange between the cytosol and the extracellular space depends on monocarboxylate transporters (MCTs), we decided to study MCTs in the developing lung. mct1 and mct8 expression levels behave similarly throughout branching (Fig. 2e, f) and recapitulate GLUTs mRNA expression behavior. mct1 displays a considerable increase in stage b1 after 48 h in culture (Fig. 2e), suggesting that it may support intracellular lactate export. Likewise, mct3 and mct4 were evaluated but displayed very low expression levels (data not shown). MCT1 has been described in both human and rat fetal lung; additionally, MCT inhibition leads to a decrease in rat lung branching [38]. More recently, Oginuma and colleagues found that mct1 is expressed, in a graded manner, in the tailbud and posterior PSM of the chicken embryo. Moreover, chicken explants cultured with MCTs inhibitor α-cyano-4-hydroxycinnamic acid (CNCn) displayed increased intracellular lactate levels [48]. On the other hand, cancer cells are equipped with MCTs, mainly MCT1 and MCT4, to exchange lactate and regulate pH homeostasis [49]. MCT4 is the predominant monocarboxylate transporter for lactate export in glycolytic cancer cells; conversely, MCT1 is a passive transporter that can operate in both directions and, in cancer systems, it is described to facilitate lactate export [50]. In fact, MCT1 lactate transport depends on the intracellular/extracellular concentrations of lactate and protons [50]. Chicken mct8, also known as thyroid hormone transporter, was previously described in the developing brain, retina, spinal cord, kidney, and testis of the chicken [51, 52]. Moreover, mct8 knockdown impacts embryonic chicken development, meaning that mct8 expression is crucial from the very early stages of organogenesis [53]. However, MCT8 does not transport lactate, but it may transport other monocarboxylates and amino acids [54]. In the embryonic chicken lung, mct1 can be the transporter responsible for cellular lactate exchange.

Pyruvate can be transformed into acetyl-CoA through an irreversible enzymatic reaction catalyzed by pyruvate dehydrogenase (PDH) and then incorporated in the TCA cycle or converted into acetate [41]. Both isoforms of pyruvate dehydrogenase, pdha and pdhb, are expressed in the embryonic chicken lung (Fig. 3e, f). In the embryonic lung explant culture medium, we did not find detectable levels of metabolites from the Krebs cycle in the 1H-NMR spectra. On the other hand, acetate was identified, and its production is increased in later stages of early pulmonary branching (Fig. 1f). Altogether, these data suggest that part of pyruvate may be diverted to acetate, through acetyl-CoA, as pulmonary branching proceeds. In mammals, glucose-derived pyruvate can generate acetate, and this phenomenon is more pronounced when under conditions of hyperactive glucose metabolism [55]. Moreover, in the adult mammalian lung, glucose carbons can incorporate acetate for lipogenesis purposes and serve for surfactant production [17, 18, 33, 56]. Thus, in the embryonic lung, acetate may also point to a hyperactive glucose metabolism and can be produced to incorporate the newly synthesized cellular membranes, a common mechanism of high proliferative cells.

Lastly, to assess whether the embryonic lung truly shows a glycolytic preference during branching morphogenesis, mitochondrial respiration was evaluated. For this purpose, the respiratory capacity of explants was assessed by measuring the basal oxygen consumption rate of the embryonic lungs (Fig. 7). After 48 h of culture, explants display a steady basal OCR in the three stages; nonetheless, there is a slight tendency to decrease as branching morphogenesis proceeds (Fig. 7). This phenomenon occurs in parallel to the progressive increase in alanine, lactate, and acetate production, previously described (Fig. 1d–f). Altogether, these results suggest that, under aerobic conditions and with functional mitochondria, the embryonic lung seems to shift to a glycolytic lactate-based metabolism throughout branching morphogenesis.

Conclusions

This study describes the temporal metabolic changes that occur during early chicken pulmonary branching. Throughout this period, we observed changes in the metabolite profile that occur concurrently with variations in the expression levels of pivotal enzymes and transporters from the glycolytic pathway and with a steady oxygen consumption rate (Fig. 8). It appears that pulmonary branching morphogenesis progressively adapts to a glycolytic lactate-based metabolic profile, suggesting a Warburg-like metabolism. This metabolic rewiring to aerobic glycolysis is also observed in other high proliferative developing systems and allows the production of energy and biomass during embryonic development [20,21,22, 24, 25, 57]. Moreover, this metabolic adaptation might exert additional signaling functions throughout pulmonary branching, and further studies are required to address this topic. We acknowledge that this study has a strong descriptive perspective, but we are addressing a subject that has not been exploited so far. This report highlights the importance of metabolic regulation during the early stages of lung development and lays the groundwork for future mechanistic studies on this topic.

Fig. 8
figure 8

Schematic representation of the metabolic profile of early pulmonary branching morphogenesis. The metabolites detected in the extracellular medium (blue labeling) can either be consumed (glucose) or produced (alanine, lactate, and acetate) by the lung tissue. The symbols refer to increase (↑) or decrease (↓) in the metabolites, considering their particular fate. Temporal metabolite changes are accompanied by variations in the expression of glut1, glut3, glut8, hk1, pfk1, mct1, mct8, ldha, ldhb, pdha, and pdhb. Decreased LDHA and increased LDHB protein levels may contribute to the observed increase in the extracellular levels of lactate. OXPHOS is active with mitochondria displaying a constant rate of oxygen consumption. The symbols refer to: does not contribute to lactate transport (*); oxidative phosphorylation is active (✓)

Full size image

Methods

Ethical statement

This work was performed at the early stages of chicken embryonic development and, therefore, does not require ethical approval following the European Parliament Directive 2010/63/EU of 22 September 2010 and the Portuguese Directive 113/2013 of 7 August 2013 on the protection of animals used for scientific purposes.

Tissue collection

Fertilized chicken eggs, Gallus gallus, were incubated between 4.5 and 5.5 days (Embryonic day 4.5–5.5) in a 49% humidified atmosphere at 37 °C (Termaks KB400, Norway). Embryonic lungs were obtained by microdissection under a stereomicroscope (Olympus SZX16, Japan) and classified in stage b1, b2, or b3 according to the number of secondary buds formed per bronchus; 1, 2, or 3, respectively [7]. Dissected lungs were processed for ex vivo lung explant culture or for in situ hybridization.

Ex vivo lung explant culture

Chicken lung explant culture was performed, as previously described [30], to study potential temporal metabolic changes while in controlled settings. Briefly, after dissection in PBS, lungs were placed on top of 8 µm nucleopore polycarbonate membranes (Whatman, USA) and incubated for 48 h in 200 µL of medium 199 (5.5 mM glucose; Sigma, USA) supplemented with 10% (V/V) chick serum (Invitrogen, USA), 5% (V/V) heat-inactivated fetal calf serum (Invitrogen), 1% (V/V) L-glutamine (Invitrogen), 1% (V/V) penicillin 5000 IU/mL plus streptomycin 5000 IU/mL (Invitrogen) and 0.25 mg/mL of ascorbic acid (Sigma). The medium was replaced by fresh supplemented medium at 24 h of culture. Lung explants were photographed at 0 h (D0), 24 h (D1), and 48 h (D2) with a camera (Olympus U-LH100HG) coupled to a stereomicroscope (Olympus SZX16). The medium was collected at D0, D1, and D2 for 1H-NMR spectroscopy analysis. D0 and D2 lung explants were collected for RNA and protein extraction; D2 lung explants were collected for EdU proliferation assay and basal oxygen consumption rate assay.

1H-NMR spectroscopy

Samples of 200 µL of medium were collected from ex vivo lung explant culture, at D0, D1, and D2, and analyzed by 1H-NMR spectroscopy (n ≥ 5/stage/condition) according to [58]. A Bruker Avance 600 MHz spectrometer with a 5 mm QXI probe and z-gradient (Bruker Biospin, Germany) was used, and spectra accessed at 25 °C. Solvent-suppressed 1H-NMR spectra were acquired with 6 kHz spectral width, 14-s interpulse, 3-s water presaturation, 45-degree pulse angle, 3.5-s acquisition time, and a minimum of 128 scans. Sodium fumarate (singlet, at 6.50 ppm) was used as an internal reference at 10 mM. The following metabolites were detected and quantified: H1-α glucose (doublet, 5.22 ppm), alanine (doublet, 1.46 ppm), lactate (doublet, 1.33 ppm), and acetate (singlet, 1.9 ppm). To quantify the relative areas of 1H-NMR resonances, the curve-fitting routine supplied with the NUTSproTM NMR spectral analysis program (Acorn NMR, USA) was used. D0 media samples were used as reference/control. The metabolite consumption/production during 48 h of explant culture was calculated following the mathematical formula | (D1-D0) + (D2-D0) |. Metabolite data were normalized to the total amount of protein.

Quantitative PCR

Total RNA was extracted from D0 and D2 lung explants using TripleXtractor directRNA kit (Grisp, Portugal). RNA integrity and purity were determined using the Experion RNA StdSens Analysis Kit (Bio-Rad, USA), and RNA concentration determined by Nanodrop 1000 Spectrophotometer (Thermo Fisher Scientific, USA). Total RNA was treated with DNase I (Thermo Fisher Scientific) and then reversely transcribed to obtain cDNA using GRS cDNA Synthesis kit (Grisp). Specific exon-exon spanning primers were designed for the amplification of target and housekeeping transcripts (Additional file 1: Table S1). Primers were optimized for annealing temperature and PCR cycles, using NZY Taq 2 × Green Master Mix (NZYTech, Portugal) and, subsequently, for efficiency range. qPCR was performed in duplicate (n ≥ 5/stage/condition) using 1 µL of cDNA and SYBR method according to manufacturer’s instructions, NZY qPCR Green Master Mix (2x) (NZYTech). mRNA expression levels were normalized for both 18s and actin-β housekeeping genes. Data on the gene expression levels were calculated following the mathematical model 2^(-ΔCt) [59].

Western blot

Pooled samples of D0 and D2 explants (3 pools per stage: 8 lungs per pool) were processed for western blot analysis, as described in [8]. Protein was obtained according to [60]. 10 µg of protein was loaded onto 10% acrylamide minigels, electrophoresed at 100 V in a Mini-PROTEAN Tetra Cell (Bio-Rad), and transferred to 0.2 µm nitrocellulose membranes using a Trans-Blot Turbo Transfer System (Bio-Rad). Blots were probed with primary antibodies for LDHA (1:20,000; #3582, Cell Signaling, USA), for total LDH that recognizes both LDHA and LDHB (1:20,000; #ab52488, Abcam, UK), and for β-tubulin (1:200,000; #ab6046, Abcam) as a loading control. Afterward, blots were incubated with anti-rabbit secondary horseradish peroxidase-conjugated antibody (1:2000; #7074, Cell Signaling). Membranes were developed with Clarity Max Western ECL substrate (Bio-Rad), and the chemiluminescent signal captured using Chemidoc XRS (Bio-Rad). Quantitative analysis was performed with Image Lab software (Bio-Rad). At least three independent experiments per pool were performed (n = 3).

RNA probes

Total RNA and cDNA from b2 lungs were obtained as previously described. Specific primers were designed for ldha (Fw: 5’-GCACTTTCCAAGTAGGTCAAATCC-3’; Rv: 5’- AGTCTTTGGTTTCACGTTGTGT-3’) and ldhb (Fw: 5’-GCAGGTTGTTGAAAGTGCCT-3’; Rv: 5’-AGTGAGTAGAGAGCCCACAT-3’). cDNA was used as a template for conventional PCR (GRS Xpert Taq 2 × MasterMix, Grisp). PCR fragments were cloned in the pCR™II-TOPO® vector (TOPO TA Cloning Kit, Invitrogen) and sequenced to determine the insert orientation (GATC Biotech, Germany). l-cam RNA probe was produced as previously described [61]. Antisense digoxigenin-labeled RNA probes were produced from PCR-amplified fragments using SP6 or T7 RNA polymerase and according to the manufacturer’s instructions (Roche, Germany).

Whole mount in situ hybridization

Dissected lungs were fixed in PBS solution containing 4% formaldehyde, 2 mM EGTA, pH 7.5, at 4 °C overnight. Then, the lungs were dehydrated in a methanol series and stored at −20 °C. Tissues were rehydrated through a methanol/PBT series and processed for whole mount in situ hybridization (n ≥ 9 per gene/stage) as previously described by [30]. Briefly, tissues were permeabilized with proteinase K solution (PBT with 0.05% proteinase K) (Roche). PBT washes were made to remove the excess of proteinase K and, tissues were incubated with a post-fixing solution (PBT with 10% formaldehyde, 0.4% glutaraldehyde). After, tissues were incubated with hybridization solution (50% formamide, 6.5% SSC, 1% EDTA 0.5 M pH 9.8, 0.25% t-RNA, 0.2% Tween 20, 0.2% heparin, 0.5% CHAPS), at 70 °C. Then, tissues were incubated with specific probes, in hybridization solution, at 70 °C, overnight. On the following day, several washes were performed with preheated hybridization solution, hybridization solution with MABT (50:50) (5.8% C4H4O4, 4.4% NaCl, 7% NaOH, 1% Tween 20, pH 7.5) and subsequently with only MABT. Next, tissues were treated with blocking solutions [MABT with 20% blocking reagent (Roche); MABT with 20% blocking reagent plus 20% goat serum (Invitrogen)]. Then, lungs were incubated in MABT, 20% blocking reagent, 20% goat serum, and 1:2000 anti-digoxigenin antibody (Roche) solution, overnight. On day 3, the tissues were washed with MABT solution. On the last day, tissues were washed in NTMT solution (0.1 M NaCl, 0.1 M Tris–HCl, 50 mM MgCl2, 1% Tween 20) and then incubated in a developing solution (NTMT with BCIP, NBT) (Roche), at 37 °C and protected from light. Each group of lungs/probes was processed simultaneously and developed for the same amount of time. All the lungs were photographed using a camera (Olympus U-LH100HG) coupled to a stereomicroscope (Olympus SZX16).

Histological sections

Hybridized chicken lungs were dehydrated through increasing ethanol series and embedded in a 2-hydroxyethyl methacrylate solution (Heraeus Kulzer, Germany). 25 µm thick histological slides were produced by a rotatory microtome (Leica RM 2155, Germany). Lung histological sections were photographed using a camera (Olympus DP70) coupled to a microscope (Olympus BX61).

Proliferation assay and confocal microscopy

After 48 h of culture, half of the explant’s media was replaced by fresh media containing EdU, at a final concentration of 150 µM. Explants were incubated with EdU for 90 min in the same culture conditions. After incubation, tissues were fixed in PBS solution containing 3.7% formaldehyde. Then, tissues were washed in PBS with 3% BSA, and permeabilized in PBS with 0.5% Triton X-100, for 90 min. Tissues were washed again and processed for Click-iT Plus EdU reaction, according to the manufacturer’s instructions (Click-iT™ Plus EdU Cell Proliferation Kit for Imaging, Invitrogen). Detection of the incorporated EdU was performed using Alexa Fluor 488. Nuclei were counterstained with Hoechst 33342 (1:2000). Images were acquired using an Olympus LPS Confocal FV3000 microscope (Olympus).

Basal oxygen consumption rate measurements

Pooled samples of D2 (2 lungs per pool/per stage) explants were processed to assess respiratory capacity (n ≥ 5/stage/condition). Lung tissue was washed in PBS and incubated, for 2 min, in respiration medium, at 37 °C. Respiration medium was composed by medium 199 (5.5 mM glucose) supplemented with 1% L-glutamine, 1% ITS (0.01 mg/mL recombinant human insulin, 0.0055 mg/mL human transferrin substantially iron-free, and 0.005 μg/mL sodium selenite) (BD Biosciences, USA) and 10 mM HEPES; pH was adjusted to 7.2–7.4, at 37 °C. To assess Oxygen Consumption Rate (OCR), a Clark-type electrode Oxytherm System (Hansatech, UK) was used [62]. Samples were transferred in 1 mL respiration medium to the pre-calibrated and thermostatized (with a water jacket, at 37ºC) electrode chamber. Oxygen recording protocol was performed for 15 min with an open-chamber plus 15 min with a closed-chamber; data acquisition refers to the last 5 min of closed-chamber mode. After each assay, tissue samples were retrieved, washed in PBS, and processed for protein extraction and quantification. Data was analyzed using Oxytrace Plus acquisition software (Hansatech) and normalized to the total amount of protein.

Statistical analysis

Statistical analysis was performed using GraphPad Prism 6 (GraphPad Software, USA). The normality of distribution was tested using the Kolmogorov–Smirnov test. One-Way ANOVA was performed and followed by Fisher’s Least Significant Difference (LSD) post hoc test for multiple comparisons. All experimental data are presented as mean ± standard error of the mean (SEM) with a statistically significant level of 5% (p < 0.05) considered.

Availability of data and materials

Not applicable.

Abbreviations

1H-NMR:

1H Nuclear Magnetic Resonance

BCIP:

5-Bromo-4-Chloro-3-Indolyl Phosphate

BSA:

Bovine Serum Albumin

cDNA:

Complementary DNA

D0:

0 Hours

D1:

24 Hours

D2:

48 Hours

EDTA:

Ethylenediamine Tetraacetic Acid

EGTA:

Ethylene Glycol Tetraacetic Acid

FGF:

Fibroblast Growth Factor

GLUT:

Glucose Transporter

GLUT1:

Glucose transporter 1

GLUT3:

Glucose transporter 3

GLUT8:

Glucose transporter 8

HEPES:

N-2-Hydroxyethylpiperazine-N'-2-Ethanesulfonic Acid

HK1:

Hexokinase 1

HK2:

Hexokinase 2

L-CAM:

Chicken E-cadherin

LDH-H:

Lactate Dehydrogenase Isoform from Heart

LDH-M:

Lactate Dehydrogenase Isoform from Muscle

LDH:

Lactate Dehydrogenase

LDHA:

Lactate Dehydrogenase Isoform A

LDHB:

Lactate Dehydrogenase Isoform B

LDHT:

Lactate Dehydrogenase Total

LSD:

Least Significant Difference

MABT:

Maleic Acid Buffer with Tween 20

MCT:

Monocarboxylate Transporter

MCT1:

Monocarboxylate Transporter 1

MCT3:

Monocarboxylate Transporter 3

MCT4:

Monocarboxylate Transporter 4

MCT8:

Monocarboxylate Transporter 8

mRNA:

Messenger RNA

NAD+ :

Nicotinamide Adenine Dinucleotide

NADH:

Nicotinamide Adenine Dinucleotide

NBT:

4-Nitro Blue Tetrazolium Chloride

OCR:

Oxygen Consumption Rate

OXPHOS:

Oxidative Phosphorylation

PBS:

Phosphate Buffered Saline

PBT:

Phosphate Buffered Saline with Tween 20

PCR:

Polymerase Chain Reaction

PDH:

Pyruvate Dehydrogenase

PDHA:

Pyruvate Dehydrogenase Isoform A

PDHB:

Pyruvate Dehydrogenase Isoform B

PFK 1:

Phosphofructokinase 1

PSM:

Presomitic Mesoderm

qPCR:

Quantitative Polymerase Chain Reaction

SEM:

Standard Error of the Mean

SHH:

Sonic Hedgehog

TCA:

Tricarboxylic Acid Cycle

WNT:

Wingless-related Integration Site

References

  1. Whitsett JA, Kalin TV, Xu Y, Kalinichenko VV. Building and regenerating the lung cell by cell. Physiol Rev. 2019. https://doi.org/10.1152/physrev.00001.2018.

    Article  PubMed  Google Scholar 

  2. Maina JN. A systematic study of the development of the airway (bronchial) system of the avian lung from days 3 to 26 of embryogenesis: a transmission electron microscopic study on the domestic fowl Gallus gallus variant domesticus. Tissue Cell. 2003. https://doi.org/10.1016/s0040-8166(03)00058-2.

    Article  PubMed  Google Scholar 

  3. Kim HY, Varner VD, Nelson CM. Apical constriction initiates new bud formation during monopodial branching of the embryonic chicken lung. Development. 2013. https://doi.org/10.1242/dev.093682.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Metzger RJ, Klein OD, Martin GR, Krasnow MA. The branching programme of mouse lung development. Nature. 2008. https://doi.org/10.1038/nature07005.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Maina JN. Comparative molecular developmental aspects of the mammalian- and the avian lungs, and the insectan tracheal system by branching morphogenesis: recent advances and future directions. Front Zool. 2012. https://doi.org/10.1186/1742-9994-9-16.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Moura RS, Correia-Pinto J. Molecular aspects of avian lung development. In: Maina JN, editor. The biology of the avian respiratory system. Cham: Springer; 2017. p. 129–46.

    Chapter  Google Scholar 

  7. Moura RS, Coutinho-Borges JP, Pacheco AP, Damota PO, Correia-Pinto J. FGF signaling pathway in the developing chick lung: expression and inhibition studies. PLoS ONE. 2011. https://doi.org/10.1371/journal.pone.0017660.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Moura RS, Carvalho-Correia E, daMota P, Correia-Pinto J. Canonical Wnt signaling activity in early stages of chick lung development. PLoS ONE. 2014. https://doi.org/10.1371/journal.pone.0112388.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Moura RS, Silva-Gonçalves C, Vaz-Cunha P, Correia-Pinto J. Expression analysis of Shh signaling members in early stages of chick lung development. Histochem Cell Biol. 2016. https://doi.org/10.1007/s00418-016-1448-1.

    Article  PubMed  Google Scholar 

  10. Fernandes-Silva H, Correia-Pinto J, Moura RS. Canonical sonic hedgehog signaling in early lung development. J Dev Biol. 2017. https://doi.org/10.3390/jdb5010003.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Fernandes-Silva H, Vaz-Cunha P, Barbosa VB, Silva-Gonçalves C, Correia-Pinto J, Moura RS. Retinoic acid regulates avian lung branching through a molecular network. Cell Mol Life Sci. 2017. https://doi.org/10.1007/s00018-017-2600-3.

    Article  PubMed  Google Scholar 

  12. Warburton D, El-Hashash A, Carraro G, Tiozzo C, Sala F, Rogers O, et al. Lung organogenesis. Curr Top Dev Biol. 2010. https://doi.org/10.1016/S0070-2153(10)90003-3.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Fernandes-Silva H, Araújo-Silva H, Correia-Pinto J, Moura RS. Retinoic acid: a key regulator of lung development. Biomolecules. 2020. https://doi.org/10.3390/biom10010152.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Engle MJ, Dooley M, Brown DJ. Evidence for lactate utilization for fetal lung glycogen synthesis. Biochem Biophys Res Commun. 1987. https://doi.org/10.1016/0006-291x(87)91335-0.

    Article  PubMed  Google Scholar 

  15. Engle MJ, Brown DJ, Dehring AF, Dooley M. Effect of lactate on glucose incorporation into fetal lung phospholipids. Exp Lung Res. 1988. https://doi.org/10.3109/01902148809062854.

    Article  PubMed  Google Scholar 

  16. Yeager HJ, Massaro D. Glucose metabolism and glycoprotein synthesis by lung slices. J Appl Physiol. 1972. https://doi.org/10.1152/jappl.1972.32.4.477.

    Article  PubMed  Google Scholar 

  17. Fisher AB. Normal and pathologic biochemistry of the lung. Environ Health Perspect. 1976. https://doi.org/10.1289/ehp.76163.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Fisher AB. Intermediary metabolism of the lung. Environ Health Perspect. 1984. https://doi.org/10.1289/ehp.8455149.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Krejci A, Tennessen JM. Metabolism in time and space - exploring the frontier of developmental biology. Development. 2017. https://doi.org/10.1242/dev.150573.

    Article  PubMed  Google Scholar 

  20. Miyazawa H, Aulehla A. Revisiting the role of metabolism during development. Development. 2018. https://doi.org/10.1242/dev.131110.

    Article  PubMed  Google Scholar 

  21. Miyazawa H, Yamaguchi Y, Sugiura Y, Honda K, Kondo K, Matsuda F, et al. Rewiring of embryonic glucose metabolism via suppression of PFK-1 and aldolase during mouse chorioallantoic branching. Development. 2017. https://doi.org/10.1242/dev.138545.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Slaninova V, Krafcikova M, Perez-Gomez R, Steffal P, Trantirek L, Bray SJ, et al. Notch stimulates growth by direct regulation of genes involved in the control of glycolysis and the tricarboxylic acid cycle. Open Biol. 2016. https://doi.org/10.1098/rsob.150155.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Homem CCF, Steinmann V, Burkard TT, Jais A, Esterbauer H, Knoblich JA. Ecdysone and mediator change energy metabolism to terminate proliferation in Drosophila neural stem cells. Cell. 2014. https://doi.org/10.1016/j.cell.2014.06.024.

    Article  PubMed  Google Scholar 

  24. Bulusu V, Prior N, Snaebjornsson MT, Kuehne A, Sonnen FK, Kress J, et al. Spatiotemporal analysis of a glycolytic activity gradient linked to mouse embryo mesoderm development. Dev Cell. 2017. https://doi.org/10.1016/j.devcel.2017.01.015.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Oginuma M, Moncuquet P, Xiong F, Karoly E, Chal J, Guevorkian K, et al. A gradient of glycolytic activity coordinates FGF and Wnt signaling during elongation of the body axis in amniote embryos. Dev Cell. 2017. https://doi.org/10.1016/j.devcel.2017.02.001.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Herriges M, Morrisey EE. Lung development: orchestrating the generation and regeneration of a complex organ. Development. 2014. https://doi.org/10.1242/dev.098186.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Vander Heiden MG, DeBerardinis RJ. Understanding the intersections between metabolism and cancer biology. Cell. 2017. https://doi.org/10.1016/j.cell.2016.12.039.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Schittny JC. Development of the lung. Cell Tissue Res. 2017. https://doi.org/10.1007/s00441-016-2545-0.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Yeganeh B, Bilodeau C, Post M. Explant culture for studying lung development. In: Delgado-Olguin P, editor. Mouse embryogenesis; methods in molecular biology. New York: Humana Press; 2018. p. 81–90.

    Chapter  Google Scholar 

  30. Moura RS. Retinoic acid as a modulator of proximal-distal patterning and branching morphogenesis of the avian lung. In: Ray S, editor. Retinoid and rexinoid signaling; methods in molecular biology. New York: Humana Press; 2019. p. 209–24.

    Chapter  Google Scholar 

  31. Spurlin JW, Siedlik MJ, Nerger BA, Pang M, Jayaraman S, Zhang R, et al. Mesenchymal proteases and tissue fluidity remodel the extracellular matrix during airway epithelial branching in the embryonic avian lung. Development. 2018. https://doi.org/10.1242/dev.175257.

    Article  Google Scholar 

  32. Carraro G, del Moral P, Warburton D. Mouse embryonic lung culture, a system to evaluate the molecular mechanisms of branching. J Vis Exp. 2010. https://doi.org/10.3791/2035.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Tierney DF. Lung metabolism and biochemistry. Annu Rev Physiol. 1974. https://doi.org/10.1146/annurev.ph.36.030174.001233.

    Article  PubMed  Google Scholar 

  34. Carver FM, Shibley IA Jr, Pennington JS, Pennington SN. Differential expression of glucose transporters during chick embryogenesis. Cell Mol Life Sci. 2001. https://doi.org/10.1007/PL00000887.

    Article  PubMed  Google Scholar 

  35. Seki Y, Sato K, Kono T, Abe H, Akiba Y. Broiler chickens (Ross strain) lack insulin-responsive glucose transporter GLUT4 and have GLUT8 cDNA. Gen Comp Endocrinol. 2003. https://doi.org/10.1016/s0016-6480(03)00145-x.

    Article  PubMed  Google Scholar 

  36. Kono T, Nishida M, Nishiki Y, Seki Y, Sato K, Akiba Y. Characterisation of glucose transporter (GLUT) gene expression in broiler chickens. Br Poult Sci. 2005. https://doi.org/10.1080/00071660500181289.

    Article  PubMed  Google Scholar 

  37. Byers MS, Howard C, Wang X. Avian and mammalian facilitative glucose transporters. Microarrays. 2017. https://doi.org/10.3390/microarrays6020007.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Granja S, Morais-Santos F, Miranda-Gonçalves V, Viana-Ferreira M, Nogueira R, Nogueira-Silva C, et al. The monocarboxylate transporter inhibitor α-cyano-4-hydroxycinnamic acid disrupts rat lung branching. Cell Physiol Biochem. 2013. https://doi.org/10.1159/000356617.

    Article  PubMed  Google Scholar 

  39. Ganapathy V, Thangaraju M, Prasad PD. Nutrient transporters in cancer: relevance to Warburg hypothesis and beyond. Pharmacol Ther. 2009. https://doi.org/10.1016/j.pharmthera.2008.09.005.

    Article  PubMed  Google Scholar 

  40. Masin M, Vazquez J, Rossi S, Groeneveld S, Samson N, Schwalie PC, et al. GLUT3 is induced during epithelial-mesenchymal transition and promotes tumor cell proliferation in non-small cell lung cancer. Cancer Metab. 2014. https://doi.org/10.1186/2049-3002-2-11.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Nelson DL, Cox MM. Lehninger principles of biochemistry. 5th ed. New York: Worth Publishers; 2008.

    Google Scholar 

  42. Tessem M, Swanson MG, Keshari KR, Albers MJ, Joun D, Tabatabai ZL, et al. Evaluation of lactate and alanine as metabolic biomarkers of prostate cancer using 1H HR-MAS spectroscopy of biopsy tissues. Magn Reson Med. 2008. https://doi.org/10.1002/mrm.21694.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Sousa CM, Biancur DE, Wang X, Halbrook CJ, Sherman MH, Zhang L, et al. Pancreatic stellate cells support tumour metabolism through autophagic alanine secretion. Nature. 2016. https://doi.org/10.1038/nature19084.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009. https://doi.org/10.1126/science.1160809.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Porporato PE, Dhup S, Dadhich RK, Copetti T, Sonveaux P. Anticancer targets in the glycolytic metabolism of tumors: a comprehensive review. Front Pharmacol. 2011. https://doi.org/10.3389/fphar.2011.00049.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012. https://doi.org/10.1038/nrg3185.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Koussounadis A, Langdon SP, Um IH, Harrison DJ, Smith VA. Relationship between differentially expressed mRNA and mRNA-protein correlations in a xenograft model system. Sci Rep. 2015. https://doi.org/10.1038/srep10775.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Oginuma M, Harima Y, Tarazona OA, Diaz-Cuadros M, Michaut A, Ishitani T, et al. Intracellular pH controls WNT downstream of glycolysis in amniote embryos. Nature. 2020. https://doi.org/10.1038/s41586-020-2428-0.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Payen VL, Hsu MY, Rädecke KS, Wyart E, Vazeille T, Bouzin C, et al. Monocarboxylate transporter MCT1 promotes tumor metastasis independently of its activity as a lactate transporter. Cancer Res. 2017. https://doi.org/10.1158/0008-5472.CAN-17-0764.

    Article  PubMed  Google Scholar 

  50. Pérez-Escuredo J, Van Hée VF, Sboarina M, Falces J, Payen VL, Pellerin L, et al. Monocarboxylate transporters in the brain and in cancer. Biochim Biophys Acta. 2016. https://doi.org/10.1016/j.bbamcr.2016.03.013.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Bourgeois NMA, Van Herck SLJ, Vancamp P, Delbaere J, Zevenbergen C, Kersseboom S, et al. Characterization of chicken thyroid hormone transporters. Endocrinology. 2016. https://doi.org/10.1210/en.2015-2025.

    Article  PubMed  Google Scholar 

  52. Vancamp P, Bourgeois NMA, Houbrechts AM, Darras VM. Knockdown of the thyroid hormone transporter MCT8 in chicken retinal precursor cells hampers early retinal development and results in a shift towards more UV/blue cones at the expense of green/red cones. Exp Eye Res. 2019. https://doi.org/10.1016/j.exer.2018.09.018.

    Article  PubMed  Google Scholar 

  53. Delbaere J, Vancamp P, Van Herck SLJ, Bourgeois NMA, Green MJ, Wingate RJT, et al. MCT8 deficiency in Purkinje cells disrupts embryonic chicken cerebellar development. J Endocrinol. 2017. https://doi.org/10.1530/JOE-16-0323.

    Article  PubMed  Google Scholar 

  54. Friesema ECH, Ganguly S, Abdalla A, Fox JEM, Halestrap AP, Visser TJ. Identification of monocarboxylate transporter 8 as a specific thyroid hormone transporter. J Biol Chem. 2003. https://doi.org/10.1074/jbc.M300909200.

    Article  PubMed  Google Scholar 

  55. Liu X, Cooper DE, Cluntun AA, Warmoes MO, Zhao S, Reid MA, et al. Acetate production from glucose and coupling to mitochondrial metabolism in mammals. Cell. 2018. https://doi.org/10.1016/j.cell.2018.08.040.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Veldhuizen R, Nag K, Orgeig S, Possmayer F. The role of lipids in pulmonary surfactant. Biochim Biophys Acta. 1998. https://doi.org/10.1016/s0925-4439(98)00061-1.

    Article  PubMed  Google Scholar 

  57. Gándara L, Wappner P. Metabo-Devo: a metabolic perspective of development. Mech Dev. 2018. https://doi.org/10.1016/j.mod.2018.02.004.

    Article  PubMed  Google Scholar 

  58. Alves MG, Oliveira PF, Martins FO, Oliveira PJ, Carvalho RA. Gender-dependent metabolic remodeling during heart preservation in cardioplegic celsior and histidine buffer solution. J Cardiovasc Pharmacol. 2012. https://doi.org/10.1097/FJC.0b013e3182391d17.

    Article  PubMed  Google Scholar 

  59. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)). Method Methods. 2001. https://doi.org/10.1006/meth.2001.1262.

    Article  PubMed  Google Scholar 

  60. Kling DE, Lorenzo HK, Trbovich AM, Kinane BT, Donahoe PK, Schnitzer JJ. MEK-1/2 inhibition reduces branching morphogenesis and causes mesenchymal cell apoptosis in fetal rat lungs. Am J Physiol Lung Cell Mol Physiol. 2002. https://doi.org/10.1152/ajplung.00200.2001.

    Article  PubMed  Google Scholar 

  61. Dady A, Blavet C, Duband JP. Timing and kinetics of E- to N-cadherin switch during neurulation in the avian embryo. Dev Dyn. 2012. https://doi.org/10.1002/dvdy.23813.

    Article  PubMed  Google Scholar 

  62. Silva AM, Oliveira PJ. Evaluation of respiration with clark-type electrode in isolated mitochondria and permeabilized animal cells. In: Palmeira C, Monero A, editors. Mitochondrial bioenergetics; methods in molecular biology. New York: Humana Press; 2018. p. 7–29.

    Chapter  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Ana Lima for the histology support, Goreti Pinto for the confocal microscopy and histology support, and Maria João Meneses and Raquel Bernardino for all the help in performing the qPCR experiments.

Funding

This work has been funded by National funds, through the Foundation for Science and Technology (FCT)—project UIDB/50026/2020 and UIDP/50026/2020; by ICVS Scientific Microscopy Platform, member of the national infrastructure PPBI—Portuguese Platform of Bioimaging (PPBI-POCI-01-0145-FEDER-022122); and by the projects NORTE-01-0145-FEDER-000013 and NORTE-01-0145-FEDER-000023, supported by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF). Hugo Fernandes-Silva was supported by a doctoral fellowship (PD/BD/137655/2018) from FCT, as part of the Inter-University Doctoral Programme in Ageing and Chronic Disease (PhDOC).

Author information

Authors and Affiliations

  1. Life and Health Sciences Research Institute (ICVS), School of Medicine, University of Minho, 4710-057, Braga, Portugal

    Hugo Fernandes-Silva, Henrique Araújo-Silva, Jorge Correia-Pinto & Rute S. Moura

  2. ICVS/3B’s - PT Government Associate Laboratory, 4710-057, Braga/Guimarães, Portugal

    Hugo Fernandes-Silva, Henrique Araújo-Silva, Jorge Correia-Pinto & Rute S. Moura

  3. PhDOC PhD Program, ICVS/3B’s, School of Medicine, University of Minho, 4710-057, Braga, Portugal

    Hugo Fernandes-Silva

  4. Unit for Multidisciplinary Research in Biomedicine (UMIB), Institute of Biomedical Sciences Abel Salazar (ICBAS), University of Porto, 4050-313, Porto, Portugal

    Marco G. Alves & Ana M. Silva

  5. Department of Microscopy, Institute of Biomedical Sciences Abel Salazar (ICBAS), University of Porto, 4050-313, Porto, Portugal

    Marco G. Alves

  6. Department of Pediatric Surgery, Hospital of Braga, 4710-243, Braga, Portugal

    Jorge Correia-Pinto

  7. QOPNA &, LAQV, Department of Chemistry, University of Aveiro, 3810-193, Aveiro, Portugal

    Pedro F. Oliveira

Authors
  1. Hugo Fernandes-Silva
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Marco G. Alves
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Henrique Araújo-Silva
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ana M. Silva
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jorge Correia-Pinto
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Pedro F. Oliveira
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Rute S. Moura
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceptualization, HFS, RSM, MGA, and PFO; Methodology, RSM, MGA, and PFO; Validation, HFS, RSM, MGA, and PFO; Formal Analysis, HFS, RSM, MGA, and PFO; Investigation, HFS, HAS, AMS, and RSM; Resources, RSM, JCP, MGA, and PFO; Writing—Original Draft preparation, HFS and RSM; Writing—Review and Editing, HFS and RSM; Supervision, RSM, MGA, and PFO; Project Administration, RSM and JCP; Funding Acquisition, RSM and JCP. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Rute S. Moura.

Ethics declarations

Ethics approval and consent to participants

This work was performed at the early stages of chicken embryonic development and, therefore, does not require ethical approval following the European Parliament Directive 2010/63/EU of 22 September 2010 and the Portuguese Directive 113/2013 of 7 August 2013 on the protection of animals used for scientific purposes.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Primers and qPCR conditions. Primer sequences forward (Fw) and reverse (Rv), corresponding PCR product size, annealing temperature and number of cycles. Figure S1. Morphometric analysis of lung explants. Figure S2. Uncropped images of LDHA and total LDH immunoblot represented in Fig. 6.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fernandes-Silva, H., Alves, M.G., Araújo-Silva, H. et al. Lung branching morphogenesis is accompanied by temporal metabolic changes towards a glycolytic preference. Cell Biosci 11, 134 (2021). https://doi.org/10.1186/s13578-021-00654-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13578-021-00654-w

Keywords

Cell & Bioscience

ISSN: 2045-3701

Contact us