Experimental Perturbations Research Articles

The molecular dimension of microbial species: 1. Ecological distinctions among, and homogeneity within, putative ecotypes of Synechococcus inhabiting the cyanobacterial mat of Mushroom Spring, Yellowstone National Park.

Based on the Stable Ecotype Model, evolution leads to the divergence of ecologically distinct populations (e.g., with different niches and/or behaviors) of ecologically interchangeable membership. In this study, pyrosequencing was used to provide deep sequence coverage of Synechococcus psaA genes and transcripts over a large number of habitat types in the Mushroom Spring microbial mat. Putative ecological species [putative ecotypes (PEs)], which were predicted by an evolutionary simulation based on the Stable Ecotype Model (Ecotype Simulation), exhibited distinct distributions relative to temperature-defined positions in the effluent channel and vertical position in the upper 1 mm-thick mat layer. Importantly, in most cases variants predicted to belong to the same PE formed unique clusters relative to temperature and depth in the mat in canonical correspondence analysis, supporting the hypothesis that while the PEs are ecologically distinct, the members of each ecotype are ecologically homogeneous. PEs responded differently to experimental perturbations of temperature and light, but the genetic variation within each PE was maintained as the relative abundances of PEs changed, further indicating that each population responded as a set of ecologically interchangeable individuals. Compared to PEs that predominate deeper within the mat photic zone, the timing of transcript abundances for selected genes differed for PEs that predominate in microenvironments closer to upper surface of the mat with spatiotemporal differences in light and O2 concentration. All of these findings are consistent with the hypotheses that Synechococcus species in hot spring mats are sets of ecologically interchangeable individuals that are differently adapted, that these adaptations control their distributions, and that the resulting distributions constrain the activities of the species in space and time.

Confrontation of fibroblasts with cancer cells in vitro: gene network analysis of transcriptome changes and differential capacity to inhibit tumor growth.

BackgroundThere is growing evidence that emerging malignancies in solid tissues might be kept under control by physical intercellular contacts with normal fibroblasts.MethodsHere we characterize transcriptional landscapes of fibroblasts that confronted cancer cells. We studied four pairs of in vitro and ex vivo fibroblast lines which, within each pair, differed in their capacity to inhibit cancer cells. The natural process was modeled in vitro by confronting the fibroblasts with PC-3 cancer cells. Fibroblast transcriptomes were recorded by Affymetrix microarrays and then investigated using network analysis.ResultsThe network enrichment analysis allowed us to separate confrontation- and inhibition-specific components of the fibroblast transcriptional response. Confrontation-specific differences were stronger and were characterized by changes in a number of pathways, including Rho, the YAP/TAZ cascade, NF-kB, and TGF-beta signaling, as well as the transcription factor RELA. Inhibition-specific differences were more subtle and characterized by involvement of Rho signaling at the pathway level and by potential individual regulators such as IL6, MAPK8, MAP2K4, PRKCA, JUN, STAT3, and STAT5A.ConclusionsWe investigated the interaction between cancer cells and fibroblasts in order to shed light on the potential mechanisms and explain the differential inhibitory capacity of the latter, which enabled both a holistic view on the process and details at the gene/protein level. The combination of our methods pointed to proteins, such as members of the Rho pathway, pro-inflammatory signature and the YAP1/TAZ cascade, that warrant further investigation via tools of experimental perturbation. We also demonstrated functional congruence between the in vitro and ex vivo models.The microarray data are made available via the Gene Expression Omnibus as GSE57199.Electronic supplementary materialThe online version of this article (doi:10.1186/s13046-015-0178-x) contains supplementary material, which is available to authorized users.

GABAB receptor cell-surface export is controlled by an endoplasmic reticulum gatekeeper.

Endoplasmic reticulum (ER) release and cell-surface export of many G protein-coupled receptors (GPCRs) are tightly regulated. For gamma-aminobutyric acid (GABA)B receptors of GABA, the major mammalian inhibitory neurotransmitter, the ligand-binding GB1 subunit is maintained in the ER by unknown mechanisms in the absence of hetero-dimerization with the GB2 subunit. We report that GB1 retention is regulated by a specific gatekeeper, PRAF2. This ER resident transmembrane protein binds to GB1, preventing its progression in the biosynthetic pathway. GB1 release occurs upon competitive displacement from PRAF2 by GB2. PRAF2 concentration, relative to that of GB1 and GB2, tightly controls cell-surface receptor density and controls GABAB function in neurons. Experimental perturbation of PRAF2 levels in vivo caused marked hyperactivity disorders in mice. These data reveal an unanticipated major impact of specific ER gatekeepers on GPCR function and identify PRAF2 as a new molecular target with therapeutic potential for psychiatric and neurological diseases involving GABAB function.

Dynamic phosphoproteomics reveals TORC1-dependent regulation of yeast nucleotide and amino acid biosynthesis.

Phosphoproteomics studies have unraveled the extent of protein phosphorylation as a key cellular regulation mechanism, but assigning functionality to specific phosphorylation events remains a major challenge. TORC1 (target of rapamycin complex 1) is a kinase-containing protein complex that transduces changes in nutrient availability into phosphorylation signaling events that alter cell growth and proliferation. To resolve the temporal sequence of phosphorylation responses to nutritionally and chemically induced changes in TORC1 signaling and to identify previously unknown kinase-substrate relationships in Saccharomyces cerevisiae, we performed quantitative mass spectrometry-based phosphoproteomic analyses after shifts in nitrogen sources and rapamycin treatment. From early phosphorylation events that were consistent over at least two experimental perturbations, we identified 51 candidate and 10 known proximal targets of TORC1 that were direct substrates of TORC1 or of one of its kinase or phosphatase substrates. By correlating these phosphoproteomics data with dynamic metabolomics data, we inferred the functional role of phosphorylation on the metabolic activity of 12 enzymes, including three candidate TORC1-proximal targets: Amd1, which is involved in nucleotide metabolism; Hom3, which is involved in amino acid metabolism; and Tsl1, which mediates carbohydrate storage. Finally, we identified the TORC1 substrates Sch9 and Atg1 as candidate kinases that phosphorylate Amd1 and Hom3, respectively.

Establishing biomechanical mechanisms in mouse models: practical guidelines for systematically evaluating phenotypic changes in the diaphyses of long bones.

Mice are widely used in studies of skeletal biology, and assessment of their bones by mechanical testing is a critical step when evaluating the functional effects of an experimental perturbation. For example, a gene knockout may target a pathway important in bone formation and result in a "low bone mass" phenotype. But how well does the skeleton bear functional loads; eg, how much do bones deform during loading and how resistant are bones to fracture? By systematic evaluation of bone morphological, densitometric, and mechanical properties, investigators can establish the "biomechanical mechanisms" whereby an experimental perturbation alters whole-bone mechanical function. The goal of this review is to clarify these biomechanical mechanisms and to make recommendations for systematically evaluating phenotypic changes in mouse bones, with a focus on long-bone diaphyses and cortical bone. Further, minimum reportable standards for testing conditions and outcome variables are suggested that will improve the comparison of data across studies. Basic biomechanical principles are reviewed, followed by a description of the cross-sectional morphological properties that best inform the net cellular effects of a given experimental perturbation and are most relevant to biomechanical function. Although morphology is critical, whole-bone mechanical properties can only be determined accurately by a mechanical test. The functional importance of stiffness, maximum load, postyield displacement, and work-to-fracture are reviewed. Because bone and body size are often strongly related, strategies to adjust whole-bone properties for body mass are detailed. Finally, a comprehensive framework is presented using real data, and several examples from the literature are reviewed to illustrate how to synthesize morphological, tissue-level, and whole-bone mechanical properties of mouse long bones.

SPINE: SParse eIgengene NEtwork linking gene expression clusters in Dehalococcoides mccartyi to perturbations in experimental conditions.

We present a statistical model designed to identify the effect of experimental perturbations on the aggregate behavior of the transcriptome expressed by the bacterium Dehalococcoides mccartyi strain 195. Strains of Dehalococcoides are used in sub-surface bioremediation applications because they organohalorespire tetrachloroethene and trichloroethene (common chlorinated solvents that contaminate the environment) to non-toxic ethene. However, the biochemical mechanism of this process remains incompletely described. Additionally, the response of Dehalococcoides to stress-inducing conditions that may be encountered at field-sites is not well understood. The constructed statistical model captured the aggregate behavior of gene expression phenotypes by modeling the distinct eigengenes of 100 transcript clusters, determining stable relationships among these clusters of gene transcripts with a sparse network-inference algorithm, and directly modeling the effect of changes in experimental conditions by constructing networks conditioned on the experimental state. Based on the model predictions, we discovered new response mechanisms for DMC, notably when the bacterium is exposed to solvent toxicity. The network identified a cluster containing thirteen gene transcripts directly connected to the solvent toxicity condition. Transcripts in this cluster include an iron-dependent regulator (DET0096-97) and a methylglyoxal synthase (DET0137). To validate these predictions, additional experiments were performed. Continuously fed cultures were exposed to saturating levels of tetrachloethene, thereby causing solvent toxicity, and transcripts that were predicted to be linked to solvent toxicity were monitored by quantitative reverse-transcription polymerase chain reaction. Twelve hours after being shocked with saturating levels of tetrachloroethene, the control transcripts (encoding for a key hydrogenase and the 16S rRNA) did not significantly change. By contrast, transcripts for DET0137 and DET0097 displayed a 46.8±11.5 and 14.6±9.3 fold up-regulation, respectively, supporting the model. This is the first study to identify transcripts in Dehalococcoides that potentially respond to tetrachloroethene solvent-toxicity conditions that may be encountered near contamination source zones in sub-surface environments.

Modulation of flagellum attachment zone protein FLAM3 and regulation of the cell shape in Trypanosoma brucei life cycle transitions.

ABSTRACTThe cell shape of Trypanosoma brucei is influenced by flagellum-to-cell-body attachment through a specialised structure – the flagellum attachment zone (FAZ). T. brucei exhibits numerous morphological forms during its life cycle and, at each stage, the FAZ length varies. We have analysed FLAM3, a large protein that localises to the FAZ region within the old and new flagellum. Ablation of FLAM3 expression causes a reduction in FAZ length; however, this has remarkably different consequences in the tsetse procyclic form versus the mammalian bloodstream form. In procyclic form cells FLAM3 RNAi results in the transition to an epimastigote-like shape, whereas in bloodstream form cells a severe cytokinesis defect associated with flagellum detachment is observed. Moreover, we demonstrate that the amount of FLAM3 and its localisation is dependent on ClpGM6 expression and vice versa. This evidence demonstrates that FAZ is a key regulator of trypanosome shape, with experimental perturbations being life cycle form dependent. An evolutionary cell biology explanation suggests that these differences are a reflection of the division process, the cytoskeleton and intrinsic structural plasticity of particular life cycle forms.

Analysis of soluble protein entry into primary cilia using semipermeabilized cells.

The primary cilium is a protrusion from the cell surface that serves as a specialized compartment for signal transduction. Many signaling factors are known to be dynamically concentrated within cilia and to require cilia for their function. Yet protein entry into primary cilia remains poorly understood. To enable a mechanistic analysis of soluble protein entry into cilia, we developed a method for semipermeabilization of mammalian cells in which the plasma membrane is permeabilized while the ciliary membrane remains intact. Using semipermeabilized cells as the basis for an in vitro diffusion-to-capture assay, we uncovered a size-dependent diffusion barrier that restricts soluble protein exchange between the cytosol and the cilium. The manipulability of this in vitro system enabled an extensive characterization of the ciliary diffusion barrier and led us to show that the barrier is mechanistically distinct from those at the axon initial segment and the nuclear pore complex. Because semipermeabilized cells enable a range of experimental perturbations that would not be easily feasible in intact cells, we believe this methodology will provide a unique resource for investigating primary cilium function in development and disease.

Macropinocytosis is the entry mechanism of amphotropic murine leukemia virus.

The entry mechanism of murine amphotropic retrovirus (A-MLV) has not been unambiguously determined. We show here that A-MLV is internalized not by caveolae or other pinocytic mechanisms but by macropinocytosis. Thus, A-MLV infection of mouse embryonic fibroblasts deficient for caveolin or dynamin, and NIH 3T3 cells knocked down for caveolin expression, was unaffected. Conversely, A-MLV infection of NIH 3T3 and HeLa cells was sensitive to amiloride analogues and actin-depolymerizing drugs that interfere with macropinocytosis. Further manipulation of the actin cytoskeleton through conditional expression of dominant positive or negative mutants of Rac1, PAK1, and RhoG, to increase or decrease macropinocytosis, similarly correlated with an augmented or inhibited infection with A-MLV, respectively. The same experimental perturbations affected the infection of viruses that use clathrin-coated-pit endocytosis or other pathways for entry only mildly or not at all. These data agree with immunofluorescence studies and cryo-immunogold labeling for electron microscopy, which demonstrate the presence of A-MLV in protrusion-rich areas of the cell surface and in cortical fluid phase (dextran)-filled macropinosomes, which also account for up to a half of the cellular uptake of the cell surface-binding lectin concanavalin A. We conclude that A-MLV use macropinocytosis as the predominant entry portal into cells. Binding and entry of virus particles into mammalian cells are the first steps of infection. Understanding how pathogens and toxins exploit or divert endocytosis pathways has advanced our understanding of membrane trafficking pathways, which benefits development of new therapeutic schemes and methods of drug delivery. We show here that amphotropic murine leukemia virus (A-MLV) pseudotyped with the amphotropic envelope protein (which expands the host range to many mammalian cells) gains entry into host cells by macropinocytosis. Macropinosomes form as large, fluid-filled vacuoles (up to 10 μm) following the collapse of cell surface protrusions and membrane scission. We used drugs or the introduction of mutant proteins that affect the actin cytoskeleton and cell surface dynamics to show that macropinocytosis and A-MLV infection are correlated, and we provide both light- and electron-microscopic evidence to show the localization of A-MLV in macropinosomes. Finally, we specifically exclude some other potential entry portals, including caveolae, previously suggested to internalize A-MLV.

How difficult is inference of mammalian causal gene regulatory networks?

Gene regulatory networks (GRNs) play a central role in systems biology, especially in the study of mammalian organ development. One key question remains largely unanswered: Is it possible to infer mammalian causal GRNs using observable gene co-expression patterns alone? We assembled two mouse GRN datasets (embryonic tooth and heart) and matching microarray gene expression profiles to systematically investigate the difficulties of mammalian causal GRN inference. The GRNs were assembled based on pieces of experimental genetic perturbation evidence from manually reading primary research articles. Each piece of perturbation evidence records the qualitative change of the expression of one gene following knock-down or over-expression of another gene. Our data have thorough annotation of tissue types and embryonic stages, as well as the type of regulation (activation, inhibition and no effect), which uniquely allows us to estimate both sensitivity and specificity of the inference of tissue specific causal GRN edges. Using these unprecedented datasets, we found that gene co-expression does not reliably distinguish true positive from false positive interactions, making inference of GRN in mammalian development very difficult. Nonetheless, if we have expression profiling data from genetic or molecular perturbation experiments, such as gene knock-out or signalling stimulation, it is possible to use the set of differentially expressed genes to recover causal regulatory relationships with good sensitivity and specificity. Our result supports the importance of using perturbation experimental data in causal network reconstruction. Furthermore, we showed that causal gene regulatory relationship can be highly cell type or developmental stage specific, suggesting the importance of employing expression profiles from homogeneous cell populations. This study provides essential datasets and empirical evidence to guide the development of new GRN inference methods for mammalian organ development.

Pattern Formation of an Epithelial Tubule by Mechanical Instability during Epididymal Development

A single epithelial tubule undergoes morphogenesis to form a functional shape during the development of internal organs; however, the mechanical processes that are directed by the molecular signals regulating tubular morphogenesis are poorly understood. Here, axial tubular buckling triggered by cell proliferation is shown to drive the morphogenesis of murine epididymal tubules through mechanical interactions between the developing epithelial tubule and its surrounding tissues. Through immunofluorescence labeling and mathematical modeling, epididymal tubule shape formation is found to depend on two factors: cell proliferation area in the tubule and mechanical resistance from the tissues surrounding the tubule. Moreover, experimental perturbations of these two factors alter the shape of the epididymal tubule as predicted by the mathematical model, suggesting that the shape of the epididymal tubule spontaneously emerges through mechanical coupling between developing tissues instead of by growing according to a predetermined fate.

The G ProteinαChaperone Ric-8 as a Potential Therapeutic Target

Resistance to inhibitors of cholinesterase (Ric-8)A and Ric-8B are essential genes that encode positive regulators of heterotrimeric G protein α subunits. Controversy persists surrounding the precise way(s) that Ric-8 proteins affect G protein biology and signaling. Ric-8 proteins chaperone nucleotide-free Gα-subunit states during biosynthetic protein folding prior to G protein heterotrimer assembly. In organisms spanning the evolutionary window of Ric-8 expression, experimental perturbation of Ric-8 genes results in reduced functional abundances of G proteins because G protein α subunits are misfolded and degraded rapidly. Ric-8 proteins also act as Gα-subunit guanine nucleotide exchange factors (GEFs) in vitro. However, Ric-8 GEF activity could strictly be an in vitro phenomenon stemming from the ability of Ric-8 to induce partial Gα unfolding, thereby enhancing GDP release. Ric-8 GEF activity clearly differs from the GEF activity of G protein-coupled receptors (GPCRs). G protein βγ is inhibitory to Ric-8 action but obligate for receptors. It remains an open question whether Ric-8 has dual functions in cells and regulates G proteins as both a molecular chaperone and GEF. Clearly, Ric-8 has a profound influence on heterotrimeric G protein function. For this reason, we propose that Ric-8 proteins are as yet untested therapeutic targets in which pharmacological inhibition of the Ric-8/Gα protein-protein interface could serve to attenuate the effects of disease-causing G proteins (constitutively active mutants) and/or GPCR signaling. This minireview will chronicle the understanding of Ric-8 function, provide a comparative discussion of the Ric-8 molecular chaperoning and GEF activities, and support the case for why Ric-8 proteins should be considered potential targets for development of new therapies.

Abstract 1543: Systematic functional assessment of driver mutations for individual ovarian cancers

Abstract Although epithelial ovarian cancer is an intractable cancer and the overall five-year survival rate is approximately 40 to 50%, clinically effective molecular-targeted therapeutics to treat ovarian cancer has not been available thus far. Recent progress of massively parallel sequencing technologies has been enabling us to identify substantial number of genomic aberrations including somatic mutations, indels and translocations in variable cancer types exemplified by epithelial ovarian cancer and non-small cell lung cancer. However, it is still very difficult to identify a ‘driver’ gene, which is functionally relevant for cancer growth, from hundreds of mutations that are usually found from a single cancer sample of an individual patient. An in vitro cultured cell line, which is originated from a clinical sample, can be utilized to assess functions of the genes that are identified to have somatic aberrations in the paired clinical tumor. We have been establishing multiple ovarian cancer cell lines derived from clinical samples and analyzing pairs of the cell lines and clinical samples with exon-capture sequencing. The genes, which are somatically mutated and shared in a pair of cell line and clinical sample, are targeted with RNA interference in the cell line, and subsequently the effect on cell proliferation and survival is examined. As such, a driver mutation for an individual cancer can be identified by a set of experimental perturbations. Identified driver mutations would be then validated with publicly available databases and with clinical archive samples. As a proof-of-principle, somatic mutations of PIK3CA and KRAS were detected as drivers to promote cell proliferation of a carcinosarcoma. We believe this approach can increase an opportunity to identify molecular targets in ovarian cancer therapeutics. Citation Format: Seiichi Mori, Takako Yokomizo, Takeshi Fujiwara, Osamu Goto, Haruko Iwase, Kazuyoshi Kato, Tokuichi Kawaguchi, Tetsuo Noda. Systematic functional assessment of driver mutations for individual ovarian cancers. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 1543. doi:10.1158/1538-7445.AM2014-1543

Optimal Design of Non-equilibrium Experiments for Genetic Network Interrogation.

Many experimental systems in biology, especially synthetic gene networks, are amenable to perturbations that are controlled by the experimenter. We developed an optimal design algorithm that calculates optimal observation times in conjunction with optimal experimental perturbations in order to maximize the amount of information gained from longitudinal data derived from such experiments. We applied the algorithm to a validated model of a synthetic Brome Mosaic Virus (BMV) gene network and found that optimizing experimental perturbations may substantially decrease uncertainty in estimating BMV model parameters.

Specifying informative experiment stimulation conditions for resolving dynamical uncertainty in biological systems.

A computationally efficient model-based design of experiments (MBDOE) strategy is developed to plan an optimal experiment by specifying the experimental stimulation magnitudes and measurement points. The strategy is extended from previous work which optimized the experimental design over a space of measurable species and time points. We include system inputs (stimulation conditions) in the experiment design search to investigate if the addition of perturbations enhances the ability of the MBDOE method to resolve uncertainties in system dynamics. The MBDOE problem is made computationally tractable by using a sparse-grid approximation of the model output dynamics, pre-specifying the time points at which the input or experimental perturbations can be applied, and creating scenario trees to explore the endogenous uncertainty. Consecutive scenario trees are used to determine the best input magnitudes and select the optimal associated measurement species and time points. We demonstrate the effectiveness of this strategy on a T-Cell Receptor (TCR) signaling pathway model.

Abraham Guz (1929-2014).

How do you describe a giant? As a person of abnormally great stature; above the usual size; or someone with much greater powers than the average person? Abraham Guz (‘Abe’) was not a tall man, but he was a scientific giant with an absolute desire to understand and find the scientific truth. This drive led him to make many original contributions in respiratory and cardiovascular physiology. Indeed, a simple literature search produces over 250 peer reviewed papers spanning some 50 years, the top 25 of which have been quoted over 2500 times. Moreover, most of his outstanding work was published in The Journal of Physiology, and we received it with gratitude because it was always insightful and ahead of its time. Many will argue that counting citations is a poor measure of scientific endeavour and I suspect Abe himself would certainly agree. However, this approach did provide us with a place to start when faced with the task of bringing together Abe's huge and varied scientific contributions.

Temperature–parasitism synergy alters intertidal soft-bottom community structure

The abundant amphipod Corophium volutator is an ecosystem engineer in soft-bottom intertidal communities due to its grazing and bioturbation activity. However, the amphipod commonly serves as second intermediate host for detrimental microphallid trematodes. This host–parasite association is potentially very sensitive to climate change as the transmission rate of larval trematodes (cercariae) from first intermediate hosts (mud snails) to amphipods, and hence infection intensity-dependent amphipod mortality, generally accelerates with increasing temperature. Given the ecosystem engineering role of C. volutator, we hypothesized that elevated temperature, indirectly through increased parasite-induced amphipod mortality, significantly affects the structure of the surrounding benthic community. To test this, we performed a three-month outdoor mesocosm experiment exposing a natural soft-bottom community of benthic plants and animals to four different treatments: relatively low mean water temperature (18°C) with low (<4%) and high (c. 31%) trematode prevalence in the snail population, and relatively high mean water temperature (22°C) with low and high trematode prevalence. Both temperature and snail parasitism had a significant impact on amphipod abundance and resulted in their almost complete eradication at the warm temperature with high parasitism. Aside from the amphipod hosts, increased temperature affected the abundance of eight faunal species out of 22 in total (36%), whereas increased level of parasitism or the parasite–temperature interaction influenced the density of three species (14%). The treatment-determined Corophium abundance played an isolated role for several species, particularly the polychaetes Hediste diversicolor (positive) and Polydora ligni (negative). Regarding primary producers, the overall frequency distribution of 21 species of benthic diatoms differed between all four treatments. Species-specific effects were few and weak, however, and microalgal abundance (chlorophyll-a) was statistically unaffected by treatments, together suggesting a rather resilient plant community towards the experimental perturbations. At the community level, the non-host macrofaunal diversity decreased with temperature in the low parasitism treatments whereas it increased with temperature in the high parasitism treatments. This suggests that sufficiently high levels of parasitism may turn a negative effect of elevated temperature on macrofaunal biodiversity into a positive effect. Hence, our experiment demonstrates that the synergy between parasitism and a relevant temperature increase, e.g. in lieu of climate oscillations or global warming, may have broad ecological consequences for the organization and function of soft-bottom communities, in part through elevated parasite-induced mortality of a central community member.

Yeast mating and image-based quantification of spatial pattern formation.

Communication between cells is a ubiquitous feature of cell populations and is frequently realized by secretion and detection of signaling molecules. Direct visualization of the resulting complex gradients between secreting and receiving cells is often impossible due to the small size of diffusing molecules and because such visualization requires experimental perturbations such as attachment of fluorescent markers, which can change diffusion properties. We designed a method to estimate such extracellular concentration profiles in vivo by using spatiotemporal mathematical models derived from microscopic analysis. This method is applied to populations of thousands of haploid yeast cells during mating in order to quantify the extracellular distributions of the pheromone α-factor and the activity of the aspartyl protease Bar1. We demonstrate that Bar1 limits the range of the extracellular pheromone signal and is critical in establishing α-factor concentration gradients, which is crucial for effective mating. Moreover, haploid populations of wild type yeast cells, but not BAR1 deletion strains, create a pheromone pattern in which cells differentially grow and mate, with low pheromone regions where cells continue to bud and regions with higher pheromone levels and gradients where cells conjugate to form diploids. However, this effect seems to be exclusive to high-density cultures. Our results show a new role of Bar1 protease regulating the pheromone distribution within larger populations and not only locally inside an ascus or among few cells. As a consequence, wild type populations have not only higher mating efficiency, but also higher growth rates than mixed MAT a bar1Δ/MAT α cultures. We provide an explanation of how a rapidly diffusing molecule can be exploited by cells to provide spatial information that divides the population into different transcriptional programs and phenotypes.

Mechanochemical regulation of oscillatory follicle cell dynamics in the developing Drosophila egg chamber.

During tissue elongation from stage 9 to stage 10 in Drosophila oogenesis, the egg chamber increases in length by ∼1.7-fold while increasing in volume by eightfold. During these stages, spontaneous oscillations in the contraction of cell basal surfaces develop in a subset of follicle cells. This patterned activity is required for elongation of the egg chamber; however, the mechanisms generating the spatiotemporal pattern have been unclear. Here we use a combination of quantitative modeling and experimental perturbation to show that mechanochemical interactions are sufficient to generate oscillations of myosin contractile activity in the observed spatiotemporal pattern. We propose that follicle cells in the epithelial layer contract against pressure in the expanding egg chamber. As tension in the epithelial layer increases, Rho kinase signaling activates myosin assembly and contraction. The activation process is cooperative, leading to a limit cycle in the myosin dynamics. Our model produces asynchronous oscillations in follicle cell area and myosin content, consistent with experimental observations. In addition, we test the prediction that removal of the basal lamina will increase the average oscillation period. The model demonstrates that in principle, mechanochemical interactions are sufficient to drive patterning and morphogenesis, independent of patterned gene expression.

Evolutionary developmental transcriptomics reveals a gene network module regulating interspecific diversity in plant leaf shape.

Despite a long-standing interest in the genetic basis of morphological diversity, the molecular mechanisms that give rise to developmental variation are incompletely understood. Here, we use comparative transcriptomics coupled with the construction of gene coexpression networks to predict a gene regulatory network (GRN) for leaf development in tomato and two related wild species with strikingly different leaf morphologies. The core network in the leaf developmental GRN contains regulators of leaf morphology that function in global cell proliferation with peripheral gene network modules (GNMs). The BLADE-ON-PETIOLE (BOP) transcription factor in one GNM controls the core network by altering effective concentration of the KNOTTED-like HOMEOBOX gene product. Comparative network analysis and experimental perturbations of BOP levels suggest that variation in BOP expression could explain the diversity in leaf complexity among these species through dynamic rewiring of interactions in the GRN. The peripheral location of the BOP-containing GNM in the leaf developmental GRN and the phenotypic mimics of evolutionary diversity caused by alteration in BOP levels identify a key role for this GNM in canalizing the leaf morphospace by modifying the maturation schedule of leaves to create morphological diversity.