Structure and dynamics of the operon map of Buchnera aphidicola sp. strain APS

OUP accepted manuscript

Operons are a hallmark of the genomic and regulatory architecture of prokaryotes. However, the mechanism by which two genes placed far apart gradually come close and form operons remains to be elucidated. Here, we propose a new model of the origin of operons: Mobile genetic elements called insertion sequences can facilitate the formation of operons by consecutive insertion–deletion–excision reactions. This mechanism barely leaves traces of insertion sequences and thus difficult to detect in nature. In this study, as a proof-of-concept, we reproducibly demonstrated operon formation in the laboratory. The insertion sequence IS3 and the insertion sequence excision enhancer are genes found in a broad range of bacterial species. We introduced these genes into insertion sequence-less Escherichia coli and found that, supporting our hypothesis, the activity of the two genes altered the expression of genes surrounding IS3, closed a 2.7 kb gap between a pair of genes, and formed new operons. This study shows how insertion sequences can facilitate the rapid formation of operons through locally increasing the structural mutation rates and highlights how coevolution with mobile elements may shape the organization of prokaryotic genomes and gene regulation.

Stable isotope studies reveal pathways for the incorporation of non-essential amino acids in Acyrthosiphon pisum (pea aphids)

Plant roots incorporate inorganic nitrogen into the amino acids glutamine, glutamic acid, asparagine and aspartic acid, which together serve as the primary metabolites of nitrogen transport to other tissues. Given the preponderance of these four amino acids, phloem sap is a nutritionally unbalanced diet for phloem-feeding insects. Therefore, aphids and other phloem feeders typically rely on microbial symbionts for the synthesis of essential amino acids. To investigate the metabolism of the four main transport amino acids by the pea aphid (Acyrthosiphon pisum), and its Buchnera aphidicola endosymbionts, aphids were fed defined diets with stable isotope-labeled glutamine, glutamic acid, asparagine or aspartic acid (U-(13)C, U-(15)N; U-(15)N; α-(15)N; or γ-(15)N). The metabolic fate of the dietary (15)N and (13)C was traced using gas chromatography-mass spectrometry (GC-MS). Nitrogen was the major contributor to the observed amino acid isotopomers with one additional unit mass (M+1). However, there was differential incorporation, with the amine nitrogen of asparagine being incorporated into other amino acids more efficiently than the amide nitrogen. Higher isotopomers (M+2, M+3 and M+4) indicated the incorporation of varying numbers of (13)C atoms into essential amino acids. GC-MS assays also showed that, even with an excess of dietary labeled glutamine, glutamic acid, asparagine or aspartic acid, the overall content of these amino acids in aphid bodies was mostly the product of catabolism of dietary amino acids and subsequent re-synthesis within the aphids. Thus, these predominant dietary amino acids are not passed directly to Buchnera endosymbionts for synthesis of essential amino acids, but are rather are produced de novo, most likely by endogenous aphid enzymes.

Novel Insights into Insect-Microbe Interactions-Role of Epigenomics and Small RNAs

It has become increasingly clear that microbes form close associations with the vast majority of animal species, especially insects. In fact, an array of diverse microbes is known to form shared metabolic pathways with their insect hosts. A growing area of research in insect-microbe interactions, notably for hemipteran insects and their mutualistic symbionts, is to elucidate the regulation of this inter-domain metabolism. This review examines two new emerging mechanisms of gene regulation and their importance in host-microbe interactions. Specifically, we highlight how the incipient areas of research on regulatory "dark matter" such as epigenomics and small RNAs, can play a pivotal role in the evolution of both insect and microbe gene regulation. We then propose specific models of how these dynamic forms of gene regulation can influence insect-symbiont-plant interactions. Future studies in this area of research will give us a systematic understanding of how these symbiotic microbes and animals reciprocally respond to and regulate their shared metabolic processes.

Cautions about the reliability of pairwise gene correlations based on expression data

BACKGROUND

Rapid growth in the availability of genome-wide transcript abundance levels through gene expression microarrays and RNAseq promises to provide deep biological insights into the complex, genome-wide transcriptional behavior of single-celled organisms. However, this promise has not yet been fully realized.

RESULTS

We find that computation of pairwise gene associations (correlation; mutual information) across a set of 2782 total genome-wide expression samples from six diverse bacteria produces unexpectedly large variation in estimates of pairwise gene association-regardless of the metric used, the organism under study, or the number and source of the samples. We pinpoint the cause to sampling bias. In particular, in repositories of expression data (e.g., Gene Expression Omnibus, GEO), many individual genes show small differences in absolute gene expression levels across the set of samples. We demonstrate that these small differences are due mainly to "noise" instead of "signal" attributable to environmental or genetic perturbations. We show that downstream analysis using gene expression levels of genes with small differences yields biased estimates of pairwise association.

CONCLUSIONS

We propose flagging genes with small differences in absolute, RMA-normalized, expression levels (e.g., standard deviation less than 0.5), as potentially yielding biased pairwise association metrics. This strategy has the potential to substantially improve the confidence in genome-wide conclusions about transcriptional behavior in bacterial organisms. Further work is needed to further refine strategies to identify genes with small difference in expression levels prior to computing gene-gene association metrics.

Natural selection for operons depends on genome size

In prokaryotes, genome size is associated with metabolic versatility, regulatory complexity, effective population size, and horizontal transfer rates. We therefore analyzed the covariation of genome size and operon conservation to assess the evolutionary models of operon formation and maintenance. In agreement with previous results, intraoperonic pairs of essential and of highly expressed genes are more conserved. Interestingly, intraoperonic pairs of genes are also more conserved when they encode proteins at similar cell concentrations, suggesting a role of cotranscription in diminishing the cost of waste and shortfall in gene expression. Larger genomes have fewer and smaller operons that are also less conserved. Importantly, lower conservation in larger genomes was observed for all classes of operons in terms of gene expression, essentiality, and balanced protein concentration. We reached very similar conclusions in independent analyses of three major bacterial clades (α- and β-Proteobacteria and Firmicutes). Operon conservation is inversely correlated to the abundance of transcription factors in the genome when controlled for genome size. This suggests a negative association between the complexity of genetic networks and operon conservation. These results show that genome size and/or its proxies are key determinants of the intensity of natural selection for operon organization. Our data fit better the evolutionary models based on the advantage of coregulation than those based on genetic linkage or stochastic gene expression. We suggest that larger genomes with highly complex genetic networks and many transcription factors endure weaker selection for operons than smaller genomes with fewer alternative tools for genetic regulation.

Genomic analysis of the regulatory elements and links with intrinsic DNA structural properties in the shrunken genome of Buchnera

Background

Buchnera aphidicola is an obligate symbiotic bacterium, associated with most of the aphididae, whose genome has drastically shrunk during intracellular evolution. Gene regulation in Buchnera has been a matter of controversy in recent years as the combination of genomic information with the experimental results has been contradictory, refuting or arguing in favour of a functional and responsive transcription regulation in Buchnera.

The goal of this study was to describe the gene transcription regulation capabilities of Buchnera based on the inventory of cis- and trans-regulators encoded in the genomes of five strains from different aphids (Acyrthosiphon pisum, Schizaphis graminum, Baizongia pistacea, Cinara cedri and Cinara tujafilina), as well as on the characterisation of some intrinsic structural properties of the DNA molecule in these bacteria.

Results

Interaction graph analysis shows that gene neighbourhoods are conserved between E. coli and Buchnera in structures called transcriptons, interactons and metabolons, indicating that selective pressures have acted on the evolution of transcriptional, protein-protein interaction and metabolic networks in Buchnera. The transcriptional regulatory network in Buchnera is composed of a few general DNA-topological regulators (Nucleoid Associated Proteins and topoisomerases), with the quasi-absence of any specific ones (except for multifunctional enzymes with a known gene expression regulatory role in Escherichia coli, such as AlaS, PepA and BolA, and the uncharacterized hypothetical regulators YchA and YrbA). The relative positioning of regulatory genes along the chromosome of Buchnera seems to have conserved its ancestral state, despite the genome erosion. Sigma-70 promoters with canonical thermodynamic sequence profiles were detected upstream of about 94% of the CDS of Buchnera in the different aphids. Based on Stress-Induced Duplex Destabilization (SIDD) measurements, unstable σ⁷⁰ promoters were found specifically associated with the regulator and transporter genes.

Conclusions

This genomic analysis provides supporting evidence of a selection of functional regulatory structures and it has enabled us to propose hypotheses concerning possible links between these regulatory elements and the DNA-topology (i.e., supercoiling, curvature, flexibility and base-pair stability) in the regulation of gene expression in the shrunken genome of Buchnera.

Multimodal dynamic response of the Buchnera aphidicola pLeu plasmid to variations in leucine demand of its host, the pea aphid Acyrthosiphon pisum

Aphids, important agricultural pests, can grow and reproduce thanks to their intimate symbiosis with the γ-proteobacterium Buchnera aphidicola that furnishes them with essential amino acids lacking in their phloem sap diet. To study how B. aphidicola, with its reduced genome containing very few transcriptional regulators, responds to variations in the metabolic requirements of its host, we concentrated on the leucine metabolic pathway. We show that leucine is a limiting factor for aphid growth and it displays a stimulatory feeding effect. Our metabolic analyses demonstrate that symbiotic aphids are able to respond to leucine starvation or excess by modulating the neosynthesis of this amino acid. At a molecular level, this response involves an early important transcriptional regulation (after 12 h of treatment) followed by a moderate change in the pLeu plasmid copy number. Both responses are no longer apparent after 7 days of treatment. These experimental data are discussed in the light of a re-annotation of the pLeu plasmid regulatory elements. Taken together, our data show that the response of B. aphidicola to the leucine demand of its host is multimodal and dynamically regulated, providing new insights concerning the genetic regulation capabilities of this bacterium in relation to its symbiotic functions.