Operons

Genome-Wide Identification of Dickeya solani Transcriptional Units Up-Regulated in Response to Plant Tissues From a Crop-Host Solanum tuberosum and a Weed-Host Solanum dulcamara

Dickeya solani is a Gram-negative bacterium able to cause disease symptoms on a variety of crop and ornamental plants worldwide. Weeds including Solanum dulcamara (bittersweet nightshade) growing near agricultural fields have been reported to support populations of soft rot bacteria in natural settings. However, little is known about the specific interaction of D. solani with such weed plants that may contribute to its success as an agricultural pathogen. The aim of this work was to assess the interaction of D. solani with its crop plant (Solanum tuberosum) and an alternative (S. dulcamara) host plant. From a collection of 10,000 Tn5 transposon mutants of D. solani IPO2222 carrying an inducible, promotorless gusA reporter gene, 210 were identified that exhibited plant tissue-dependent expression of the gene/operon into which the Tn5 insertion had occurred. Thirteen Tn5 mutants exhibiting the greatest plant tissue induction of such transcriptional units in S. tuberosum or S. dulcamara as measured by qRT-PCR were assessed for plant host colonization, virulence, and ability to macerate plant tissue, as well as phenotypes likely to contribute to the ecological fitness of D. solani, including growth rate, carbon and nitrogen source utilization, motility, chemotaxis toward plant extracts, biofilm formation, growth under anaerobic conditions and quorum sensing. These 13 transcriptional units encode proteins involved in bacterial interactions with plants, with functions linked to cell envelope structure, chemotaxis and carbon metabolism. The selected 13 genes/operons were differentially expressed in, and thus contributed preferentially to D. solani fitness in potato and/or S. dulcamara stem, leaf, and root tissues.

Phylogenetic comparison between Type IX Secretion System (T9SS) protein components suggests evidence of horizontal gene transfer

Porphyromonas gingivalis is one of the major bacteria that causes periodontitis. Chronic periodontitis is a severe form of periodontal disease that ultimately leads to tooth loss. Virulence factors that contribute to periodontitis are secreted by Type IX Secretion System (T9SS). There are aspects of T9SS protein components that have yet to be characterised. Thus, the aim of this study is to investigate the phylogenetic relationship between members of 20 T9SS component protein families. The Bayesian Inference (BI) trees for 19 T9SS protein components exhibit monophyletic clades for all major classes under Bacteroidetes with strong support for the monophyletic clades or its subclades that is consistent with phylogeny exhibited by the constructed BI tree of 16S rRNA. The BI tree of PorR is different from the 19 BI trees of T9SS protein components as it does not exhibit monophyletic clades for all major classes under Bacteroidetes. There is strong support for the phylogeny exhibited by the BI tree of PorR which deviates from the phylogeny based on 16S rRNA. Hence, it is possible that the porR gene is subjected to horizontal transfer as it is known that virulence factor genes could be horizontally transferred. Seven genes (porR included) that are involved in the biosynthesis of A-LPS are found to be flanked by insertion sequences (IS5 family transposons). Therefore, the intervening DNA segment that contains the porR gene might be transposed and subjected to conjugative transfer. Thus, the seven genes can be co-transferred via horizontal gene transfer. The BI tree of UgdA does not exhibit monophyletic clades for all major classes under Bacteroidetes which is similar to the BI tree of PorR (both are a part of the seven genes). Both BI trees also exhibit similar topology as the four identified clusters with strong support and have similar relative positions to each other in both BI trees. This reinforces the possibility that porR and the other six genes might be horizontally transferred. Other than the BI tree of PorR, the 19 other BI trees of T9SS protein components also exhibit evidence of horizontal gene transfer. However, their genes might undergo horizontal gene transfer less frequently compared to porR because the intervening DNA segment that contains porR is easily exchanged between bacteria under Bacteroidetes due to the presence of insertion sequences (IS5 family transposons) that flank it. In conclusion, this study can provide a better understanding about the phylogeny of T9SS protein components.

Active and repressed biosynthetic gene clusters have spatially distinct chromosome states

While colocalization within a bacterial operon enables coexpression of the constituent genes, the mechanistic logic of clustering of nonhomologous monocistronic genes in eukaryotes is not immediately obvious. Biosynthetic gene clusters that encode pathways for specialized metabolites are an exception to the classical eukaryote rule of random gene location and provide paradigmatic exemplars with which to understand eukaryotic cluster dynamics and regulation. Here, using 3C, Hi-C, and Capture Hi-C (CHi-C) organ-specific chromosome conformation capture techniques along with high-resolution microscopy, we investigate how chromosome topology relates to transcriptional activity of clustered biosynthetic pathway genes in Arabidopsis thaliana Our analyses reveal that biosynthetic gene clusters are embedded in local hot spots of 3D contacts that segregate cluster regions from the surrounding chromosome environment. The spatial conformation of these cluster-associated domains differs between transcriptionally active and silenced clusters. We further show that silenced clusters associate with heterochromatic chromosomal domains toward the periphery of the nucleus, while transcriptionally active clusters relocate away from the nuclear periphery. Examination of chromosome structure at unrelated clusters in maize, rice, and tomato indicates that integration of clustered pathway genes into distinct topological domains is a common feature in plant genomes. Our results shed light on the potential mechanisms that constrain coexpression within clusters of nonhomologous eukaryotic genes and suggest that gene clustering in the one-dimensional chromosome is accompanied by compartmentalization of the 3D chromosome.

Targeted Inhibition of Gut Microbial Trimethylamine N-Oxide Production Reduces Renal Tubulointerstitial Fibrosis and Functional Impairment in a Murine Model of Chronic Kidney Disease

OBJECTIVE

Gut microbial metabolism of dietary choline, a nutrient abundant in a Western diet, produces trimethylamine (TMA) and the atherothrombosis- and fibrosis-promoting metabolite TMA-N-oxide (TMAO). Recent clinical and animal studies reveal that elevated TMAO levels are associated with heightened risks for both cardiovascular disease and incident chronic kidney disease development. Despite this, studies focusing on therapeutically targeting gut microbiota-dependent TMAO production and its impact on preserving renal function are limited. Approach and Results: Herein we examined the impact of pharmacological inhibition of choline diet-induced gut microbiota-dependent production of TMA, and consequently TMAO, on renal tubulointerstitial fibrosis and functional impairment in a model of chronic kidney disease. Initial studies with a gut microbial choline TMA-lyase mechanism-based inhibitor, iodomethylcholine, confirmed both marked suppression of TMA generation, and consequently TMAO levels, and selective targeting of the gut microbial compartment (ie, both accumulation of the drug in intestinal microbes and limited systemic exposure in the host). Dietary supplementation of either choline or TMAO significantly augmented multiple indices of renal functional impairment and fibrosis associated with chronic subcutaneous infusion of isoproterenol. However, the presence of the gut microbiota-targeting inhibitor iodomethylcholine blocked choline diet-induced elevation in TMAO, and both significantly improved decline in renal function, and significantly attenuated multiple indices of tubulointerstitial fibrosis. Iodomethylcholine treatment also reversed many choline diet-induced changes in cecal microbial community composition associated with TMAO and renal functional impairment.

CONCLUSIONS

Selective targeting of gut microbiota-dependent TMAO generation may prevent adverse renal structural and functional alterations in subjects at risk for chronic kidney disease.

Evolutionary and functional patterns of shared gene neighbourhood in fungi

Gene clusters comprise genomically co-localized and potentially co-regulated genes that tend to be conserved across species. In eukaryotes, multiple examples of metabolic gene clusters are known, particularly among fungi and plants. However, little is known about how gene clustering patterns vary among taxa or with respect to functional roles. Furthermore, mechanisms of the formation, maintenance and evolution of gene clusters remain unknown. We surveyed 341 fungal genomes to discover gene clusters shared by different species, independently of their functions. We inferred 12,120 cluster families, which comprised roughly one third of the gene space and were enriched in genes associated with diverse cellular functions. Additionally, most clusters did not encode transcription factors, suggesting that they are regulated distally. We used phylogenomics to characterize the evolutionary history of these clusters. We found that most clusters originated once and were transmitted vertically, coupled to differential loss. However, convergent evolution-that is, independent appearance of the same cluster-was more prevalent than anticipated. Finally, horizontal gene transfer of entire clusters was somewhat restricted, with the exception of those associated with secondary metabolism. Altogether, our results provide insights on the evolution of gene clustering as well as a broad catalogue of evolutionarily conserved gene clusters whose function remains to be elucidated.

Unexpected diversity in eukaryotic transcription revealed by the retrotransposon hotspot family of Trypanosoma brucei

The path from DNA to RNA to protein in eukaryotes is guided by a series of factors linking transcription, mRNA export and translation. Many of these are conserved from yeast to humans. Trypanosomatids, which diverged early in the eukaryotic lineage, exhibit unusual features such as polycistronic transcription and trans-splicing of all messenger RNAs. They possess basal transcription factors, but lack recognisable orthologues of many factors required for transcription elongation and mRNA export. We show that retrotransposon hotspot (RHS) proteins fulfil some of these functions and that their depletion globally impairs nascent RNA synthesis by RNA polymerase II. Three sub-families are part of a coordinated process in which RHS6 is most closely associated with chromatin, RHS4 is part of the Pol II complex and RHS2 connects transcription with the translation machinery. In summary, our results show that the components of eukaryotic transcription are far from being universal, and reveal unsuspected plasticity in the course of evolution.

Start-Stop Assembly: a functionally scarless DNA assembly system optimized for metabolic engineering

DNA assembly allows individual DNA constructs or libraries to be assembled quickly and reliably. Most methods are either: (i) Modular, easily scalable and suitable for combinatorial assembly, but leave undesirable 'scar' sequences; or (ii) bespoke (non-modular), scarless but less suitable for construction of combinatorial libraries. Both have limitations for metabolic engineering. To overcome this trade-off we devised Start-Stop Assembly, a multi-part, modular DNA assembly method which is both functionally scarless and suitable for combinatorial assembly. Crucially, 3 bp overhangs corresponding to start and stop codons are used to assemble coding sequences into expression units, avoiding scars at sensitive coding sequence boundaries. Building on this concept, a complete DNA assembly framework was designed and implemented, allowing assembly of up to 15 genes from up to 60 parts (or mixtures); monocistronic, operon-based or hybrid configurations; and a new streamlined assembly hierarchy minimizing the number of vectors. Only one destination vector is required per organism, reflecting our optimization of the system for metabolic engineering in diverse organisms. Metabolic engineering using Start-Stop Assembly was demonstrated by combinatorial assembly of carotenoid pathways in Escherichia coli resulting in a wide range of carotenoid production and colony size phenotypes indicating the intended exploration of design space.

A Two-Component System That Modulates Cyclic di-GMP Metabolism Promotes Legionella pneumophila Differentiation and Viability in Low-Nutrient Conditions

During its life cycle, the environmental pathogen Legionella pneumophila alternates between a replicative and transmissive cell type when cultured in broth, macrophages, or amoebae. Within a protozoan host, L. pneumophila further differentiates into the hardy cell type known as the mature infectious form (MIF). The second messenger cyclic di-GMP coordinates lifestyle changes in many bacterial species, but its role in the L. pneumophila life cycle is less understood. Using an in vitro broth culture model that approximates the intracellular transition from the replicative to the transmissive form, here we investigate the contribution to L. pneumophila differentiation of a two-component system (TCS) that regulates cyclic di-GMP metabolism. The TCS is encoded by lpg0278-lpg0277 and is cotranscribed with lpg0279, which encodes a protein upregulated in MIF cells. The promoter for this operon is RpoS dependent and induced in nutrient-limiting conditions that do not support replication, as demonstrated using a gfp reporter and quantitative PCR (qPCR). The response regulator of the TCS (Lpg0277) is a bifunctional enzyme that both synthesizes and degrades cyclic di-GMP. Using a panel of site-directed point mutants, we show that cyclic di-GMP synthesis mediated by a conserved GGDEF domain promotes growth arrest of replicative L. pneumophila, accumulation of pigment and poly-3-hydroxybutyrate storage granules, and viability in nutrient-limiting conditions. Genetic epistasis tests predict that the MIF protein Lpg0279 acts as a negative regulator of the TCS. Thus, L. pneumophila is equipped with a regulatory network in which cyclic di-GMP stimulates the switch from a replicative to a resilient state equipped to survive in low-nutrient environments.IMPORTANCE Although an intracellular pathogen, L. pneumophila has developed mechanisms to ensure long-term survival in low-nutrient aqueous conditions. Eradication of L. pneumophila from contaminated water supplies has proven challenging, as outbreaks have been traced to previously remediated systems. Understanding the genetic determinants that support L. pneumophila persistence in low-nutrient environments can inform design and assessment of remediation strategies. Here we characterize a genetic locus that encodes a two-component signaling system (lpg0278-lpg0277) and a putative regulator protein (lpg0279) that modulates the production of the messenger molecule cyclic di-GMP. We show that this locus promotes both L. pneumophila cell differentiation and survival in nutrient-limiting conditions, thus advancing the understanding of the mechanisms that contribute to L. pneumophila environmental resilience.

Stronger When Together: Clustering of Plant NLR Disease resistance Genes

Gene clustering is rare in eukaryotes. However, nucleotide-binding leucine-rich repeat (NLR)-encoding disease resistance (R) genes show consistent clustering in plant genomes. These arrangements are likely to provide coregulatory benefits, as suggested by growing evidence that the gene products of both paired and larger clusters of NLRs act together in triggering immunity. Head-to-head gene pairs where one of the encoded NLRs includes an integrated decoy domain appear to behave differently than clusters evolved from closely related typical NLRs. These patterns may help to explain the broad resistance that most plants have despite their finite number of R genes. By taking into consideration the relationship between genomic arrangement and function, we can improve our understanding of and ability to predict plant immune detection.

Multifaceted Stoichiometry Control of Bacterial Operons Revealed by Deep Proteome Quantification

More than half of the protein-coding genes in bacteria are organized in polycistronic operons composed of two or more genes. It remains under debate whether the operon organization maintains the stoichiometric expression of the genes within an operon. In this study, we performed a label-free data-independent acquisition hyper reaction monitoring mass-spectrometry (HRM-MS) experiment to quantify the Escherichia coli proteome in exponential phase and quantified 93.6% of the cytosolic proteins, covering 67.9% and 56.0% of the translating polycistronic operons in BW25113 and MG1655 strains, respectively. We found that the translational regulation contributes largely to the proteome complexity: the shorter operons tend to be more tightly controlled for stoichiometry than longer operons; the operons which mainly code for complexes is more tightly controlled for stoichiometry than the operons which mainly code for metabolic pathways. The gene interval (distance between adjacent genes in one operon) may serve as a regulatory factor for stoichiometry. The catalytic efficiency might be a driving force for differential expression of enzymes encoded in one operon. These results illustrated the multifaceted nature of the operon regulation: the operon unified transcriptional level and gene-specific translational level. This multi-level regulation benefits the host by optimizing the efficiency of the productivity of metabolic pathways and maintenance of different types of protein complexes.

Suppressing photorespiration for the improvement in photosynthesis and crop yields: A review on the role of S-allantoin as a nitrogen source

Environmental variations resulting in biotic and abiotic stresses demand adaptive changes in the photosynthetic machinery. To cope with these challenges, plant scientists are constantly striving to enhance photosynthetic activity. The photorespiration pathway, which fixes O₂ and releases CO₂ in C₃ plants, competes with photosynthesis. One method to increase yield would be to enhance photosynthesis by engineering the photorespiratory pathway. To date, three engineered photorespiratory pathways have been produced, of which two have been proven experimentally in the model plant, Arabidopsis thaliana. These approaches might be helpful in enhancing crop resilience to future environmental challenges. In partially photorespiratory suppressed plants, it is hypothesized that a gene cluster may have formed between bacterial glycolate dehydrogenase (GDH), glyoxylate carboligase (GCL), and tartronic semi aldehyde (TSR) genes with Arabidopsis allantoin degradation genes like Arabidopsis allantoinase (AtALN) to utilize S-allantoin as a source of nitrogen. Observations of the use of allantoin as an exclusive source of nitrogen or energy by Arabidopsis and Escherichia coli led us to propose a genetic switch control model between nitrogen assimilation and energy producing pathways in partially photorespiratory suppressed plants.

Recombination drives the evolution of mutational robustness

Recombination can impose fitness costs as beneficial parental combinations of alleles are broken apart, a phenomenon known as recombination load. Computational models suggest that populations may evolve a reduced recombination load by reducing either the likelihood of recombination events (bring interacting loci in physical proximity) or the strength of interactions between loci (make loci more independent of one another). We review evidence for each of these possibilities and their consequences for the genotype-fitness relationship. In particular, we expect that reducing interaction strengths between loci will lead to genomes that are also robust to mutational perturbations, but reducing recombination rates alone will not. We note that both mechanisms most likely played a role in the evolution of extant populations, and that both can result in the frequently-observed pattern of physical linkage between interacting loci.

A genomics perspective on natural product biosynthesis in plant pathogenic bacteria

This review summarizes findings from genomics-inspired natural product research in plant pathogenic bacteria and discusses emerging trends in this field.

Image-based promoter prediction: a promoter prediction method based on evolutionarily generated patterns

Prediction of promoter regions is crucial for studying gene function and regulation. The well-accepted position weight matrix method for this purpose relies on predefined motifs, which would hinder application across different species. Here, we introduce image-based promoter prediction (IBPP) as a method that creates an “image” from training promoter sequences using an evolutionary approach and predicts promoters by matching with the “image”. We used Escherichia coli σ70 promoter sequences to test the performance of IBPP and the combination of IBPP and a support vector machine algorithm (IBPP-SVM). The “images” generated with IBPP could effectively distinguish promoter from non-promoter sequences. Compared with IBPP, IBPP-SVM showed a substantial improvement in sensitivity. Furthermore, both methods showed good performance for sequences of up to 2,000 nt in length. The performances of IBPP and IBPP-SVM were largely affected by the threshold and dimension of vectors, respectively. The source code and documentation are freely available at https://github.com/hahatcdg/IBPP.

Development of a gut microbe-targeted nonlethal therapeutic to inhibit thrombosis potential

Trimethylamine N-oxide (TMAO) is a gut microbiota-derived metabolite that enhances both platelet responsiveness and in vivo thrombosis potential in animal models, and TMAO plasma levels predict incident atherothrombotic event risks in human clinical studies. TMAO is formed by gut microbe-dependent metabolism of trimethylamine (TMA) moiety-containing nutrients, which are abundant in a Western diet. Here, using a mechanism-based inhibitor approach targeting a major microbial TMA-generating enzyme pair, CutC and CutD (CutC/D), we developed inhibitors that are potent, time-dependent, and irreversible and that do not affect commensal viability. In animal models, a single oral dose of a CutC/D inhibitor significantly reduced plasma TMAO levels for up to 3 d and rescued diet-induced enhanced platelet responsiveness and thrombus formation, without observable toxicity or increased bleeding risk. The inhibitor selectively accumulated within intestinal microbes to millimolar levels, a concentration over 1-million-fold higher than needed for a therapeutic effect. These studies reveal that mechanism-based inhibition of gut microbial TMA and TMAO production reduces thrombosis potential, a critical adverse complication in heart disease. They also offer a generalizable approach for the selective nonlethal targeting of gut microbial enzymes linked to host disease limiting systemic exposure of the inhibitor in the host.

Investigation on the Evolutionary Relation of Diverse Polyhydroxyalkanoate Gene Clusters in Betaproteobacteria

Products of numerous genes (phaC, phaA, phaB, phaP, phaR, and phaZ) are involved in the synthesis and degradation processes of the ubiquitous prokaryotic polyhydroxyalkanoate (PHA) intracellular reserve storage system. In this study, we performed a bioinformatics analysis to identify PHA-related genes and proteins in the genome of 66 selected organisms (class: Betaproteobacteria) that occur in various habitats; besides, evolutionary trajectories of the PHA system are reported here. The identified PHA-related genes were organized into clusters, and the gene arrangement was highly diverse. The occurrence and distribution of PHA-related clusters revealed that a single cluster was primarily segmented into small gene groups among various genomes, which were further reorganized as novel clusters based on various functional genes. The individual phylogenies of gene and protein sequences supported that the clusters were assembled through the relocation of native orthologous genes that underwent insertion, deletion, and elongation events. Furthermore, the neighboring genes provided valuable evolutionary and functional cues regarding the conservation and maintenance of PHA-related genes in the genome. Overall, the aforementioned results strongly indicate the influence of horizontal gene transfer on the organization of PHA-related gene clusters. Therefore, our results reveal new insights into the organization, evolutionary history, and cluster conservation of the PHA-related gene inventories among Betaproteobacterial organisms.

Computational operon prediction in whole-genomes and metagenomes

Microbial diversity in unique environmental settings enables abrupt responses catalysed by altering the gene regulation and formation of gene clusters called operons. Operons increases bacterial adaptability, which in turn increases their survival. This review article presents the emergence of computational operon prediction methods for whole microbial genomes and metagenomes, and discusses their strengths and limitations. Most of the whole-genome operon prediction methods struggle to generalize on unrelated genomes. The applicability of universal whole-genome operon prediction methods to metagenomic data is an interesting yet less investigated question. We have evaluated the potential of various operon prediction features for genomic and metagenomic data. Most of operon prediction methods with high accuracy have been compiled into databases. Despite of the high predictive performance, the data among many databases are not completely consistent for similar species. We performed a correlation analysis between the computationally predicted operon databases and experimentally validated data for Escherichia coli, Bacillus subtilis and Mycobacterium tuberculosis. Operon prediction for most of the less characterized microbes cannot be verified due to absence of experimentally validated operons. The generation of validated information for other microbes would test the authenticity of operon databases for other less annotated microbes as well. Advances in sequencing technologies and development of better analysis methods will help researchers to overcome the technological hurdles (such as long sequencing reads and improved contig size) and further improve operon predictions and better utilize operonic information.

Phylogenetic diversification and developmental implications of poly-(R)-3-hydroxyalkanoate gene cluster assembly in prokaryotes

Many polyhydroxyalkanoates (PHAs) system genes, such as phaC, phaA, phaB, phaR, phaP and phaZ, are often found to be organised in the form of operon-like clusters. In this study, a genome survey was performed to identify such clustered PHA systems among 256 prokaryotic organisms. These data were then used to generate a comprehensive 16S rRNA species tree depicting the phylogenetic distribution of the observed clusters with diverse gene arrangements. In addition, the gene occurrences and physical linkages between PHA system genes were quantitatively estimated. From this, we identified a centrally connected hub gene, i.e. the phaC gene of PHA. Furthermore, a comparative investigation was performed between the clusters of PHA and glycogen, which decoded the role of the hub gene in the cluster organisation of both systems. Together, these findings suggest that the highly connected hub gene might contribute substantively towards the organisation and maintenance of the gene network connectivity in the clusters, particularly in the storage reserve systems.

Mapping Causal Variants with Single-Nucleotide Resolution Reveals Biochemical Drivers of Phenotypic Change

Understanding the sequence determinants that give rise to diversity among individuals and species is the central challenge of genetics. However, despite ever greater numbers of sequenced genomes, most genome-wide association studies cannot distinguish causal variants from linked passenger mutations spanning many genes. We report that this inherent challenge can be overcome in model organisms. By pushing the advantages of inbred crossing to its practical limit in Saccharomyces cerevisiae, we improved the statistical resolution of linkage analysis to single nucleotides. This "super-resolution" approach allowed us to map 370 causal variants across 26 quantitative traits. Missense, synonymous, and cis-regulatory mutations collectively gave rise to phenotypic diversity, providing mechanistic insight into the basis of evolutionary divergence. Our data also systematically unmasked complex genetic architectures, revealing that multiple closely linked driver mutations frequently act on the same quantitative trait. Single-nucleotide mapping thus complements traditional deletion and overexpression screening paradigms and opens new frontiers in quantitative genetics.

RNAseq analysis of α-proteobacterium Gluconobacter oxydans 621H

Background

The acetic acid bacterium Gluconobacter oxydans 621H is characterized by its exceptional ability to incompletely oxidize a great variety of carbohydrates in the periplasm. The metabolism of this α-proteobacterium has been characterized to some extent, yet little is known about its transcriptomes and related data. In this study, we applied two different RNAseq approaches. Primary transcriptomes enriched for 5′-ends of transcripts were sequenced to detect transcription start sites, which allow subsequent analysis of promoter motifs, ribosome binding sites, and 5´-UTRs. Whole transcriptomes were sequenced to identify expressed genes and operon structures.

Results

Sequencing of primary transcriptomes of G. oxydans revealed 2449 TSSs, which were classified according to their genomic context followed by identification of promoter and ribosome binding site motifs, analysis of 5´-UTRs including validation of predicted cis-regulatory elements and correction of start codons. 1144 (41%) of all genes were found to be expressed monocistronically, whereas 1634 genes were organized in 571 operons. Together, TSSs and whole transcriptome data were also used to identify novel intergenic (18), intragenic (328), and antisense transcripts (313).

Conclusions

This study provides deep insights into the transcriptional landscapes of G. oxydans. The comprehensive transcriptome data, which we made publicly available, facilitate further analysis of promoters and other regulatory elements. This will support future approaches for rational strain development and targeted gene expression in G. oxydans. The corrections of start codons further improve the high quality genome reference and support future proteome analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4415-x) contains supplementary material, which is available to authorized users.

Structural characterization of the bacterial proteasome homolog BPH reveals a tetradecameric double-ring complex with unique inner cavity properties

Eukaryotic and archaeal proteasomes are paradigms for self-compartmentalizing proteases. To a large extent, their function requires interplay with hexameric ATPases associated with diverse cellular activities (AAA+) that act as substrate unfoldases. Bacteria have various types of self-compartmentalizing proteases; in addition to the proteasome itself, these include the proteasome homolog HslV, which functions together with the AAA+ HslU; the ClpP protease with its partner AAA+ ClpX; and Anbu, a recently characterized ancestral proteasome variant. Previous bioinformatic analysis has revealed a novel bacterial member of the proteasome family Betaproteobacteria proteasome homolog (BPH). Using cluster analysis, we here affirmed that BPH evolutionarily descends from HslV. Crystal structures of the Thiobacillus denitrificans and Cupriavidus metallidurans BPHs disclosed a homo-oligomeric double-ring architecture in which the active sites face the interior of the cylinder. Using small-angle X-ray scattering (SAXS) and electron microscopy averaging, we found that BPH forms tetradecamers in solution, unlike the dodecamers seen in HslV. Although the highly acidic inner surface of BPH was in striking contrast to the cavity characteristics of the proteasome and HslV, a classical proteasomal reaction mechanism could be inferred from the covalent binding of the proteasome-specific inhibitor epoxomicin to BPH. A ligand-bound structure implied that the elongated BPH inner pore loop may be involved in substrate recognition. The apparent lack of a partner unfoldase and other unique features, such as Ser replacing Thr as the catalytic residue in certain BPH subfamilies, suggest a proteolytic function for BPH distinct from those of known bacterial self-compartmentalizing proteases.

Genome-wide identification of physically clustered genes suggests chromatin-level co-regulation in male reproductive development in Arabidopsis thaliana

Co-expression of physically linked genes occurs surprisingly frequently in eukaryotes. Such chromosomal clustering may confer a selective advantage as it enables coordinated gene regulation at the chromatin level. We studied the chromosomal organization of genes involved in male reproductive development in Arabidopsis thaliana. We developed an in-silico tool to identify physical clusters of co-regulated genes from gene expression data. We identified 17 clusters (96 genes) involved in stamen development and acting downstream of the transcriptional activator MS1 (MALE STERILITY 1), which contains a PHD domain associated with chromatin re-organization. The clusters exhibited little gene homology or promoter element similarity, and largely overlapped with reported repressive histone marks. Experiments on a subset of the clusters suggested a link between expression activation and chromatin conformation: qRT-PCR and mRNA in situ hybridization showed that the clustered genes were up-regulated within 48 h after MS1 induction; out of 14 chromatin-remodeling mutants studied, expression of clustered genes was consistently down-regulated only in hta9/hta11, previously associated with metabolic cluster activation; DNA fluorescence in situ hybridization confirmed that transcriptional activation of the clustered genes was correlated with open chromatin conformation. Stamen development thus appears to involve transcriptional activation of physically clustered genes through chromatin de-condensation.

DprA from Neisseria meningitidis: properties and role in natural competence for transformation

DNA processing chain A (DprA) is a DNA-binding protein that is ubiquitous in bacteria and expressed in some archaea. DprA is active in many bacterial species that are competent for transformation of DNA, but its role in Neisseriameningitidis (Nm) is not well characterized. An Nm mutant lacking DprA was constructed, and the phenotypes of the wild-type and ΔdprA mutant were compared. The salient feature of the phenotype of dprA null cells is the total lack of competence for genetic transformation shown by all of the donor DNA substrates tested in this study. Here, Nm wild-type and dprA null cells appeared to be equally resistant to genotoxic stress. The gene encoding DprA_Nm was cloned and overexpressed, and the biological activities of DprA_Nm were further investigated. DprA_Nm binds ssDNA more strongly than dsDNA, but lacks DNA uptake sequence-specific DNA binding. DprA_Nm dimerization and interaction with the C-terminal part of the single-stranded binding protein SSB_Nmwere demonstrated. dprA is co-expressed with smg, a downstream gene of unknown function, and the gene encoding topoisomerase 1, topA.

BsmR degrades c-di-GMP to modulate biofilm formation of nosocomial pathogen Stenotrophomonas maltophilia

c-di-GMP is a cellular second messenger that regulates diverse bacterial processes, including swimming, biofilm formation and virulence. However, in Stenotrophomonas maltophilia, a nosocomial pathogen that frequently infects immunodeficient or immunoincompetent patients, the regulatory function of c-di-GMP remains unclear. Here we show that BsmR is a negative regulator of biofilm development that degrades c-di-GMP through its EAL domain. Increasing BsmR expression resulted in significant increase in bacterial swimming and decrease in cell aggregation. BsmR regulates the expression of at least 349 genes. Among them, 34 involved in flagellar assembly and a flagellar-assembly-related transcription factor (fsnR) are positively regulated. Although BsmR is a response regulator of the two-component signaling system, its role in biofilm formation depends on the expression level of its respective gene (bsmR), not on the protein’s phosphorylation level. A transcription factor, BsmT, whose coding gene is located in the same tetra-cistronic operon as bsmR, was shown to directly bind to the promoter region of the operon and, through a positive regulatory loop, modulate bsmR transcription. Thus, our results revealed that the c-di-GMP signaling pathway controls biofilm formation and swimming in S. maltophilia, suggesting c-di-GMP signaling as a target in the development of novel antibacterial agents to resist this pathogen.

Plant metabolic clusters - from genetics to genomics

Contents 771 I. 771 II. 772 III. 780 IV. 781 V. 786 786 References 786 SUMMARY: Plant natural products are of great value for agriculture, medicine and a wide range of other industrial applications. The discovery of new plant natural product pathways is currently being revolutionized by two key developments. First, breakthroughs in sequencing technology and reduced cost of sequencing are accelerating the ability to find enzymes and pathways for the biosynthesis of new natural products by identifying the underlying genes. Second, there are now multiple examples in which the genes encoding certain natural product pathways have been found to be grouped together in biosynthetic gene clusters within plant genomes. These advances are now making it possible to develop strategies for systematically mining multiple plant genomes for the discovery of new enzymes, pathways and chemistries. Increased knowledge of the features of plant metabolic gene clusters - architecture, regulation and assembly - will be instrumental in expediting natural product discovery. This review summarizes progress in this area.

Evolutionary trajectory of phytoalexin biosynthetic gene clusters in rice

Plants frequently possess operon-like gene clusters for specialized metabolism. Cultivated rice, Oryza sativa, produces antimicrobial diterpene phytoalexins represented by phytocassanes and momilactones, and the majority of their biosynthetic genes are clustered on chromosomes 2 and 4, respectively. These labdane-related diterpene phytoalexins are biosynthesized from geranylgeranyl diphosphate via ent-copalyl diphosphate or syn-copalyl diphosphate. The two gene clusters consist of genes encoding diterpene synthases and chemical-modification enzymes including P450s. In contrast, genes for the biosynthesis of gibberellins, which are labdane-related phytohormones, are scattered throughout the rice genome similar to other plant genomes. The mechanism of operon-like gene cluster formation remains undefined despite previous studies in other plant species. Here we show an evolutionary insight into the rice gene clusters by a comparison with wild Oryza species. Comparative genomics and biochemical studies using wild rice species from the AA genome lineage, including Oryza barthii, Oryza glumaepatula, Oryza meridionalis and the progenitor of Asian cultivated rice Oryza rufipogon indicate that gene clustering for biosynthesis of momilactones and phytocassanes had already been accomplished before the domestication of rice. Similar studies using the species Oryza punctata from the BB genome lineage, the distant FF genome lineage species Oryza brachyantha and an outgroup species Leersia perrieri suggest that the phytocassane biosynthetic gene cluster was present in the common ancestor of the Oryza species despite the different locations, directions and numbers of their member genes. However, the momilactone biosynthetic gene cluster evolved within Oryza before the divergence of the BB genome via assembly of ancestral genes.

The Role of Transcription Factors at Antisense-Expressing Gene Pairs in Yeast

Genes encoded close to one another on the chromosome are often coexpressed, by a mechanism and regulatory logic that remain poorly understood. We surveyed the yeast genome for tandem gene pairs oriented tail-to-head at which expression antisense to the upstream gene was conserved across species. The intergenic region at most such tandem pairs is a bidirectional promoter, shared by the downstream gene mRNA and the upstream antisense transcript. Genomic analyses of these intergenic loci revealed distinctive patterns of transcription factor regulation. Mutation of a given transcription factor verified its role as a regulator in trans of tandem gene pair loci, including the proximally initiating upstream antisense transcript and downstream mRNA and the distally initiating upstream mRNA. To investigate cis-regulatory activity at such a locus, we focused on the stress-induced NAD(P)H dehydratase YKL151C and its downstream neighbor, the metabolic enzyme GPM1. Previous work has implicated the region between these genes in regulation of GPM1 expression; our mutation experiments established its function in rich medium as a repressor in cis of the distally initiating YKL151C sense RNA, and an activator of the proximally initiating YKL151C antisense RNA. Wild-type expression of all three transcripts required the transcription factor Gcr2. Thus, at this locus, the intergenic region serves as a focal point of regulatory input, driving antisense expression and mediating the coordinated regulation of YKL151C and GPM1. Together, our findings implicate transcription factors in the joint control of neighboring genes specialized to opposing conditions and the antisense transcripts expressed between them.

A complete annotation of the chromosomes of the cellulase producer Trichoderma reesei provides insights in gene clusters, their expression and reveals genes required for fitness

BACKGROUND

Investigations on a few eukaryotic model organisms showed that many genes are non-randomly distributed on chromosomes. In addition, chromosome ends frequently possess genes that are important for the fitness of the organisms. Trichoderma reesei is an industrial producer of enzymes for food, feed and biorefinery production. Its seven chromosomes have recently been assembled, thus making an investigation of its chromosome architecture possible.

RESULTS

We manually annotated and mapped 9194 ORFs on their respective chromosomes and investigated the clustering of the major gene categories and of genes encoding carbohydrate-active enzymes (CAZymes), and the relationship between clustering and expression. Genes responsible for RNA processing and modification, amino acid metabolism, transcription, translation and ribosomal structure and biogenesis indeed showed loose clustering, but this had no impact on their expression. A third of the genes encoding CAZymes also occurred in loose clusters that also contained a high number of genes encoding small secreted cysteine-rich proteins. Five CAZyme clusters were located less than 50 kb apart from the chromosome ends. These genes exhibited the lowest basal (but not induced) expression level, which correlated with an enrichment of H3K9 methylation in the terminal 50 kb areas indicating gene silencing. No differences were found in the expression of CAZyme genes present in other parts of the chromosomes. The putative subtelomeric areas were also enriched in genes encoding secreted proteases, amino acid permeases, enzyme clusters for polyketide synthases (PKS)-non-ribosomal peptide synthase (NRPS) fusion proteins (PKS-NRPS) and proteins involved in iron scavenging. They were strongly upregulated during conidiation and interaction with other fungi.

CONCLUSIONS

Our findings suggest that gene clustering on the T. reesei chromosomes occurs but generally has no impact on their expression. CAZyme genes, located in subtelomers, however, exhibited a much lower basal expression level. The gene inventory of the subtelomers suggests a major role of competition for nitrogen and iron supported by antibiosis for the fitness of T. reesei. The availability of fully annotated chromosomes will facilitate the use of genetic crossings in identifying still unknown genes responsible for specific traits of T. reesei.

Gene selection and cloning approaches for co-expression and production of recombinant protein-protein complexes

Multiprotein complexes play essential roles in all cells and X-ray crystallography can provide unparalleled insight into their structure and function. Many of these complexes are believed to be sufficiently stable for structural biology studies, but the production of protein-protein complexes using recombinant technologies is still labor-intensive. We have explored several strategies for the identification and cloning of heterodimers and heterotrimers that are compatible with the high-throughput (HTP) structural biology pipeline developed for single proteins. Two approaches are presented and compared which resulted in co-expression of paired genes from a single expression vector. Native operons encoding predicted interacting proteins were selected from a repertoire of genomes, and cloned directly to expression vector. In an alternative approach, Helicobacter pylori proteins predicted to interact strongly were cloned, each associated with translational control elements, then linked into an artificial operon. Proteins were then expressed and purified by standard HTP protocols, resulting to date in the structure determination of two H. pylori complexes.

Modulation of genetic clusters for synthesis of bioactive molecules in fungal endophytes: A review

Novel drugs with unique and targeted mode of action are very much need of the hour to treat and manage severe multidrug infections and other life-threatening complications. Though natural molecules have proved to be effective and environmentally safe, the relative paucity of discovery of new drugs has forced us to lean towards synthetic chemistry for developing novel drug molecules. Plants and microbes are the major resources that we rely upon in our pursuit towards discovery of novel compounds of pharmacological importance with less toxicity. Endophytes, an eclectic group of microbes having the potential to chemically bridge the gap between plants and microbes, have attracted the most attention due to their relatively high metabolic versatility. Since continuous large scale supply of major metabolites from microfungi and especially endophytes is severely impeded by the phenomenon of attenuation in axenic cultures, the major challenge is to understand the regulatory mechanisms in operation that drive the expression of metabolic gene clusters of pharmaceutical importance. This review is focused on the major regulatory elements that operate in filamentous fungi and various combinatorial multi-disciplinary approaches involving bioinformatics, molecular biology, and metabolomics that could aid in large scale synthesis of important lead molecules.

Footprints of optimal protein assembly strategies in the operonic structure of prokaryotes

In this work, we investigate optimality principles behind synthesis strategies for protein complexes using a dynamic optimization approach. We show that the cellular capacity of protein synthesis has a strong influence on optimal synthesis strategies reaching from a simultaneous to a sequential synthesis of the subunits of a protein complex. Sequential synthesis is preferred if protein synthesis is strongly limited, whereas a simultaneous synthesis is optimal in situations with a high protein synthesis capacity. We confirm the predictions of our optimization approach through the analysis of the operonic organization of protein complexes in several hundred prokaryotes. Thereby, we are able to show that cellular protein synthesis capacity is a driving force in the dissolution of operons comprising the subunits of a protein complex. Thus, we also provide a tested hypothesis explaining why the subunits of many prokaryotic protein complexes are distributed across several operons despite the presumably less precise co-regulation.

Co-regulation of translation in protein complexes

Background

Co-regulation of gene expression has been known for many years, and studied widely both globally and for individual genes. Nevertheless, most analyses concerned transcriptional control, which in case of physically interacting proteins and protein complex subunits may be of secondary importance. This research is the first quantitative analysis that provides global-scale evidence for translation co-regulation among associated proteins.

Results

By analyzing the results of our previous quantitative model of translation, we have demonstrated that protein production rates plus several other translational parameters, such as mRNA and protein abundance, or number of produced proteins from a gene, are well concerted between stable complex subunits and party hubs. This may be energetically favorable during synthesis of complex building blocks and ensure their accurate production in time. In contrast, for connections with regulatory particles and date hubs translational co-regulation is less visible, indicating that in these cases maintenance of accurate levels of interacting particles is not necessarily beneficial.

Conclusions

Similar results obtained for distantly related model organisms, Saccharomyces cerevisiae and Homo sapiens, suggest that the phenomenon of translational co-regulation applies to the variety of living organisms and concerns many complex constituents. This phenomenon was also observed among the set of functionally linked proteins from Escherichia coli operons. This leads to the conclusion that translational regulation of a protein should always be studied with respect to the expression of its primary interacting partners.

Reviewers

This article was reviewed by Sandor Pongor and Claus Wilke.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0048-7) contains supplementary material, which is available to authorized users.

Regulation of metabolic gene clusters in Arabidopsis thaliana

Recent discoveries have revealed that the genes for the biosynthesis of a variety of plant specialized metabolites are organized in operon-like clusters within plant genomes. Here we identify a regulatory process that is required for normal expression of metabolic gene clusters in Arabidopsis thaliana. Comparative gene expression analysis of a representative clustered gene was performed in a set of chromatin mutant lines. Subsequently, metabolite levels were analysed by GC-MS and the local chromatin structure was investigated by chromatin immunoprecipitation and nucleosome positioning. We show that the transcript levels of genes within two metabolic clusters are coordinately reduced in an arp6 and h2a.z background. We demonstrate that H2A.Z enrichment in the clusters is positively correlated with active cluster expression. We further show that nucleosome stability within the cluster regions is higher in the arp6 background compared with the wild-type. These results implicate ARP6 and H2A.Z in the regulation of metabolic clusters in Arabidopsis thaliana through localized chromatin modifications that enable the coordinate expression of groups of contiguous genes. These findings shed light on the complex process of cluster regulation, an area that could in the future open up new opportunities for the discovery and manipulation of specialized metabolic pathways in plants.

Metabolic engineering to enhance the value of plants as green factories

The promise of plants to serve as the green factories of the future is ever increasing. Plants have been used traditionally for construction, energy, food and feed. Bioactive compounds primarily derived from specialized plant metabolism continue to serve as important scaffold molecules for pharmaceutical drug production. Yet, the past few years have witnessed a growing interest on plants as the ultimate harvesters of carbon and energy from the sun, providing carbohydrate and lipid biofuels that would contribute to balancing atmospheric carbon. How can the metabolic output from plants be increased even further, and what are the bottlenecks? Here, we present what we perceive to be the main opportunities and challenges associated with increasing the efficiency of plants as chemical factories. We offer some perspectives on when it makes sense to use plants as production systems because the amount of biomass needed makes any other system unfeasible. However, there are other instances in which plants serve as great sources of biological catalysts, yet are not necessarily the best-suited systems for production. We also present emerging opportunities for manipulating plant genomes to make plant synthetic biology a reality.

Computational identification of a new SelD-like family that may participate in sulfur metabolism in hyperthermophilic sulfur-reducing archaea

BACKGROUND

Selenium (Se) and sulfur (S) are closely related elements that exhibit similar chemical properties. Some genes related to S metabolism are also involved in Se utilization in many organisms. However, the evolutionary relationship between the two utilization traits is unclear.

RESULTS

In this study, we conducted a comparative analysis of the selenophosphate synthetase (SelD) family, a key protein for all known Se utilization traits, in all sequenced archaea. Our search showed a very limited distribution of SelD and Se utilization in this kingdom. Interestingly, a SelD-like protein was detected in two orders of Crenarchaeota: Sulfolobales and Thermoproteales. Sequence and phylogenetic analyses revealed that SelD-like protein contains the same domain and conserved functional residues as those of SelD, and might be involved in S metabolism in these S-reducing organisms. Further genome-wide analysis of patterns of gene occurrence in different thermoproteales suggested that several genes, including SirA-like, Prx-like and adenylylsulfate reductase, were strongly related to SelD-like gene. Based on these findings, we proposed a simple model wherein SelD-like may play an important role in the biosynthesis of certain thiophosphate compound.

CONCLUSIONS

Our data suggest novel genes involved in S metabolism in hyperthermophilic S-reducing archaea, and may provide a new window for understanding the complex relationship between Se and S metabolism in archaea.

Overlapping genes: a window on gene evolvability

Background

The forces underlying genome architecture and organization are still only poorly understood in detail. Overlapping genes (genes partially or entirely overlapping) represent a genomic feature that is shared widely across biological organisms ranging from viruses to multi-cellular organisms. In bacteria, a third of the annotated genes are involved in an overlap. Despite the widespread nature of this arrangement, its evolutionary origins and biological ramifications have so far eluded explanation.

Results

Here we present a comparative approach using information from 699 bacterial genomes that sheds light on the evolutionary dynamics of overlapping genes. We show that these structures exhibit high levels of plasticity.

Conclusions

We propose a simple model allowing us to explain the observed properties of overlapping genes based on the importance of initiation and termination of transcriptional and translational processes. We believe that taking into account the processes leading to the expression of protein-coding genes hold the key to the understanding of overlapping genes structures.

Clustering of gene ontology terms in genomes

Although protein coding genes occupy only a small fraction of genomes in higher species, they are not randomly distributed within or between chromosomes. Clustering of genes with related function(s) and/or characteristics has been evident at several different levels. To study how common the clustering of functionally related genes is and what kind of functions the end products of these genes are involved, we collected gene ontology (GO) terms for complete genomes and developed a method to detect previously undefined gene clustering. Exhaustive analysis was performed for seven widely studied species ranging from human to Escherichia coli. To overcome problems related to varying gene lengths and densities, a novel method was developed and a fixed number of genes were analyzed irrespective of the genome span covered. Statistically very significant GO term clustering was apparent in all the investigated genomes. The analysis window, which ranged from 5 to 50 consecutive genes, revealed extensive GO term clusters for genes with widely varying functions. Here, the most interesting and significant results are discussed and the complete dataset for each analyzed species is available at the GOme database at http://bioinf.uta.fi/GOme. The results indicated that clusters of genes with related functions are very common, not only in bacteria, in which operons are frequent, but also in all the studied species irrespective of how complex they are. There are some differences between species but in all of them GO term clusters are common and of widely differing sizes. The presented method can be applied to analyze any genome or part of a genome for which descriptive features are available, and thus is not restricted to ontology terms. This method can also be applied to investigate gene and protein expression patterns. The results pave a way for further studies of mechanisms that shape genome structure and evolutionary forces related to them.

The rise of operon-like gene clusters in plants

Gene clusters are common features of prokaryotic genomes also present in eukaryotes. Most clustered genes known are involved in the biosynthesis of secondary metabolites. Although horizontal gene transfer is a primary source of prokaryotic gene cluster (operon) formation and has been reported to occur in eukaryotes, the predominant source of cluster formation in eukaryotes appears to arise de novo or through gene duplication followed by neo- and sub-functionalization or translocation. Here we aim to provide an overview of the current knowledge and open questions related to plant gene cluster functioning, assembly, and regulation. We also present potential research approaches and point out the benefits of a better understanding of gene clusters in plants for both fundamental and applied plant science.

The precarious prokaryotic chromosome

Evolutionary selection for optimal genome preservation, replication, and expression should yield similar chromosome organizations in any type of cells. And yet, the chromosome organization is surprisingly different between eukaryotes and prokaryotes. The nuclear versus cytoplasmic accommodation of genetic material accounts for the distinct eukaryotic and prokaryotic modes of genome evolution, but it falls short of explaining the differences in the chromosome organization. I propose that the two distinct ways to organize chromosomes are driven by the differences between the global-consecutive chromosome cycle of eukaryotes and the local-concurrent chromosome cycle of prokaryotes. Specifically, progressive chromosome segregation in prokaryotes demands a single duplicon per chromosome, while other "precarious" features of the prokaryotic chromosomes can be viewed as compensations for this severe restriction.

Comprehensive analysis of the Corynebacterium glutamicum transcriptome using an improved RNAseq technique

Background

The use of RNAseq to resolve the transcriptional organization of an organism was established in recent years and also showed the complexity and dynamics of bacterial transcriptomes. The aim of this study was to comprehensively investigate the transcriptome of the industrially relevant amino acid producer and model organism Corynebacterium glutamicum by RNAseq in order to improve its genome annotation and to describe important features for transcription and translation.

Results

RNAseq data sets were obtained by two methods, one that focuses on 5′-ends of primary transcripts and another that provides the overall transcriptome with an improved resolution of 3′-ends of transcripts. Subsequent data analysis led to the identification of more than 2,000 transcription start sites (TSSs), the definition of 5′-UTRs (untranslated regions) for annotated protein-coding genes, operon structures and many novel transcripts located between or in antisense orientation to protein-coding regions. Interestingly, a high number of mRNAs (33%) is transcribed as leaderless transcripts. From the data, consensus promoter and ribosome binding site (RBS) motifs were identified and it was shown that the majority of genes in C. glutamicum are transcribed monocistronically, but operons containing up to 16 genes are also present.

Conclusions

The comprehensive transcriptome map of C. glutamicum established in this study represents a major step forward towards a complete definition of genetic elements (e.g. promoter regions, gene starts and stops, 5′-UTRs, RBSs, transcript starts and ends) and provides the ideal basis for further analyses on transcriptional regulatory networks in this organism. The methods developed are easily applicable for other bacteria and have the potential to be used also for quantification of transcriptomes, replacing microarrays in the near future.

A metabolic gene cluster in Lotus japonicus discloses novel enzyme functions and products in triterpene biosynthesis

Genes for triterpene biosynthetic pathways exist as metabolic gene clusters in oat and Arabidopsis thaliana plants. We characterized the presence of an analogous gene cluster in the model legume Lotus japonicus. In the genomic regions flanking the oxidosqualene cyclase AMY2 gene, genes for two different classes of cytochrome P450 and a gene predicted to encode a reductase were identified. Functional characterization of the cluster genes was pursued by heterologous expression in Nicotiana benthamiana. The gene expression pattern was studied under different developmental and environmental conditions. The physiological role of the gene cluster in nodulation and plant development was studied in knockdown experiments. A novel triterpene structure, dihydrolupeol, was produced by AMY2. A new plant cytochrome P450, CYP71D353, which catalyses the formation of 20-hydroxybetulinic acid in a sequential three-step oxidation of 20-hydroxylupeol was characterized. The genes within the cluster are highly co-expressed during root and nodule development, in hormone-treated plants and under various environmental stresses. A transcriptional gene silencing mechanism that appears to be involved in the regulation of the cluster genes was also revealed. A tightly co-regulated cluster of functionally related genes is involved in legume triterpene biosynthesis, with a possible role in plant development.

CluGene: A Bioinformatics Framework for the Identification of Co-Localized, Co-Expressed and Co-Regulated Genes Aimed at the Investigation of Transcriptional Regulatory Networks from High-Throughput Expression Data

The full understanding of the mechanisms underlying transcriptional regulatory networks requires unravelling of complex causal relationships. Genome high-throughput technologies produce a huge amount of information pertaining gene expression and regulation; however, the complexity of the available data is often overwhelming and tools are needed to extract and organize the relevant information. This work starts from the assumption that the observation of co-occurrent events (in particular co-localization, co-expression and co-regulation) may provide a powerful starting point to begin unravelling transcriptional regulatory networks. Co-expressed genes often imply shared functional pathways; co-expressed and functionally related genes are often co-localized, too; moreover, co-expressed and co-localized genes are also potential targets for co-regulation; finally, co-regulation seems more frequent for genes mapped to proximal chromosome regions. Despite the recognized importance of analysing co-occurrent events, no bioinformatics solution allowing the simultaneous analysis of co-expression, co-localization and co-regulation is currently available. Our work resulted in developing and valuating CluGene, a software providing tools to analyze multiple types of co-occurrences within a single interactive environment allowing the interactive investigation of combined co-expression, co-localization and co-regulation of genes. The use of CluGene will enhance the power of testing hypothesis and experimental approaches aimed at unravelling transcriptional regulatory networks. The software is freely available at http://bioinfolab.unipg.it/.

Genomic positions of co-expressed genes: echoes of chromosome organisation in gene expression data

Background

The relationships between gene expression and nuclear structure, chromosome territories in particular, are currently being elucidated experimentally. Each chromosome occupies an individual, spatially-limited space with a preferential position relative to the nuclear centre that may be specific to the cell and tissue type. We sought to discover whether patterns in gene expression databases might exist that would mirror prevailing or recurring nuclear structure patterns, chromosome territory interactions in particular.

Results

We used human gene expression datasets, both from a tissue expression atlas and from a large set including diverse types of perturbations. We identified groups of positional gene clusters over-represented in gene expression clusters. We show that some pairs of chromosomes and pairs of 10 Mbp long chromosome regions are significantly enriched in the expression clusters. The functions of genes involved in inter-chromosome co-expression relationships are non-random and predominantly related to cell-cell communication and reaction to external stimuli.

Conclusions

We suggest that inter-chromosomal gene co-expression can be interpreted in the context of nuclear structure, and that even expression datasets that include very diverse conditions and cell types show consistent relationships.

Evolution of a complex locus for terpene biosynthesis in solanum

Functional gene clusters, containing two or more genes encoding different enzymes for the same pathway, are sometimes observed in plant genomes, most often when the genes specify the synthesis of specialized defensive metabolites. Here, we show that a cluster of genes in tomato (Solanum lycopersicum; Solanaceae) contains genes for terpene synthases (TPSs) that specify the synthesis of monoterpenes and diterpenes from cis-prenyl diphosphates, substrates that are synthesized by enzymes encoded by cis-prenyl transferase (CPT) genes also located within the same cluster. The monoterpene synthase genes in the cluster likely evolved from a diterpene synthase gene in the cluster by duplication and divergence. In the orthologous cluster in Solanum habrochaites, a new sesquiterpene synthase gene was created by a duplication event of a monoterpene synthase followed by a localized gene conversion event directed by a diterpene synthase gene. The TPS genes in the Solanum cluster encoding cis-prenyl diphosphate-utilizing enzymes are closely related to a tobacco (Nicotiana tabacum; Solanaceae) diterpene synthase encoding Z-abienol synthase (Nt-ABS). Nt-ABS uses the substrate copal-8-ol diphosphate, which is made from the all-trans geranylgeranyl diphosphate by copal-8-ol diphosphate synthase (Nt-CPS2). The Solanum gene cluster also contains an ortholog of Nt-CPS2, but it appears to encode a nonfunctional protein. Thus, the Solanum functional gene cluster evolved by duplication and divergence of TPS genes, together with alterations in substrate specificity to utilize cis-prenyl diphosphates and through the acquisition of CPT genes.

Phylogenetic analysis of erythritol catabolic loci within the Rhizobiales and proteobacteria

BACKGROUND

The ability to use erythritol as a sole carbon source is not universal among the Rhizobiaceae. Based on the relatedness to the catabolic genes in Brucella it has been suggested that the eryABCD operon may have been horizontally transferred into Rhizobium. During work characterizing a locus necessary for the transport and catabolism of erythritol, adonitol and L-arabitol in Sinorhizobium meliloti, we became interested in the differences between the erythritol loci of S. meliloti and R. leguminosarum. Utilizing the Ortholog Neighborhood Viewer from the DOE Joint Genome Institute database it appeared that loci for erythritol and polyol utilization had distinct arrangements that suggested these loci may have undergone genetic rearrangements.

RESULTS

A data set was established of genetic loci containing erythritol/polyol orthologs for 19 different proteobacterial species. These loci were analyzed for genetic content and arrangement of genes associated with erythritol, adonitol and L-arabitol catabolism. Phylogenetic trees were constructed for core erythritol catabolic genes and contrasted with the species phylogeny. Additionally, phylogenetic trees were constructed for genes that showed differences in arrangement among the putative erythritol loci in these species.

CONCLUSIONS

Three distinct erythritol/polyol loci arrangements have been identified that reflect metabolic need or specialization. Comparison of the phylogenetic trees of core erythritol catabolic genes with species phylogeny provides evidence that is consistent with these loci having been horizontally transferred from the alpha-proteobacteria into both the beta and gamma-proteobacteria. ABC transporters within these loci adopt 2 unique genetic arrangements, and although biological data suggests they are functional erythritol transporters, phylogenetic analysis suggests they may not be orthologs and probably should be considered analogs. Finally, evidence for the presence of paralogs, and xenologs of erythritol catabolic genes in some of the genomes included in the analysis is provided.

Regulation of internal promoters in a zinc-responsive operon is influenced by transcription from upstream promoters

In the cyanobacterium Anabaena sp. strain PCC 7120 (also known as Nostoc sp. strain PCC 7120), a zinc-responsive operon (all4725-all4721) has been described, which contains 4 distinct promoters. The two most upstream ones bind Zur with high affinity, whereas the other two do not or do so with a very low affinity. In this paper, a detailed characterization of the four promoters is presented, showing that all four were induced by metal depletion, and they were constitutively derepressed in a zur mutant, despite the two downstream promoters not being direct targets for this regulator. Crucially, induction by metal depletion of the two downstream promoters was abrogated when transcription initiated at the upstream promoters was interrupted by a polar insertion midway in the operon. In contrast, insertion of a nitrogen-responsive promoter at a roughly similar position provoked the two downstream promoters to adopt a regulatory pattern mimicking that of the inserted promoter. Thus, regulation of the two downstream promoters is apparently influenced by transcription from promoters upstream. Evidence is presented indicating that the activity of the two downstream promoters is kept basal in Anabaena by repression. A regulatory model compatible with these results is proposed, where promoters controlled by repression in bacterial operons may be subjected to a hierarchical regulation depending on their position in the operon. According to this model, internal promoters may respond to stimuli governing the activity of promoters upstream by an indirect regulation and to specific stimuli by a direct regulation.