Multiple Copies of flhDC in Paraburkholderia unamae Regulate Flagellar Gene Expression, Motility, and Biofilm Formation.

FlhDC is a heterohexameric complex that acts as a master regulator of flagellar biosynthesis genes in numerous bacteria. Previous studies have identified a single flhDC operon encoding this complex. However, we found that two flhDC loci are present throughout Paraburkholderia, and two additional flhC copies are also present in Paraburkholderia unamae. Systematic deletion analysis in P. unamae of the different flhDC copies showed that one of the operons, flhDC1, plays the predominant role, with deletion of its genes resulting in a severe inhibition of motility and biofilm formation. Expression analysis using promoter-lacZ fusions and real-time quantitative PCR support the primary role of flhDC1 in flagellar gene regulation, with flhDC2 a secondary contributor. Phylogenetic analysis shows the presence of the flhDC1 and flhDC2 operons throughout Paraburkholderia. In contrast, Burkholderia and other bacteria only carry the copy syntenous with flhDC2. The variations in impact each copy of flhDC has on downstream processes indicate that regulation of FlhDC in P. unamae, and likely other Paraburkholderia species, is regulated at least in part by the presence of multiple copies of these genes. IMPORTANCE Motility is important in the colonization of plant roots by beneficial and pathogenic bacteria, with flagella playing essential roles in host cell adhesion, entrance, and biofilm formation. Flagellar biosynthesis is energetically expensive. Its complex regulation by the FlhDC master regulator is well studied in peritrichous flagella expressing enterics. We report the unique presence throughout Paraburkholderia of multiple copies of flhDC. In P. unamae, the flhDC1 copy showed higher expression and a greater effect on swim motility, flagellar development, and regulation of downstream genes, than the flhDC2 copy that is syntenous to flhDC in Escherichia coli and pathogenic Burkholderia spp. The flhDC genes have evolved differently in these plant-growth-promoting bacteria, giving an additional layer of complexity in gene regulation by FlhDC.

Genetic interaction profiles of regulatory kinases differ between environmental conditions and cellular states.

Cell growth and quiescence in eukaryotic cells is controlled by an evolutionarily conserved network of signaling pathways. Signal transduction networks operate to modulate a wide range of cellular processes and physiological properties when cells exit proliferative growth and initiate a quiescent state. How signaling networks function to respond to diverse signals that result in cell cycle exit and establishment of a quiescent state is poorly understood. Here, we studied the function of signaling pathways in quiescent cells using global genetic interaction mapping in the model eukaryotic cell, Saccharomyces cerevisiae (budding yeast). We performed pooled analysis of genotypes using molecular barcode sequencing (Bar-seq) to test the role of ~4,000 gene deletion mutants and ~12,000 pairwise interactions between all non-essential genes and the protein kinase genes TOR1, RIM15, and PHO85 in three different nutrient-restricted conditions in both proliferative and quiescent cells. We detect up to 10-fold more genetic interactions in quiescent cells than proliferative cells. We find that both individual gene effects and genetic interaction profiles vary depending on the specific pro-quiescence signal. The master regulator of quiescence, RIM15, shows distinct genetic interaction profiles in response to different starvation signals. However, vacuole-related functions show consistent genetic interactions with RIM15 in response to different starvation signals, suggesting that RIM15 integrates diverse signals to maintain protein homeostasis in quiescent cells. Our study expands genome-wide genetic interaction profiling to additional conditions, and phenotypes, and highlights the conditional dependence of epistasis.

A complete statistical model for calibration of RNA-seq counts using external spike-ins and maximum likelihood theory.

A fundamental assumption, common to the vast majority of high-throughput transcriptome analyses, is that the expression of most genes is unchanged among samples and that total cellular RNA remains constant. As the number of analyzed experimental systems increases however, different independent studies demonstrate that this assumption is often violated. We present a calibration method using RNA spike-ins that allows for the measurement of absolute cellular abundance of RNA molecules. We apply the method to pooled RNA from cell populations of known sizes. For each transcript, we compute a nominal abundance that can be converted to absolute by dividing by a scale factor determined in separate experiments: the yield coefficient of the transcript relative to that of a reference spike-in measured with the same protocol. The method is derived by maximum likelihood theory in the context of a complete statistical model for sequencing counts contributed by cellular RNA and spike-ins. The counts are based on a sample from a fixed number of cells to which a fixed population of spike-in molecules has been added. We illustrate and evaluate the method with applications to two global expression data sets, one from the model eukaryote Saccharomyces cerevisiae, proliferating at different growth rates, and differentiating cardiopharyngeal cell lineages in the chordate Ciona robusta. We tested the method in a technical replicate dilution study, and in a k-fold validation study.

Single-cell copy number variant detection reveals the dynamics and diversity of adaptation.

Copy number variants (CNVs) are a pervasive source of genetic variation and evolutionary potential, but the dynamics and diversity of CNVs within evolving populations remain unclear. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics with which they are generated, selected, and maintained. Here, we developed a fluorescent CNV reporter to detect de novo gene amplifications and deletions in individual cells. We used the CNV reporter in Saccharomyces cerevisiae to study CNV formation at the GAP1 locus, which encodes the general amino acid permease, in different nutrient-limited chemostat conditions. We find that under strong selection, GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution, resulting in predictable dynamics. Molecular characterization of CNV-containing lineages shows that the CNV reporter detects different classes of CNVs, including aneuploidies, nonreciprocal translocations, tandem duplications, and complex CNVs. Despite GAP1's proximity to repeat sequences that facilitate intrachromosomal recombination, breakpoint analysis revealed that short inverted repeat sequences mediate formation of at least 50% of GAP1 CNVs. Inverted repeat sequences are also found at breakpoints at the DUR3 locus, where CNVs are selected in urea-limited chemostats. Analysis of 28 CNV breakpoints indicates that inverted repeats are typically 8 nucleotides in length and separated by 40 bases. The features of these CNVs are consistent with origin-dependent inverted-repeat amplification (ODIRA), suggesting that replication-based mechanisms of CNV formation may be a common source of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 102-104 independent CNV-containing lineages initially compete within populations, resulting in extreme clonal interference. However, only a small number (18-21) of CNV lineages ever constitute more than 1% of the CNV subpopulation, and as selection progresses, the diversity of CNV lineages declines. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into their dynamics, diversity, and formation mechanisms in the context of adaptive evolution.

Systematic identification of factors mediating accelerated mRNA degradation in response to changes in environmental nitrogen.

Cellular responses to changing environments frequently involve rapid reprogramming of the transcriptome. Regulated changes in mRNA degradation rates can accelerate reprogramming by clearing or stabilizing extant transcripts. Here, we measured mRNA stability using 4-thiouracil labeling in the budding yeast Saccharomyces cerevisiae during a nitrogen upshift and found that 78 mRNAs are subject to destabilization. These transcripts include Nitrogen Catabolite Repression (NCR) and carbon metabolism mRNAs, suggesting that mRNA destabilization is a mechanism for targeted reprogramming of the transcriptome. To explore the molecular basis of destabilization we implemented a SortSeq approach to screen the pooled deletion collection library for trans factors that mediate rapid GAP1 mRNA repression. We combined low-input multiplexed Barcode sequencing with branched-DNA single-molecule mRNA FISH and Fluorescence-activated cell sorting (BFF) to identify the Lsm1-7p/Pat1p complex and general mRNA decay machinery as important for GAP1 mRNA clearance. We also find that the decapping modulators EDC3 and SCD6, translation factor eIF4G2, and the 5' UTR of GAP1 are factors that mediate rapid repression of GAP1 mRNA, suggesting that translational control may impact the post-transcriptional fate of mRNAs in response to environmental changes.

An incoherent feedforward loop facilitates adaptive tuning of gene expression.

We studied adaptive evolution of gene expression using long-term experimental evolution of Saccharomyces cerevisiae in ammonium-limited chemostats. We found repeated selection for non-synonymous variation in the DNA binding domain of the transcriptional activator, GAT1, which functions with the repressor, DAL80 in an incoherent type-1 feedforward loop (I1-FFL) to control expression of the high affinity ammonium transporter gene, MEP2. Missense mutations in the DNA binding domain of GAT1 reduce its binding to the GATAA consensus sequence. However, we show experimentally, and using mathematical modeling, that decreases in GAT1 binding result in increased expression of MEP2 as a consequence of properties of I1-FFLs. Our results show that I1-FFLs, one of the most commonly occurring network motifs in transcriptional networks, can facilitate adaptive tuning of gene expression through modulation of transcription factor binding affinities. Our findings highlight the importance of gene regulatory architectures in the evolution of gene expression.

Steady-state and dynamic gene expression programs in Saccharomyces cerevisiae in response to variation in environmental nitrogen.

Cell growth rate is regulated in response to the abundance and molecular form of essential nutrients. InSaccharomyces cerevisiae(budding yeast), the molecular form of environmental nitrogen is a major determinant of cell growth rate, supporting growth rates that vary at least threefold. Transcriptional control of nitrogen use is mediated in large part by nitrogen catabolite repression (NCR), which results in the repression of specific transcripts in the presence of a preferred nitrogen source that supports a fast growth rate, such as glutamine, that are otherwise expressed in the presence of a nonpreferred nitrogen source, such as proline, which supports a slower growth rate. Differential expression of the NCR regulon and additional nitrogen-responsive genes results in >500 transcripts that are differentially expressed in cells growing in the presence of different nitrogen sources in batch cultures. Here we find that in growth rate-controlled cultures using nitrogen-limited chemostats, gene expression programs are strikingly similar regardless of nitrogen source. NCR expression is derepressed in all nitrogen-limiting chemostat conditions regardless of nitrogen source, and in these conditions, only 34 transcripts exhibit nitrogen source-specific differential gene expression. Addition of either the preferred nitrogen source, glutamine, or the nonpreferred nitrogen source, proline, to cells growing in nitrogen-limited chemostats results in rapid, dose-dependent repression of the NCR regulon. Using a novel means of computational normalization to compare global gene expression programs in steady-state and dynamic conditions, we find evidence that the addition of nitrogen to nitrogen-limited cells results in the transient overproduction of transcripts required for protein translation. Simultaneously, we find that that accelerated mRNA degradation underlies the rapid clearing of a subset of transcripts, which is most pronounced for the highly expressed NCR-regulated permease genesGAP1,MEP2,DAL5,PUT4, andDIP5 Our results reveal novel aspects of nitrogen-regulated gene expression and highlight the need for a quantitative approach to study how the cell coordinates protein translation and nitrogen assimilation to optimize cell growth in different environments.

The use of chemostats in microbial systems biology.

Cells regulate their rate of growth in response to signals from the external world. As the cell grows, diverse cellular processes must be coordinated including macromolecular synthesis, metabolism and ultimately, commitment to the cell division cycle. The chemostat, a method of experimentally controlling cell growth rate, provides a powerful means of systematically studying how growth rate impacts cellular processes - including gene expression and metabolism - and the regulatory networks that control the rate of cell growth. When maintained for hundreds of generations chemostats can be used to study adaptive evolution of microbes in environmental conditions that limit cell growth. We describe the principle of chemostat cultures, demonstrate their operation and provide examples of their various applications. Following a period of disuse after their introduction in the middle of the twentieth century, the convergence of genome-scale methodologies with a renewed interest in the regulation of cell growth and the molecular basis of adaptive evolution is stimulating a renaissance in the use of chemostats in biological research.

System-level analysis of genes and functions affecting survival during nutrient starvation in Saccharomyces cerevisiae.

An essential property of all cells is the ability to exit from active cell division and persist in a quiescent state. For single-celled microbes this primarily occurs in response to nutrient deprivation. We studied the genetic requirements for survival of Saccharomyces cerevisiae when starved for either of two nutrients: phosphate or leucine. We measured the survival of nearly all nonessential haploid null yeast mutants in mixed populations using a quantitative sequencing method that estimates the abundance of each mutant on the basis of frequency of unique molecular barcodes. Starvation for phosphate results in a population half-life of 337 hr whereas starvation for leucine results in a half-life of 27.7 hr. To measure survival of individual mutants in each population we developed a statistical framework that accounts for the multiple sources of experimental variation. From the identities of the genes in which mutations strongly affect survival, we identify genetic evidence for several cellular processes affecting survival during nutrient starvation, including autophagy, chromatin remodeling, mRNA processing, and cytoskeleton function. In addition, we found evidence that mitochondrial and peroxisome function is required for survival. Our experimental and analytical methods represent an efficient and quantitative approach to characterizing genetic functions and networks with unprecedented resolution and identified genotype-by-environment interactions that have important implications for interpretation of studies of aging and quiescence in yeast.

Rett syndrome: clinical manifestations in males with MECP2 mutations.

Rett syndrome is a neurodevelopmental disorder characterized by cognitive and adaptive regression with autistic features, loss of acquired skills, and stereotypic hand movements that almost exclusively affects females. It is an X-linked dominant disorder, with presumed lethality in males. Nonetheless, there are a few descriptions of males suspected of having Rett syndrome. With the recent discovery that the MECP2 gene is responsible for most cases of Rett syndrome, it is possible to molecularly assess cases of affected males by direct sequencing analysis. We describe an Israeli family consisting of a female having classic Rett syndrome and a male sibling with severe neonatal encephalopathy. Molecular analysis revealed that both sister and brother have the same MECP2 gene mutation; however, their mother does not. This case, as well as other published studies of males with MECP2 mutations, reveals that the clinical manifestations in viable males vary from neonates with severe encephalopathy to adults with mental retardation and demonstrate genotype-phenotype correlations.