Multiphoton excited polymerized biomimetic models of collagen fiber morphology to study single cell and collective migration dynamics in pancreatic cancer.

The respective roles of aligned collagen fiber morphology found in the extracellular matrix (ECM) of pancreatic cancer patients and cellular migration dynamics have been gaining attention because of their connection with increased aggressive phenotypes and poor prognosis. To better understand how collagen fiber morphology influences cell-matrix interactions associated with metastasis, we used Second Harmonic Generation (SHG) images from patient biopsies with Pancreatic ductal adenocarcinoma (PDAC) as models to fabricate collagen scaffolds to investigate processes associated with motility. Using the PDAC BxPC-3 metastatic cell line, we investigated single and collective cell dynamics on scaffolds of varying collagen alignment. Collective or clustered cells grown on the scaffolds with the highest collagen fiber alignment had increased E-cadherin expression and larger focal adhesion sites compared to single cells, consistent with metastatic behavior. Analysis of single cell motility revealed that the dynamics were characterized by random walk on all substrates. However, examining collective motility over different time points showed that the migration was super-diffusive and enhanced on highly aligned fibers, whereas it was hindered and sub-diffusive on un-patterned substrates. This was further supported by the more elongated morphology observed in collectively migrating cells on aligned collagen fibers. Overall, this approach allows the decoupling of single and collective cell behavior as a function of collagen alignment and shows the relative importance of collective cell behavior as well as fiber morphology in PDAC metastasis. We suggest these scaffolds can be used for further investigations of PDAC cell biology. STATEMENT OF SIGNIFICANCE: Pancreatic ductal adenocarcinoma (PDAC) has a high mortality rate, where aligned collagen has been associated with poor prognosis. Biomimetic models representing this architecture are needed to understand complex cellular interactions. The SHG image-based models based on stromal collagen from human biopsies afford the measurements of cell morphology, cadherin and focal adhesion expression as well as detailed motility dynamics. Using a metastatic cell line, we decoupled the roles of single cell and collective cell behavior as well as that arising from aligned collagen. Our data suggests that metastatic characteristics are enhanced by increased collagen alignment and that collective cell behavior is more relevant to metastatic processes. These scaffolds provide new insight in this disease and can be a platform for further experiments such as testing drug efficacy.

Cellular Energy Cycle Mediates an Advection-Like Forward Cell Flow to Support Collective Invasion.

Collective cell migration is a model for nonequilibrium biological dynamics, which is important for morphogenesis, pattern formation, and cancer metastasis. The current understanding of cellular collective dynamics is based primarily on cells moving within a 2D epithelial monolayer. However, solid tumors often invade surrounding tissues in the form of a stream-like 3D structure, and how biophysical cues are integrated at the cellular level to give rise to this collective streaming remains unclear. Here, it is shown that cell cycle-mediated bioenergetics drive a forward advective flow of cells and energy to the front to support 3D collective invasion. The cell division cycle mediates a corresponding energy cycle such that cellular adenosine triphosphate (ATP) energy peaks just before division. A reaction-advection-diffusion (RAD) type model coupled with experimental measurements further indicates that most cells enter an active division cycle at rear positions during 3D streaming. Once the cells progress to a later stage toward division, the high intracellular energy allows them to preferentially stream toward the tip and become leader cells. This energy-driven cellular flow may be a fundamental characteristic of 3D collective dynamics based on thermodynamic principles important for not only cancer invasion but also tissue morphogenesis.

The highly metastatic 4T1 breast carcinoma model possesses features of a hybrid epithelial/mesenchymal phenotype.

Epithelial-mesenchymal transitions (EMTs) are thought to promote metastasis via downregulation of E-cadherin (also known as Cdh1) and upregulation of mesenchymal markers such as N-cadherin (Cdh2) and vimentin (Vim). Contrary to this, E-cadherin is retained in many invasive carcinomas and promotes collective cell invasion. To investigate how E-cadherin regulates metastasis, we examined the highly metastatic, E-cadherin-positive murine 4T1 breast cancer model, together with the less metastatic, 4T1-related cell lines 4T07, 168FARN and 67NR. We found that 4T1 cells display a hybrid epithelial/mesenchymal phenotype with co-expression of epithelial and mesenchymal markers, whereas 4T07, 168FARN, and 67NR cells display progressively more mesenchymal phenotypes in vitro that relate inversely to their metastatic capacity in vivo. Using RNA interference and constitutive expression, we demonstrate that the expression level of E-cadherin does not determine 4T1 or 4T07 cell metastatic capacity in mice. Mechanistically, 4T1 cells possess highly dynamic, unstable cell-cell junctions and can undergo collective invasion without E-cadherin downregulation. However, 4T1 orthotopic tumors in vivo also contain subregions of EMT-like loss of E-cadherin. Thus, 4T1 cells function as a model for carcinomas with a hybrid epithelial/mesenchymal phenotype that promotes invasion and metastasis.

Self-extinguishing relay waves enable homeostatic control of human neutrophil swarming.

Neutrophils collectively migrate to sites of injury and infection. How these swarms are coordinated to ensure the proper level of recruitment is unknown. Using an ex vivo model of infection, we show that human neutrophil swarming is organized by multiple pulsatile chemoattractant waves. These waves propagate through active relay in which stimulated neutrophils trigger their neighbors to release additional swarming cues. Unlike canonical active relays, we find these waves to be self-terminating, limiting the spatial range of cell recruitment. We identify an NADPH-oxidase-based negative feedback loop that is needed for this self-terminating behavior. We observe near-constant levels of neutrophil recruitment over a wide range of starting conditions, revealing surprising robustness in the swarming process. This homeostatic control is achieved by larger and more numerous swarming waves at lower cell densities. We link defective wave termination to a broken recruitment homeostat in the context of human chronic granulomatous disease.

Optimal Control of Collective Electrotaxis in Epithelial Monolayers.

Epithelial monolayers are some of the best-studied models for collective cell migration due to their abundance in multicellular systems and their tractability. Experimentally, the collective migration of epithelial monolayers can be robustly steered e.g. using electric fields, via a process termed electrotaxis. Theoretically, however, the question of how to design an electric field to achieve a desired spatiotemporal movement pattern is underexplored. In this work, we construct and calibrate an ordinary differential equation model to predict the average velocity of the centre of mass of a cellular monolayer in response to stimulation with an electric field. We use this model, in conjunction with optimal control theory, to derive physically realistic optimal electric field designs to achieve a variety of aims, including maximising the total distance travelled by the monolayer, maximising the monolayer velocity, and keeping the monolayer velocity constant during stimulation. Together, this work is the first to present a unified framework for optimal control of collective monolayer electrotaxis and provides a blueprint to optimally steer collective migration using other external cues.

New Opportunities for Electric Fields in Promoting Wound Healing: Collective Electrotaxis.

Significance: It has long been hypothesized that naturally occurring electric fields (EFs) aid wound healing by guiding cell migration. Consequently, the application of EFs has significant potential for promoting wound healing. However, the mechanisms underlying the cellular response to EFs remain unclear. Recent Advances: Although the directed migration of isolated single cells under EFs has been studied for decades, only recently has experimental evidence demonstrated the distinct collective migration of large sheets of keratinocytes and corneal epithelial cells in response to applied EFs. Accumulating evidence suggests that the emergent properties of cell groups in response to EF guidance offer new opportunities for EF-assisted directional migration. Critical Issues: In this review, we provide an overview of the field of collective electrotaxis, highlighting key advances made in recent years. We also discuss advanced engineering strategies utilized to manipulate collective electrotaxis. Future Directions: We outline a series of unanswered questions in this field and propose potential applications of collective electrotaxis in developing electrical stimulation technologies for wound healing.

Experience reduces route selection for conspecifics by the collectively migrating white stork.

Migration can be an energetically costly behavior with strong fitness consequences in terms of mortality and reproduction.1,2,3,4,5,6,7,8,9,10,11 Migrants should select migratory routes to minimize their costs, but both costs and benefits may change with experience.12,13,14 This raises the question of whether experience changes how individuals select their migratory routes. Here, we investigate the effect of age on route selection criteria in a collectively migrating soaring bird, the white stork (Ciconia ciconia). We perform step-selection analysis on a longitudinal dataset tracking 158 white storks over up to 9 years to quantify how they select their routes based on the social and atmospheric environments and to examine how this selection changes with age. We find clear ontogenetic shifts in route selection criteria. Juveniles choose routes that have good atmospheric conditions and high conspecific densities. Yet, as they gain experience, storks' selection on the availability of social information reduces-after their fifth migration, experienced birds also choose routes with low conspecific densities. Thus, our results suggest that as individuals age, they gradually replace information gleaned from other individuals with information gained from experience, allowing them to shift their migration timing and increasing the timescale at which they select their routes.

Emergence, Pattern, and Frequency of Spontaneous Waves in Spreading Epithelial Monolayers.

A striking phenomenon of collective cell motion is that they can exhibit a spontaneously emerging wave during epithelia expansions. However, the fundamental mechanism, governing the emergence and its crucial characteristics (e.g., the eigenfrequency and the pattern), remains an enigma. By introducing a mechanochemical feedback loop, we develop a highly efficient discrete vertex model to investigate the spatiotemporal evolution of spreading epithelia. We find both numerically and analytically that expanding cell monolayers display a power-law dependence of wave frequency on the local heterogeneities (i.e., cell density) with a scaling exponent of -1/2. Moreover, our study demonstrates the quantitative capability of the proposed model in capturing distinct X-, W-, and V-mode wave patterns. We unveil that the phase transition between these modes is governed by the distribution of active self-propulsion forces. Our work provides an avenue for rigorous quantitative investigations into the collective motion and pattern formation of cell groups.

A simple agent-based hybrid model to simulate the biophysics of glioblastoma multiforme cells and the concomitant evolution of the oxygen field.

BACKGROUND AND OBJECTIVES: Glioblastoma multiforme (GBM) is one of the most aggressive cancers of the central nervous system. It is characterized by a high mitotic activity and an infiltrative ability of the glioma cells, neovascularization and necrosis. GBM evolution entails the continuous interplay between heterogeneous cell populations, chemotaxis, and physical cues through different scales. In this work, an agent-based hybrid model is proposed to simulate the coupling of the multiscale biological events involved in the GBM invasion, specifically the individual and collective migration of GBM cells and the concurrent evolution of the oxygen field and phenotypic plasticity. An asset of the formulation is that it is conceptually and computationally simple but allows to reproduce the complexity and the progression of the GBM micro-environment at cell and tissue scales simultaneously. METHODS: The migration is reproduced as the result of the interaction between every single cell and its micro-environment. The behavior of each individual cell is formulated through genotypic variables whereas the cell micro-environment is modeled in terms of the oxygen concentration and the cell density surrounding each cell. The collective behavior is formulated at a cellular scale through a flocking model. The phenotypic plasticity of the cells is induced by the micro-environment conditions, considering five phenotypes. RESULTS: The model has been contrasted by benchmark problems and experimental tests showing the ability to reproduce different scenarios of glioma cell migration. In all cases, the individual and collective cell migration and the coupled evolution of both the oxygen field and phenotypic plasticity have been properly simulated. This simple formulation allows to mimic the formation of relevant hallmarks of glioblastoma multiforme, such as the necrotic cores, and to reproduce experimental evidences related to the mitotic activity in pseudopalisades. CONCLUSIONS: In the collective migration, the survival of the clusters prevails at the expense of cell mitosis, regardless of the size of the groups, which delays the formation of necrotic foci and reduces the rate of oxygen consumption.

E-cadherin biomaterials reprogram collective cell migration and cell cycling by forcing homeostatic conditions.

Cells attach to the world through either cell-extracellular matrix adhesion or cell-cell adhesion, and traditional biomaterials imitate the matrix for integrin-based adhesion. However, materials incorporating cadherin proteins that mimic cell-cell adhesion offer an alternative to program cell behavior and integrate into living tissues. We investigated how cadherin substrates affect collective cell migration and cell cycling in epithelia. Our approach involved biomaterials with matrix proteins on one-half and E-cadherin proteins on the other, forming a "Janus" interface across which we grew a single sheet of cells. Tissue regions over the matrix side exhibited normal collective dynamics, but an abrupt behavior shift occurred across the Janus boundary onto the E-cadherin side, where cells attached to the substrate via E-cadherin adhesions, resulting in stalled migration and slowing of the cell cycle. E-cadherin surfaces disrupted long-range mechanical coordination and nearly doubled the length of the G0/G1 phase of the cell cycle, linked to the lack of integrin focal adhesions on the E-cadherin surface.

A biased random walk approach for modeling the collective chemotaxis of neural crest cells.

Collective cell migration is a multicellular phenomenon that arises in various biological contexts, including cancer and embryo development. 'Collectiveness' can be promoted by cell-cell interactions such as co-attraction and contact inhibition of locomotion. These mechanisms act on cell polarity, pivotal for directed cell motility, through influencing the intracellular dynamics of small GTPases such as Rac1. To model these dynamics we introduce a biased random walk model, where the bias depends on the internal state of Rac1, and the Rac1 state is influenced by cell-cell interactions and chemoattractive cues. In an extensive simulation study we demonstrate and explain the scope and applicability of the introduced model in various scenarios. The use of a biased random walk model allows for the derivation of a corresponding partial differential equation for the cell density while still maintaining a certain level of intracellular detail from the individual based setting.

N-cadherin dynamically regulates pediatric glioma cell migration in complex environments.

Pediatric high-grade gliomas are highly invasive and essentially incurable. Glioma cells migrate between neurons and glia, along axon tracts, and through extracellular matrix surrounding blood vessels and underlying the pia. Mechanisms that allow adaptation to such complex environments are poorly understood. N-cadherin is highly expressed in pediatric gliomas and associated with shorter survival. We found that inter-cellular homotypic N-cadherin interactions differentially regulate glioma migration according to the microenvironment, stimulating migration on cultured neurons or astrocytes but inhibiting invasion into reconstituted or astrocyte-deposited extracellular matrix. N-cadherin localizes to filamentous connections between migrating leader cells but to epithelial-like junctions between followers. Leader cells have more surface and recycling N-cadherin, increased YAP1/TAZ signaling, and increased proliferation relative to followers. YAP1/TAZ signaling is dynamically regulated as leaders and followers change position, leading to altered N-cadherin levels and organization. Together, the results suggest that pediatric glioma cells adapt to different microenvironments by regulating N-cadherin dynamics and cell-cell contacts.

The myocardium utilizes a platelet-derived growth factor receptor alpha (Pdgfra)-phosphoinositide 3-kinase (PI3K) signaling cascade to steer toward the midline during zebrafish heart tube formation.

Coordinated cell movement is a fundamental process in organ formation. During heart development, bilateral myocardial precursors collectively move toward the midline (cardiac fusion) to form the primitive heart tube. Extrinsic influences such as the adjacent anterior endoderm are known to be required for cardiac fusion. We previously showed however, that the platelet-derived growth factor receptor alpha (Pdgfra) is also required for cardiac fusion (Bloomekatz et al., 2017). Nevertheless, an intrinsic mechanism that regulates myocardial movement has not been elucidated. Here, we show that the phosphoinositide 3-kinase (PI3K) intracellular signaling pathway has an essential intrinsic role in the myocardium directing movement toward the midline. In vivo imaging further reveals midline-oriented dynamic myocardial membrane protrusions that become unpolarized in PI3K-inhibited zebrafish embryos where myocardial movements are misdirected and slower. Moreover, we find that PI3K activity is dependent on and interacts with Pdgfra to regulate myocardial movement. Together our findings reveal an intrinsic myocardial steering mechanism that responds to extrinsic cues during the initiation of cardiac development.

A toolbox to analyze collective cell migration, proliferation and cellular organization simultaneously.

BACKGROUND: Analyses of collective cell migration and orientation phenomena are needed to assess the behavior of multicellular clusters. While some tools to the authors' knowledge none is capable to analyze collective migration, cellular orientation and proliferation in phase contrast images simultaneously. METHODS: We provide a tool based to analyze phase contrast images of dense cell layers. PIV is used to calculatevelocity fields, while the structure tensor provides cellular orientation. An artificial neural network is used to identify cell division events, allowing to correlate migratory and organizational phenomena with cell density. CONCLUSION: The presented tool allows the simultaneous analysis of collective cell behavior from phase contrast images in terms of migration, (self-)organization and proliferation.

Phenotypically sorted highly and weakly migratory triple negative breast cancer cells exhibit migratory and metastatic commensalism.

BACKGROUND: Intratumor heterogeneity is a well-established hallmark of cancer that impedes cancer research, diagnosis, and treatment. Previously, we phenotypically sorted human breast cancer cells based on migratory potential. When injected into mice, highly migratory cells were weakly metastatic and weakly migratory cells were highly metastatic. The purpose of this study was to determine whether these weakly and highly migratory cells interact with each other in vitro or in vivo. METHODS: To assess the relationship between heterogeneity in cancer cell migration and metastatic fitness, MDA-MB-231 and SUM159PT triple negative breast cancer cells were phenotypically sorted into highly migratory and weakly migratory subpopulations and assayed separately and in a 1:1 mixture in vitro and in vivo for metastatic behaviors. Unpaired, two-tailed Student's t-tests, Mann-Whitney tests, ordinary, one-way ANOVAs, and Kruskal-Wallis H tests were performed as appropriate with p < 0.05 as the cutoff for statistical significance. RESULTS: When highly and weakly migratory cells are co-seeded in mixed spheroids, the weakly migratory cells migrated farther than weakly migratory only spheroids. In mixed spheroids, leader-follower behavior occurred with highly migratory cells leading the weakly migratory cells in migration strands. When cell suspensions of highly migratory, weakly migratory, or a 1:1 mixture of both subpopulations were injected orthotopically into mice, both the mixed cell suspensions and weakly migratory cells showed significant distal metastasis, but the highly migratory cells did not metastasize significantly to any location. Notably, significantly more distal metastasis was observed in mice injected with the 1:1 mixture compared to either subpopulation alone. CONCLUSIONS: This study suggests that weakly migratory cells interact with highly migratory cells in a commensal fashion resulting in increased migration and metastasis. Together, these findings indicate that cancer cell subpopulation migration ability does not correlate with metastatic potential and that cooperation between highly migratory and weakly migratory subpopulations can enhance overall metastatic fitness.

E-cadherin biointerfaces reprogram collective cell migration and cell cycling by forcing homeostatic conditions.

Cells attach to the world around them in two ways-cell:extracellular-matrix adhesion and cell:cell adhesion-and conventional biomaterials are made to resemble the matrix to encourage integrin-based cell adhesion. However, interest is growing for cell-mimetic interfaces that mimic cell-cell interactions using cadherin proteins, as this offers a new way to program cell behavior and design synthetic implants and objects that can integrate directly into living tissues. Here, we explore how these cadherin-based materials affect collective cell behaviors, focusing specifically on collective migration and cell cycle regulation in cm-scale epithelia. We built culture substrates where half of the culture area was functionalized with matrix proteins and the contiguous half was functionalized with E-cadherin proteins, and we grew large epithelia across this 'Janus' interface. Parts of the tissues in contact with the matrix side of the Janus interface exhibited normal collective dynamics, but an abrupt shift in behaviors happened immediately across the Janus boundary onto the E-cadherin side, where cells formed hybrid E-cadherin junctions with the substrate, migration effectively froze in place, and cell-cycling significantly decreased. E-cadherin materials suppressed long-range mechanical correlations in the tissue and mechanical information reflected off the substrate interface. These effects could not be explained by conventional density, shape index, or contact inhibition explanations. E-cadherin surfaces nearly doubled the length of the G0/G1 phase of the cell cycle, which we ultimately connected to the exclusion of matrix focal adhesions induced by the E-cadherin culture surface.

CD9 negatively regulates collective electrotaxis of the epidermal monolayer by controlling and coordinating the polarization of leader cells.

Background: Endogenous electric fields (EFs) play an essential role in guiding the coordinated collective migration of epidermal cells to the wound centre during wound healing. Although polarization of leadercells is essential for collective migration, the signal mechanisms responsible for the EF-induced polarization of leader cells under electrotactic collective migration remain unclear. This study aims to determine how the leader cells are polarized and coordinated during EF-guided collective migration of epidermal cell sheets. Methods: Collective migration of the human epidermal monolayer (human immortalized keratinocytes HaCaT) under EFs was observed via time-lapse microscopy. The involvement of tetraspanin-29 (CD9) in EF-induced fibrous actin (F-actin) polarization of leader cells as well as electrotactic migration of the epidermal monolayer was evaluated by genetic manipulation. Blocking, rescue and co-culture experiments were conducted to explore the downstream signalling of CD9. Results: EFs guided the coordinated collective migration of the epithelial monolayer to the anode, with dynamic formation of pseudopodia in leader cells at the front edge of the monolayer along the direction of migration. F-actin polarization, as expected, played an essential role in pseudopod formation in leader cells under EFs. By confocal microscopy, we found that CD9 was colocalized with F-actin on the cell surface and was particularly downregulated in leader cells by EFs. Interestingly, genetic overexpression of CD9 abolished EF-induced F-actin polarization in leader cells as well as collective migration in the epidermal monolayer. Mechanistically, CD9 determined the polarization of F-actin in leader cells by downregulating a disintegrin and metalloprotease 17/heparin-binding epidermal growth factor-like growth factor/epidermal growth factor receptor (ADAM17/HB-EGF/EGFR) signalling. The abolished polarization of leader cells due to CD9 overexpression could be restored in a co-culture monolayer where normal cells and CD9-overexpressing cells were mixed; however, this restoration was eliminated again by the addition of the HB-EGF-neutralizing antibody. Conclusion: CD9 functions as a key regulator in the EF-guided collective migration of the epidermal monolayer by controlling and coordinating the polarization of leader cells through ADAM17/HB-EGF/EGFR signalling.

Migration and division in cell monolayers on substrates with topological defects.

Collective movement and organization of cell monolayers are important for wound healing and tissue development. Recent experiments highlighted the importance of liquid crystal order within these layers, suggesting that +1 topological defects have a role in organizing tissue morphogenesis. We study fibroblast organization, motion, and proliferation on a substrate with micron-sized ridges that induce +1 and -1 topological defects using simulation and experiment. We model cells as self-propelled deformable ellipses that interact via a Gay-Berne potential. Unlike earlier work on other cell types, we see that density variation near defects is not explained by collective migration. We propose instead that fibroblasts have different division rates depending on their area and aspect ratio. This model captures key features of our previous experiments: the alignment quality worsens at high cell density and, at the center of the +1 defects, cells can adopt either highly anisotropic or primarily isotropic morphologies. Experiments performed with different ridge heights confirm a prediction of this model: Suppressing migration across ridges promotes higher cell density at the +1 defect. Our work enables a mechanism for tissue patterning using topological defects without relying on cell migration.

Self-extinguishing relay waves enable homeostatic control of human neutrophil swarming.

Neutrophils exhibit self-amplified swarming to sites of injury and infection. How swarming is controlled to ensure the proper level of neutrophil recruitment is unknown. Using an ex vivo model of infection, we find that human neutrophils use active relay to generate multiple pulsatile waves of swarming signals. Unlike classic active relay systems such as action potentials, neutrophil swarming relay waves are self-extinguishing, limiting the spatial range of cell recruitment. We identify an NADPH-oxidase-based negative feedback loop that is needed for this self-extinguishing behavior. Through this circuit, neutrophils adjust the number and size of swarming waves for homeostatic levels of cell recruitment over a wide range of initial cell densities. We link a broken homeostat to neutrophil over-recruitment in the context of human chronic granulomatous disease.

Electron irradiation effects in transmission electron microscopy: Random displacements and collective migrations.

Electron beam damage in transmission electron microscopy (TEM) is complicated because the damage phenomena can be the result of random atomic displacements or collective migrations. The former is categorized as the primary beam effects and the latter is the secondary beam effects. The mechanisms for these two distinguishing atomic processes of damage are different. The primary beam effects can be caused by the mechanisms of knock-on and/or radiolysis, while the secondary effects must be driven by a field that is induced by electron irradiation. One such field has been identified to be the electric field produced by the accumulated charges due to the ejection of secondary and Auger electrons from the irradiated region. One convincing example is the electron irradiation-induced domain switch in ferroelectric materials, in which the collective cation displacements are driven by the induced electric field. A detailed interpretation is given in this review. The sintering of metal NPs under electron irradiation is a secondary beam effect and is most likely also caused by the induced electric fields. The interactions between the charged NP and substrate, and between charged NPs, result in NP motion. Interchanging atoms between NPs during the sintering may also be driven by the electric fields. Although many beam-damage phenomena in C nanotubes and layered materials, such as graphene, BN, and transition metal dichalcogenides, are caused by the primary beam effects and have been well studied experimentally and theoretically in the literature, some phenomena from the secondary beam effects have also been identified in this review. These phenomena are sensitive to electron current density, the shape and orientation of the specimen, and even the illumination mode (i.e., TEM or STEM). Unfortunately, the mechanisms responsible for these phenomena still need to be clarified.