Diversification of small RNA pathways underlies germline RNA interference incompetence in wild Caenorhabditis elegans strains.

The discovery that experimental delivery of dsRNA can induce gene silencing at target genes revolutionized genetics research, by both uncovering essential biological processes and creating new tools for developmental geneticists. However, the efficacy of exogenous RNA interference (RNAi) varies dramatically within the Caenorhabditis elegans natural population, raising questions about our understanding of RNAi in the lab relative to its activity and significance in nature. Here, we investigate why some wild strains fail to mount a robust RNAi response to germline targets. We observe diversity in mechanism: in some strains, the response is stochastic, either on or off among individuals, while in others, the response is consistent but delayed. Increased activity of the Argonaute PPW-1, which is required for germline RNAi in the laboratory strain N2, rescues the response in some strains but dampens it further in others. Among wild strains, genes known to mediate RNAi exhibited very high expression variation relative to other genes in the genome as well as allelic divergence and strain-specific instances of pseudogenization at the sequence level. Our results demonstrate functional diversification in the small RNA pathways in C. elegans and suggest that RNAi processes are evolving rapidly and dynamically in nature.

Beyond the reference: gene expression variation and transcriptional response to RNA interference in Caenorhabditis elegans.

Though natural systems harbor genetic and phenotypic variation, research in model organisms is often restricted to a reference strain. Focusing on a reference strain yields a great depth of knowledge but potentially at the cost of breadth of understanding. Furthermore, tools developed in the reference context may introduce bias when applied to other strains, posing challenges to defining the scope of variation within model systems. Here, we evaluate how genetic differences among 5 wild Caenorhabditis elegans strains affect gene expression and its quantification, in general and after induction of the RNA interference (RNAi) response. Across strains, 34% of genes were differentially expressed in the control condition, including 411 genes that were not expressed at all in at least 1 strain; 49 of these were unexpressed in reference strain N2. Reference genome mapping bias caused limited concern: despite hyperdiverse hotspots throughout the genome, 92% of variably expressed genes were robust to mapping issues. The transcriptional response to RNAi was highly strain- and target-gene-specific and did not correlate with RNAi efficiency, as the 2 RNAi-insensitive strains showed more differentially expressed genes following RNAi treatment than the RNAi-sensitive reference strain. We conclude that gene expression, generally and in response to RNAi, differs across C. elegans strains such that the choice of strain may meaningfully influence scientific inferences. Finally, we introduce a resource for querying gene expression variation in this dataset at https://wildworm.biosci.gatech.edu/rnai/.

Beyond the reference: gene expression variation and transcriptional response to RNAi in C. elegans.

A universal feature of living systems is that natural variation in genotype underpins variation in phenotype. Yet, research in model organisms is often constrained to a single genetic background, the reference strain. Further, genomic studies that do evaluate wild strains typically rely on the reference strain genome for read alignment, leading to the possibility of biased inferences based on incomplete or inaccurate mapping; the extent of reference bias can be difficult to quantify. As an intermediary between genome and organismal traits, gene expression is well positioned to describe natural variability across genotypes generally and in the context of environmental responses, which can represent complex adaptive phenotypes. C. elegans sits at the forefront of investigation into small-RNA gene regulatory mechanisms, or RNA interference (RNAi), and wild strains exhibit natural variation in RNAi competency following environmental triggers. Here, we examine how genetic differences among five wild strains affect the C. elegans transcriptome in general and after inducing RNAi responses to two germline target genes. Approximately 34% of genes were differentially expressed across strains; 411 genes were not expressed at all in at least one strain despite robust expression in others, including 49 genes not expressed in reference strain N2. Despite the presence of hyper-diverse hotspots throughout the C. elegans genome, reference mapping bias was of limited concern: over 92% of variably expressed genes were robust to mapping issues. Overall, the transcriptional response to RNAi was strongly strain-specific and highly specific to the target gene, and the laboratory strain N2 was not representative of the other strains. Moreover, the transcriptional response to RNAi was not correlated with RNAi phenotypic penetrance; the two germline RNAi incompetent strains exhibited substantial differential gene expression following RNAi treatment, indicating an RNAi response despite failure to reduce expression of the target gene. We conclude that gene expression, both generally and in response to RNAi, differs across C. elegans strains such that choice of strain may meaningfully influence scientific conclusions. To provide a public, easily accessible resource for querying gene expression variation in this dataset, we introduce an interactive website at https://wildworm.biosci.gatech.edu/rnai/ .

Single-molecule FISH in C. elegans embryos reveals early embryonic expression dynamics of par-2 , lgl-1 and chin-1 and possible differences between hyper-diverged strains.

Wild C. elegans strains harbor natural variation in developmental pathways, but investigating these differences requires precise and well-powered phenotyping methods. Here we employ a microfluidics platform for single-molecule FISH to simultaneously visualize the transcripts of three genes in embryos of two distinct strains. We capture transcripts at high resolution by developmental stage in over one hundred embryos of each strain and observe wide-scale conservation of expression, but subtle differences in par-2 and chin-1 abundance and rate of change. As both genes reside in a genomic interval of hyper-divergence, these results may reflect consequences of pathway evolution over long timescales.

High-Temporal-Resolution smFISH Method for Gene Expression Studies in Caenorhabditis elegans Embryos.

Recent development in fluorescence-based molecular tools has contributed significantly to developmental studies, including embryogenesis. Many of these tools rely on multiple steps of sample manipulation, so obtaining large sample sizes presents a major challenge as it can be labor-intensive and time-consuming. However, large sample sizes are required to uncover critical aspects of embryogenesis, for example, subtle phenotypic differences or gene expression dynamics. This problem is particularly relevant for single-molecule fluorescence in situ hybridization (smFISH) studies in Caenorhabditis elegans embryogenesis. Microfluidics can help address this issue by allowing a large number of samples and parallelization of experiments. However, performing efficient reagent exchange on chip for large numbers of embryos remains a bottleneck. Here, we present a microfluidic pipeline for large-scale smFISH imaging of C. elegans embryos with minimized labor. We designed embryo traps and engineered a protocol allowing for efficient chemical exchange for hundreds of C. elegans embryos simultaneously. Furthermore, the device design and small footprint optimize imaging throughput by facilitating spatial registration and enabling minimal user input. We conducted the smFISH protocol on chip and demonstrated that image quality is preserved. With one device replacing the equivalent of 10 glass slides of embryos mounted manually, our microfluidic approach greatly increases throughput. Finally, to highlight the capability of our platform to perform longitudinal studies with high temporal resolution, we conducted a temporal analysis of par-1 gene expression in early C. elegans embryos. The method demonstrated here paves the way for systematic high-temporal-resolution studies that will benefit large-scale RNAi and drug screens and in systems beyond C. elegans embryos.

Molecular characterization of lysR-lysXE, gcdR-gcdHG and amaR-amaAB operons for lysine export and catabolism: a comprehensive lysine catabolic network in Pseudomonas aeruginosa PAO1.

Among multiple interconnected pathways for l-Lysine catabolism in pseudomonads, it has been reported that Pseudomonas aeruginosa PAO1 employs the decarboxylase and the transaminase pathways. However, up until now, knowledge of several genes involved in operation and regulation of these pathways was still missing. Transcriptome analyses coupled with promoter activity measurements and growth phenotype analyses led us to identify new members in l-Lys and d-Lys catabolism and regulation, including gcdR-gcdHG for glutarate utilization, dpkA, amaR-amaAB and PA2035 for d-Lys catabolism, lysR-lysXE for putative l-Lys efflux and lysP for putative l-Lys uptake. The gcdHG operon encodes an acyl-CoA transferase (gcdG) and glutaryl-CoA dehydrogenase (gcdH) and is under the control of the transcriptional activator GcdR. Growth on l-Lys was enhanced in the mutants of lysX and lysE, supporting the operation of l-Lys efflux. The transcriptional activator LysR is responsible for l-Lys specific induction of lysXE and the PA4181-82 operon of unknown function. The putative operator sites of GcdR and LysR were deduced from serial deletions and comparative genomic sequence analyses, and the formation of nucleoprotein complexes was demonstrated with purified His-tagged GcdR and LysR. The amaAB operon encodes two enzymes to convert pipecolate to 2-aminoadipate. Induction of the amaAB operon by l-Lys, d-Lys and pipecolate requires a functional AmaR, supporting convergence of Lys catabolic pathways to pipecolate. Growth on pipecolate was retarded in the gcdG and gcdH mutants, suggesting the importance of glutarate in pipecolate and 2-aminoadipate utilization. Furthermore, this study indicated links in the control of interconnected networks of lysine and arginine catabolism in P. aeruginosa.

HES-Mediated Repression of Pten in Caenorhabditis elegans.

The hairy/enhancer-of-split (HES) group of transcription factors controls embryonic development, often by acting downstream of the Notch signaling pathway; however, little is known about postembryonic roles of these proteins. In Caenorhabditis elegans, the six proteins that make up the REF-1 family are considered to be HES orthologs that act in both Notch-dependent and Notch-independent pathways to regulate embryonic events. To further our understanding of how the REF-1 family works to coordinate postembryonic cellular events, we performed a functional characterization of the REF-1 family member, HLH-25. We show that, after embryogenesis, hlh-25 expression persists throughout every developmental stage, including dauer, into adulthood. Like animals that carry loss-of-function alleles in genes required for normal cell-cycle progression, the phenotypes of hlh-25 animals include reduced brood size, unfertilized oocytes, and abnormal gonad morphology. Using gene expression microarray, we show that the HLH-25 transcriptional network correlates with the phenotypes of hlh-25 animals and that the C. elegans Pten ortholog, daf-18, is one major hub in the network. Finally, we show that HLH-25 regulates C. elegans lifespan and dauer recovery, which correlates with a role in the transcriptional repression of daf-18 activity. Collectively, these data provide the first genetic evidence that HLH-25 may be a functional ortholog of mammalian HES1, which represses PTEN activity in mice and human cells.

Functional characterization of the agtABCD and agtSR operons for 4-aminobutyrate and 5-aminovalerate uptake and regulation in Pseudomonas aeruginosa PAO1.

Growth of Pseudomonas aeruginosa on diamines cadaverine, putrescine, and diaminopropane requires the γ-glutamylation pathway to convert these diamines into δ-aminovalerate (AMV), γ-aminobutyrate (GABA), and ß-alanine. From DNA microarrays experiments the agtABCD operon (PA0603-0606) encoding components for an ABC transporter system was found inducible by exogenous AMV, GABA, and ß-alanine, but not by diamines. Induction of the agtABCD operon was abolished in the mutants of upstream agtS (PA0600) or agtR (PA0601) genes encoding the membrane-anchored sensor and the response regulator of a two-component regulatory system, respectively. Growth phenotype analysis supports the physiological functions of these agt genes on utilization of AMV and GABA. Through measurements of ß-galactosidase activities from an agtA::lacZ fusion, the requirement of a functional AgtS in control of the induction effect by exogenous AMV and GABA was further substantiated. The recombinant hexa-hisidine tagged agtR was constructed and purified to demonstrate its specific interactions with the agtA promoter region by electrophoretic mobility shift assays. In summary, this study establishes the functions of agtSR and agtABCD operons in AMV and GABA uptake, and provides a potential linkage between AMV/GABA metabolism and polymicrobial infection through the recently reported function of agtR in sensing of peptidoglycan shed by gram-positive bacteria (Korgaonkar et al., Proc Natl Acad Sci USA 110:1059-1064, 2013).

Molecular characterization of PauR and its role in control of putrescine and cadaverine catabolism through the γ-glutamylation pathway in Pseudomonas aeruginosa PAO1.

Pseudomonas aeruginosa PAO1 grows on a variety of polyamines as the sole source of carbon and nitrogen. Catabolism of polyamines is mediated by the γ-glutamylation pathway, which is complicated by the existence of multiple homologous enzymes with redundant specificities toward different polyamines for a more diverse metabolic capacity in this organism. Through a series of markerless gene knockout mutants and complementation tests, specific combinations of pauABCD (polyamine utilization) genes were deciphered for catabolism of different polyamines. Among six pauA genes, expression of pauA1, pauA2, pauA4, and pauA5 was found to be inducible by diamines putrescine (PUT) and cadaverine (CAD) but not by diaminopropane. Activation of these promoters was regulated by the PauR repressor, as evidenced by constitutively active promoters in the pauR mutant. The activities of these promoters were further enhanced by exogenous PUT or CAD in the mutant devoid of all six pauA genes. The recombinant PauR protein with a hexahistidine tag at its N terminus was purified, and specific bindings of PauR to the promoter regions of most pau operons were demonstrated by electromobility shift assays. Potential interactions of PUT and CAD with PauR were also suggested by chemical cross-linkage analysis with glutaraldehyde. In comparison, growth on PUT was more proficient than that on CAD, and this observed growth phenotype was reflected in a strong catabolite repression of pauA promoter activation by CAD but was completely absent as reflected by activation by PUT. In summary, this study clearly establishes the function of PauR in control of pau promoters in response to PUT and CAD for their catabolism through the γ-glutamylation pathway.

Genome-wide microarrray analysis reveals roles for the REF-1 family member HLH-29 in ferritin synthesis and peroxide stress response.

In Caenorhabditis elegans, the six proteins that make up the REF-1 family have been identified as functional homologs of the Hairy/Enhancer of Split (HES) proteins. These transcription factors act in both Notch dependent and Notch-independent pathways to regulate embryonic events during development; however, their post-embryonic functions are not well defined. As a first step toward understanding how the REF-1 family works together to coordinate post-embryonic events, we used gene expression microarray analysis to identify transcriptional targets of HLH-29 in L4/young adult stage animals. Here we show that HLH-29 targets are genes needed for the regulation of growth and lifespan, including genes required for oxidative stress response and fatty acid metabolism, and the ferritin genes, ftn-1 and ftn-2. We show that HLH-29 regulates ftn-1 expression via promoter sequences upstream of the iron-dependent element that is recognized by the hypoxia inducible factor, HIF-1. Additionally, hlh-29 mutants are more resistant to peroxide stress than wild-type animals and ftn-1(RNAi) animals, even in the presence of excess iron. Finally we show that HLH-29 acts parallel to DAF-16 but upstream of the microphthalmia transcription factor ortholog, HLH-30, to regulate ftn-1 expression under normal growth conditions.

Promoter recognition and activation by the global response regulator CbrB in Pseudomonas aeruginosa.

In Pseudomonas aeruginosa, the CbrA/CbrB two-component system is instrumental in the maintenance of the carbon-nitrogen balance and for growth on carbon sources that are energetically less favorable than the preferred dicarboxylate substrates. The CbrA/CbrB system drives the expression of the small RNA CrcZ, which antagonizes the repressing effects of the catabolite repression control protein Crc, an RNA-binding protein. Dicarboxylates appear to cause carbon catabolite repression by inhibiting the activity of the CbrA/CbrB system, resulting in reduced crcZ expression. Here we have identified a conserved palindromic nucleotide sequence that is present in upstream activating sequences (UASs) of promoters under positive control by CbrB and σ(54) RNA polymerase, especially in the UAS of the crcZ promoter. Evidence for recognition of this palindromic sequence by CbrB was obtained in vivo from mutational analysis of the crcZ promoter and in vitro from electrophoretic mobility shift assays using crcZ promoter fragments and purified CbrB protein truncated at the N terminus. Integration host factor (IHF) was required for crcZ expression. CbrB also activated the lipA (lipase) promoter, albeit less effectively, apparently by interacting with a similar but less conserved palindromic sequence in the UAS of lipA. As expected, succinate caused CbrB-dependent catabolite repression of the lipA promoter. Based on these results and previously published data, a consensus CbrB recognition sequence is proposed. This sequence has similarity to the consensus NtrC recognition sequence, which is relevant for nitrogen control.

L-lysine catabolism is controlled by L-arginine and ArgR in Pseudomonas aeruginosa PAO1.

In comparison to other pseudomonads, Pseudomonas aeruginosa grows poorly in L-lysine as a sole source of nutrient. In this study, the ldcA gene (lysine decarboxylase A; PA1818), previously identified as a member of the ArgR regulon of L-arginine metabolism, was found essential for L-lysine catabolism in this organism. LdcA was purified to homogeneity from a recombinant strain of Escherichia coli, and the results of enzyme characterization revealed that this pyridoxal-5-phosphate-dependent decarboxylase takes L-lysine, but not L-arginine, as a substrate. At an optimal pH of 8.5, cooperative substrate activation by L-lysine was depicted from kinetics studies, with calculated K(m) and V(max) values of 0.73 mM and 2.2 µmole/mg/min, respectively. Contrarily, the ldcA promoter was induced by exogenous L-arginine but not by L-lysine in the wild-type strain PAO1, and the binding of ArgR to this promoter region was demonstrated by electromobility shift assays. This peculiar arginine control on lysine utilization was also noted from uptake experiments in which incorporation of radioactively labeled L-lysine was enhanced in cells grown in the presence of L-arginine but not L-lysine. Rapid growth on L-lysine was detected in a mutant devoid of the main arginine catabolic pathway and with a higher basal level of the intracellular L-arginine pool and hence elevated ArgR-responsive regulons, including ldcA. Growth on L-lysine as a nitrogen source can also be enhanced when the aruH gene encoding an arginine/lysine:pyruvate transaminase was expressed constitutively from plasmids; however, no growth of the ldcA mutant on L-lysine suggests a minor role of this transaminase in L-lysine catabolism. In summary, this study reveals a tight connection of lysine catabolism to the arginine regulatory network, and the lack of lysine-responsive control on lysine uptake and decarboxylation provides an explanation of L-lysine as a poor nutrient for P. aeruginosa.

Inhibition of quorum sensing in Vibrio harveyi by boronic acids.

Bacterial quorum sensing refers to the ability of bacteria to control gene expression through the detection of a threshold concentration of certain chemicals called autoinducer(s), which are secreted by self and/or other bacteria. Quorum sensing is implicated in the regulation of pathologically relevant events such as biofilm formation, virulence, conjugation, sporulation, and swarming mobility. Inhibitors of bacterial quorum sensing are valuable research tools and potential antimicrobial agents. In this paper, we describe the discovery of several boronic acid inhibitors of bacterial quorum sensing in Vibrio harveyi with IC(50) values in the low to sub-micromolar range in whole cell assays.

Synthesis and evaluation of new antagonists of bacterial quorum sensing in Vibrio harveyi.

Bacterial quorum sensing has received much attention in recent years because of its relevance to pathological events such as biofilm formation. Based on the structures of two lead inhibitors (IC50: 35-55 microM) against autoinducer-2-mediated quorum sensing identified through virtual screening, we synthesized 39 analogues and examined their inhibitory activities. Twelve of these new analogues showed equal or better inhibitory activities than the lead inhibitors. The best compound showed an IC50 value of approximately 6 microM in a whole-cell assay using Vibrio harveyi as the model organism. The structure-activity relationship is discussed herein.

Structure-based discovery and experimental verification of novel AI-2 quorum sensing inhibitors against Vibrio harveyi.

Quorum sensing has been implicated in the control of pathologically relevant bacterial behavior such as secretion of virulence factors, biofilm formation, sporulation, and swarming motility. The AI-2 quorum sensing pathway is found in both gram-positive and gram-negative bacteria. Therefore, antagonizing AI-2 quorum sensing is a possible approach to modifying bacterial behaviour. However, efforts in developing inhibitors of AI-2-mediated quorum sensing are especially lacking. High-throughput virtual screening using the V. harveyi LuxP crystal structure identified two compounds that were found to antagonize AI-2-mediated quorum sensing in V. harveyi without cytotoxicity. The sulfone functionality of these inhibitors was identified as critical to their ability to mimic the natural ligand in their interactions with Arg 215 and Arg 310 of the active site.

Identification of boronic acids as antagonists of bacterial quorum sensing in Vibrio harveyi.

Bacterial quorum sensing plays a very important role in the regulation of biofilm formation, virulence, conjugation, sporulation, and swarming mobility. Inhibitors of bacterial quorum sensing are important research tools and potential therapeutic agents. In this paper, we describe for the first time the discovery of several boronic acids as single digit micromolar inhibitors of bacterial quorum sensing in Vibrio harveyi.

Transcriptome analysis of agmatine and putrescine catabolism in Pseudomonas aeruginosa PAO1.

Polyamines (putrescine, spermidine, and spermine) are major organic polycations essential for a wide spectrum of cellular processes. The cells require mechanisms to maintain homeostasis of intracellular polyamines to prevent otherwise severe adverse effects. We performed a detailed transcriptome profile analysis of Pseudomonas aeruginosa in response to agmatine and putrescine with an emphasis in polyamine catabolism. Agmatine serves as the precursor compound for putrescine (and hence spermidine and spermine), which was proposed to convert into 4-aminobutyrate (GABA) and succinate before entering the tricarboxylic acid cycle in support of cell growth, as the sole source of carbon and nitrogen. Two acetylpolyamine amidohydrolases, AphA and AphB, were found to be involved in the conversion of agmatine into putrescine. Enzymatic products of AphA were confirmed by mass spectrometry analysis. Interestingly, the alanine-pyruvate cycle was shown to be indispensable for polyamine utilization. The newly identified dadRAX locus encoding the regulator alanine transaminase and racemase coupled with SpuC, the major putrescine-pyruvate transaminase, were key components to maintaining alanine homeostasis. Corresponding mutant strains were severely hampered in polyamine utilization. On the other hand, an alternative gamma-glutamylation pathway for the conversion of putrescine into GABA is present in some organisms. Subsequently, GabD, GabT, and PA5313 were identified for GABA utilization. The growth defect of the PA5313 gabT double mutant in GABA suggested the importance of these two transaminases. The succinic-semialdehyde dehydrogenase activity of GabD and its induction by GABA were also demonstrated in vitro. Polyamine utilization in general was proven to be independent of the PhoPQ two-component system, even though a modest induction of this operon was induced by polyamines. Multiple potent catabolic pathways, as depicted in this study, could serve pivotal roles in the control of intracellular polyamine levels.

CDK-dependent activation of poly(ADP-ribose) polymerase member 10 (PARP10).

Proteins of the poly(ADP-ribose) polymerase (PARP) family play a wide array of functions, covering virtually every aspect of DNA metabolism and function, most notably with the response to DNA damage, transcription, and the maintenance of genomic stability. Here we report the identification and characterization of a novel PARP family member, PARP10 (FLJ14464 or hypothetical protein LOC84875). Overexpression of PARP10 results in loss of cell viability, although down-expression by short hairpin RNA leads to delayed G1 progression and concomitant cell death. PARP10 exists in both cytoplasm and nucleus, but only nucleolar PARP10 acquires CDK-dependent phosphorylation through late-G1 to S phase, and from prometaphase to cytokinesis in the nucleolar organizing regions. The PARP activity of PARP10 depends on phosphorylation by CDK2-cyclin E in vitro. CDK-phosphorylated PARP10 is absent in growth-arrested cells. These results suggest that PARP10 functions in cell proliferation and may serve as a marker for proliferating cells.