Programmed and self-organized flow of information during morphogenesis.

How the shape of embryos and organs emerges during development is a fundamental question that has fascinated scientists for centuries. Tissue dynamics arise from a small set of cell behaviours, including shape changes, cell contact remodelling, cell migration, cell division and cell extrusion. These behaviours require control over cell mechanics, namely active stresses associated with protrusive, contractile and adhesive forces, and hydrostatic pressure, as well as material properties of cells that dictate how cells respond to active stresses. In this Review, we address how cell mechanics and the associated cell behaviours are robustly organized in space and time during tissue morphogenesis. We first outline how not only gene expression and the resulting biochemical cues, but also mechanics and geometry act as sources of morphogenetic information to ultimately define the time and length scales of the cell behaviours driving morphogenesis. Next, we present two idealized modes of how this information flows - how it is read out and translated into a biological effect - during morphogenesis. The first, akin to a programme, follows deterministic rules and is hierarchical. The second follows the principles of self-organization, which rests on statistical rules characterizing the system's composition and configuration, local interactions and feedback. We discuss the contribution of these two modes to the mechanisms of four very general classes of tissue deformation, namely tissue folding and invagination, tissue flow and extension, tissue hollowing and, finally, tissue branching. Overall, we suggest a conceptual framework for understanding morphogenetic information that encapsulates genetics and biochemistry as well as mechanics and geometry as information modules, and the interplay of deterministic and self-organized mechanisms of their deployment, thereby diverging considerably from the traditional notion that shape is fully encoded and determined by genes.

Insights into gene expression changes under conditions that facilitate horizontal gene transfer (mating) of a model archaeon.

Horizontal gene transfer is a means by which bacteria, archaea, and eukaryotes are able to trade DNA within and between species. While there are a variety of mechanisms through which this genetic exchange can take place, one means prevalent in the archaeon Haloferax volcanii involves the transient formation of cytoplasmic bridges between cells and is referred to as mating. This process can result in the exchange of very large fragments of DNA between the participating cells. Genes governing the process of mating, including triggers to initiate mating, mechanisms of cell fusion, and DNA exchange, have yet to be characterized. We used a transcriptomic approach to gain a more detailed knowledge of how mating might transpire. By examining the differential expression of genes expressed in cells harvested from mating conditions on a filter over time and comparing them to those expressed in a shaking culture, we were able to identify genes and pathways potentially associated with mating. These analyses provide new insights into both the mechanisms and barriers of mating in Hfx. volcanii.

Tunable Gene Expression System Independent of Downstream Coding Sequence.

Fine control of the expression levels of proteins constitutes a major challenge in synthetic biology and metabolic engineering. However, the dependence of translation initiation on the downstream coding sequence (CDS) obscures accurate prediction of the protein expression levels from mRNA sequences. Here, we present a tunable gene-expression system comprising 24 expression cassettes that produce predefined relative expression levels of proteins ranging from 0.001 to 1 without being influenced by the downstream CDS. To validate the practical utility of the tunable expression system, it was applied to a synthetic circuit displaying three states of fluorescence depending on the difference in protein expression levels. To demonstrate the suitability of application to metabolic engineering, this system was used to diversify the levels of key metabolic enzymes. As a result, expression-optimized strains were capable of producing 2.25 g/L of cadaverine, 2.59 g/L of L-proline, and 95.7 mg/L of 1-propanol. Collectively, the tunable expression system could be utilized to optimize genetic circuits for desired operation and to optimize metabolic fluxes through biosynthetic pathways for enhancing production yields of bioproducts. This tunable system will be useful for studying basic and applied biological sciences in addition to applications in synthetic biology and metabolic engineering.

Expansion of Intrinsically Disordered Proteins Increases the Range of Stability of Liquid-Liquid Phase Separation.

Proteins containing intrinsically disordered regions (IDRs) are ubiquitous within biomolecular condensates, which are liquid-like compartments within cells formed through liquid-liquid phase separation (LLPS). The sequence of amino acids of a protein encodes its phase behaviour, not only by establishing the patterning and chemical nature (e.g., hydrophobic, polar, charged) of the various binding sites that facilitate multivalent interactions, but also by dictating the protein conformational dynamics. Besides behaving as random coils, IDRs can exhibit a wide-range of structural behaviours, including conformational switching, where they transition between alternate conformational ensembles. Using Molecular Dynamics simulations of a minimal coarse-grained model for IDRs, we show that the role of protein conformation has a non-trivial effect in the liquid-liquid phase behaviour of IDRs. When an IDR transitions to a conformational ensemble enriched in disordered extended states, LLPS is enhanced. In contrast, IDRs that switch to ensembles that preferentially sample more compact and structured states show inhibited LLPS. This occurs because extended and disordered protein conformations facilitate LLPS-stabilising multivalent protein-protein interactions by reducing steric hindrance; thereby, such conformations maximize the molecular connectivity of the condensed liquid network. Extended protein configurations promote phase separation regardless of whether LLPS is driven by homotypic and/or heterotypic protein-protein interactions. This study sheds light on the link between the dynamic conformational plasticity of IDRs and their liquid-liquid phase behaviour.

Species identity and diversity effects on invasion resistance of tropical freshwater plant communities.

Biotic resistance mediated by native plant diversity has long been hypothesized to reduce the success of invading plant species in terrestrial systems in temperate regions. However, still little is known about the mechanisms driving invasion patterns in other biomes or latitudes. We help to fill this gap by investigating how native plant community presence and diversity, and the presence of native phylogenetically closely related species to an invader, would affect invader Hydrilla verticillata establishment success in tropical freshwater submerged plant communities. The presence of a native community suppressed the growth of H. verticillata, but did not prevent its colonisation. Invader growth was negatively affected by native plant productivity, but independent of native species richness and phylogenetic relatedness to the invader. Native plant production was not related to native species richness in our study. We show that resistance in these tropical aquatic submerged plant communities is mainly driven by the presence and biomass of a native community independent of native species diversity. Our study illustrates that resistance provided by these tropical freshwater submerged plant communities to invasive species contrasts to resistance described for other ecosystems. This emphasizes the need to include understudied systems when predicting patterns of species invasiveness and ecosystem invasibility across biomes.

N-Terminal Region of Vibrio parahemolyticus Thermostable Direct Hemolysin Regulates the Membrane-Damaging Action of the Toxin.

Thermostable direct hemolysin (TDH) of Vibrio parahemolyticus is a membrane-damaging pore-forming toxin with potent cytolytic/cytotoxic activity. TDH exists as a tetramer consisting of protomers with a core ß-sandwich domain, flanked by an 11-amino acid long N-terminal region (NTR). This NTR could not be modeled in the previously determined crystal structure of TDH. Moreover, the functional implication of NTR for the membrane-damaging action of TDH remains unknown. In the present study, we have explored the implications of NTR for the structure-function mechanism of TDH. Our data show that the presence of NTR modulates the physicochemical property of TDH in terms of augmenting the amyloidogenic propensity of the protein. Deletion of NTR compromises the binding of TDH toward target cell membranes and drastically affects the membrane-damaging cytolytic/cytotoxic activity of the toxin. Mutations of aromatic/hydrophobic residues within NTR also confer compromised cell-killing activity. Moreover, covalent trapping of NTR, via an engineered disulfide bond, against the core ß-sandwich domain also abrogates the cytolytic/cytotoxic activity of TDH. This observation suggests that an unrestrained configuration of NTR is crucial for the membrane-damaging action of TDH. On the basis of our study, we propose a model explaining the role of NTR in the membrane-damaging function of TDH.

A microphysiological model of the bronchial airways reveals the interplay of mechanical and biochemical signals in bronchospasm.

In asthma, the contraction of the airway smooth muscle and the subsequent decrease in airflow involve a poorly understood set of mechanical and biochemical events. Organ-level and molecular-scale models of the airway are frequently based on purely mechanical or biochemical considerations and do not account for physiological mechanochemical couplings. Here, we present a microphysiological model of the airway that allows for the quantitative analysis of the interactions between mechanical and biochemical signals triggered by compressive stress on epithelial cells. We show that a mechanical stimulus mimicking a bronchospastic challenge triggers the marked contraction and delayed relaxation of airway smooth muscle, and that this is mediated by the discordant expression of cyclooxygenase genes in epithelial cells and regulated by the mechanosensor and transcriptional co-activator Yes-associated protein. A mathematical model of the intercellular feedback interactions recapitulates aspects of obstructive disease of the airways, which include pathognomonic features of severe difficult-to-treat asthma. The microphysiological model could be used to investigate the mechanisms of asthma pathogenesis and to develop therapeutic strategies that disrupt the positive feedback loop that leads to persistent airway constriction.

Enzymatic origin and various curvatures of metabolic scaling in microbes.

The famous and controversial power law is a basal metabolic scaling model mainly derived from the "surface rule" or a fractal transport network. However, this law neglects biological mechanisms in the important active state. Here, we hypothesized that the relative metabolic rate and growth rate of actively growing microbes are driven by the changeable rate of their rate-limiting enzymes and concluded that natural logarithmic microbial metabolism (lnλ) and growth (or biomass) (lnM) are both dependent on limiting resources, and then developed novel models with interdependence between lnλ and lnM. We tested the models using the data obtained from the literature. We explain how and why the scaling is usually curved with the difference between microbial metabolic and growth (or biomass's) half-saturation constants (KM, Kλ) in the active state and agree that the linear relationship of the power law is a particular case under the given condition: KM = Kλ, which means that the enzyme dynamics may drive active and basal metabolic scaling relationships. Our interdependent model is more general than the power law, which is important for integrating the ecology and biochemical processes.

A brief review of bioluminescent systems (2019).

Despite being widely used in reporter technologies, bioluminescent systems are largely understudied. Of at least forty different bioluminescent systems thought to exist in nature, molecular components of only seven light-emitting reactions are known, and the full biochemical pathway leading to light emission is only understood for two of them. Here, we provide a succinct overview of currently known bioluminescent systems highlighting available tools for research and discussing future applications.

SimCAL: a flexible tool to compute biochemical reaction similarity.

BACKGROUND: Computation of reaction similarity is a pre-requisite for several bioinformatics applications including enzyme identification for specific biochemical reactions, enzyme classification and mining for specific inhibitors. Reaction similarity is often assessed at either two levels: (i) comparison across all the constituent substrates and products of a reaction, reaction level similarity, (ii) comparison at the transformation center with various degrees of neighborhood, transformation level similarity. Existing reaction similarity computation tools are designed for specific applications and use different features and similarity measures. A single system integrating these diverse features enables comparison of the impact of different molecular properties on similarity score computation. RESULTS: To address these requirements, we present SimCAL, an integrated system to calculate reaction similarity with novel features and capability to perform comparative assessment. SimCAL provides reaction similarity computation at both whole reaction level and transformation level. Novel physicochemical features such as stereochemistry, mass, volume and charge are included in computing reaction fingerprint. Users can choose from four different fingerprint types and nine molecular similarity measures. Further, a comparative assessment of these features is also enabled. The performance of SimCAL is assessed on 3,688,122 reaction pairs with Enzyme Commission (EC) number from MetaCyc and achieved an area under the curve (AUC) of > 0.9. In addition, SimCAL results showed strong correlation with state-of-the-art EC-BLAST and molecular signature based reaction similarity methods. CONCLUSIONS: SimCAL is developed in java and is available as a standalone tool, with intuitive, user-friendly graphical interface and also as a console application. With its customizable feature selection and similarity calculations, it is expected to cater a wide audience interested in studying and analyzing biochemical reactions and metabolic networks.

Galectin-3 is essential for proper bone cell differentiation and activity, bone remodeling and biomechanical competence in mice.

OBJECTIVE: Galectin-3 is constitutively expressed in bone cells and was recently shown to modulate osteogenic transdifferentiation of vascular smooth muscle cells and atherosclerotic calcification. However, the role of galectin-3 in bone physiology is largely undefined. To address this issue, we analyzed (1) the skeletal features of 1-, 3- and 6-month-old galectin-3 null (Lgals3-/-) and wild type (WT) mice and (2) the differentiation and function of osteoblasts and osteoclasts derived from these animals. METHODS: Long bone phenotype, gene expression profile, and remodeling were investigated by micro-computed tomography, real time-PCR, static and dynamic histomorphometry, and assessment of biochemical markers of bone resorption and formation. Bone competence was also evaluated by biomechanical testing at 3 months. In vitro, the effects of galectin-3 deficiency on bone cell differentiation and function were investigated by assessing (a) gene expression of osteoblast markers, alkaline phosphatase activity, mineralization assay, and WNT/ß-catenin signaling (of which galectin-3 is a known regulator) in osteoblasts; and (b) tartrate-resistant acid phosphatase activity and bone resorption activity in osteoclasts. RESULTS: Lgals3-/- mice revealed a wide range of age-dependent alterations including lower bone formation and higher bone resorption, accelerated age-dependent trabecular bone loss (p < 0.01 vs. WT at 3 months) and reduced bone strength (p < 0.01 vs. WT at 3 months). These abnormalities were accompanied by a steady inflammatory state, as revealed by higher bone expression of the pro-inflammatory cytokines interleukin (IL)-1ß and IL-6 (p < 0.001 vs. WT at 3 months), increased content of osteal macrophages (p < 0.01 vs. WT at 3 months), and reduced expression of markers of alternative (M2) macrophage activation. Lgals3-/- osteoblasts and osteoclasts showed impaired terminal differentiation, reduced mineralization capacity (p < 0.01 vs. WT cells) and resorption activity (p < 0.01 vs. WT cells). Mechanistically, impaired differentiation and function of Lgals3-/- osteoblasts was associated with altered WNT/ß-catenin signaling (p < 0.01 vs. WT cells). CONCLUSIONS: These data provide evidence for a contribution of galectin-3 to bone cell maturation and function, bone remodeling, and biomechanical competence, thus identifying galectin-3 as a promising therapeutic target for age-related disorders of bone remodeling.

Ocean biogeochemistry modeled with emergent trait-based genomics.

Marine ecosystem models have advanced to incorporate metabolic pathways discovered with genomic sequencing, but direct comparisons between models and "omics" data are lacking. We developed a model that directly simulates metagenomes and metatranscriptomes for comparison with observations. Model microbes were randomly assigned genes for specialized functions, and communities of 68 species were simulated in the Atlantic Ocean. Unfit organisms were replaced, and the model self-organized to develop community genomes and transcriptomes. Emergent communities from simulations that were initialized with different cohorts of randomly generated microbes all produced realistic vertical and horizontal ocean nutrient, genome, and transcriptome gradients. Thus, the library of gene functions available to the community, rather than the distribution of functions among specific organisms, drove community assembly and biogeochemical gradients in the model ocean.

Molecular mechanisms of thermal resistance of the insect trypanosomatid Crithidia thermophila.

In the present work, we investigated molecular mechanisms governing thermal resistance of a monoxenous trypanosomatid Crithidia luciliae thermophila, which we reclassified as a separate species C. thermophila. We analyzed morphology, growth kinetics, and transcriptomic profiles of flagellates cultivated at low (23°C) and elevated (34°C) temperature. When maintained at high temperature, they grew significantly faster, became shorter, with genes involved in sugar metabolism and mitochondrial stress protection significantly upregulated. Comparison with another thermoresistant monoxenous trypanosomatid, Leptomonas seymouri, revealed dramatic differences in transcription profiles of the two species with only few genes showing the same expression pattern. This disparity illustrates differences in the biology of these two parasites and distinct mechanisms of their thermotolerance, a prerequisite for living in warm-blooded vertebrates.

An efficient algorithm for computing fixed length attractors based on bounded model checking in synchronous Boolean networks with biochemical applications.

Genetic regulatory networks are the key to understanding biochemical systems. One condition of the genetic regulatory network under different living environments can be modeled as a synchronous Boolean network. The attractors of these Boolean networks will help biologists to identify determinant and stable factors. Existing methods identify attractors based on a random initial state or the entire state simultaneously. They cannot identify the fixed length attractors directly. The complexity of including time increases exponentially with respect to the attractor number and length of attractors. This study used the bounded model checking to quickly locate fixed length attractors. Based on the SAT solver, we propose a new algorithm for efficiently computing the fixed length attractors, which is more suitable for large Boolean networks and numerous attractors' networks. After comparison using the tool BooleNet, empirical experiments involving biochemical systems demonstrated the feasibility and efficiency of our approach.

Modulation of microRNAs by phytochemicals in cancer: underlying mechanisms and translational significance.

MicroRNAs (miRNAs) are small, endogenous noncoding RNAs that regulate a variety of biological processes such as differentiation, development, and survival. Recent studies suggest that miRNAs are dysregulated in cancer and play critical roles in cancer initiation, progression, and chemoresistance. Therefore, exploitation of miRNAs as targets for cancer prevention and therapy could be a promising approach. Extensive evidence suggests that many naturally occurring phytochemicals regulate the expression of numerous miRNAs involved in the pathobiology of cancer. Therefore, an understanding of the regulation of miRNAs by phytochemicals in cancer, their underlying molecular mechanisms, and functional consequences on tumor pathophysiology may be useful in formulating novel strategies to combat this devastating disease. These aspects are discussed in this review paper with an objective of highlighting the significance of these observations from the translational standpoint.

The mechanisms underlying α-amanitin resistance in Drosophila melanogaster: a microarray analysis.

The rapid evolution of toxin resistance in animals has important consequences for the ecology of species and our economy. Pesticide resistance in insects has been a subject of intensive study; however, very little is known about how Drosophila species became resistant to natural toxins with ecological relevance, such as α-amanitin that is produced in deadly poisonous mushrooms. Here we performed a microarray study to elucidate the genes, chromosomal loci, molecular functions, biological processes, and cellular components that contribute to the α-amanitin resistance phenotype in Drosophila melanogaster. We suggest that toxin entry blockage through the cuticle, phase I and II detoxification, sequestration in lipid particles, and proteolytic cleavage of α-amanitin contribute in concert to this quantitative trait. We speculate that the resistance to mushroom toxins in D. melanogaster and perhaps in mycophagous Drosophila species has evolved as cross-resistance to pesticides, other xenobiotic substances, or environmental stress factors.

Unveiling microbial life in new deep-sea hypersaline Lake Thetis. Part I: Prokaryotes and environmental settings.

In September 2008, an expedition of the RV Urania was devoted to exploration of the genomic richness of deep hypersaline anoxic lakes (DHALs) located in the Western part of the Mediterranean Ridge. Approximately 40 nautical miles SE from Urania Lake, the presence of anoxic hypersaline lake, which we named Thetis, was confirmed by swath bathymetry profiling and through immediate sampling casts. The brine surface of the Thetis Lake is located at a depth of 3258 m with a thickness of ≈ 157 m. Brine composition was found to be thalassohaline, saturated by NaCl with a total salinity of 348‰, which is one of highest value reported for DHALs. Similarly to other Mediterranean DHALs, seawater-brine interface of Thetis represents a steep pycno- and chemocline with gradients of salinity, electron donors and acceptors and posseses a remarkable stratification of prokaryotic communities, observed to be more metabolically active in the upper interface where redox gradient was sharper. [(14) C]-bicarbonate fixation analysis revealed that microbial communities are sustained by sulfur-oxidizing chemolithoautotrophic primary producers that thrive within upper interface. Besides microaerophilic autotrophy, heterotrophic sulfate reduction, methanogenesis and anaerobic methane oxidation are likely the predominant processes driving the ecosystem of Thetis Lake.

Hyperglycemia perturbs biochemical networks in human trophoblast BeWo cells.

Determining the effects of hyperglycemia on gene expression in placental trophoblast is important to gain a better understanding of how diabetes adversely affects pregnancy. In this study, we examined whether exposure to high glucose during forskolin-induced differentiation affects gene expression in differentiated trophoblasts. Human trophoblast BeWo cells were differentiated under low glucose (LG: 11 mM) or high glucose (HG: 25 mM) conditions. Gene expression was analyzed using a GeneChip system and the obtained data were analyzed using Ingenuity Pathways Analysis. In HG conditions, there were marked alterations in gene expression in differentiated BeWo cells compared with LG conditions. In particular, BeWo cells responded to HG with major changes in the expression levels of cell cycle- and metabolism-related genes. We selected the aromatase gene for further investigation of the molecular mechanisms. Mannitol or 3-O-methylglucose did not mimic the expression changes caused by HG, indicating that the effect of glucose was not due to a difference in osmotic pressure, and that glucose metabolism plays an essential role in inducing the HG effects. Cotreatment with N-acetylcysteine reduced the effect of HG on aromatase gene expression, suggesting that hyperglycemia may perturb biochemical networks because of the elevation of oxidative stress. Overall, our results will aid further understanding of the effect of diabetes on the regulation of trophoblast differentiation and function.

Quantitative global studies of reactomes and metabolomes using a vectorial representation of reactions and chemical compounds.

BACKGROUND: Global studies of the protein repertories of organisms are providing important information on the characteristics of the protein space. Many of these studies entail classification of the protein repertory on the basis of structure and/or sequence similarities. The situation is different for metabolism. Because there is no good way of measuring similarities between chemical reactions, there is a barrier to the development of global classifications of "metabolic space" and subsequent studies comparable to those done for protein sequences and structures. RESULTS: In this work, we propose a vectorial representation of chemical reactions, which allows them to be compared and classified. In this representation, chemical compounds, reactions and pathways may be represented in the same vectorial space. We show that the representation of chemical compounds reflects their physicochemical properties and can be used for predictive purposes. We use the vectorial representations of reactions to perform a global classification of the reactome of the model organism E. coli. CONCLUSIONS: We show that this unsupervised clustering results in groups of enzymes more coherent in biological terms than equivalent groupings obtained from the EC hierarchy. This hierarchical clustering produces an optimal set of 21 groups which we analyzed for their biological meaning.