Complex regulation in a Comamonas platform for diverse aromatic carbon metabolism.

Critical to a sustainable energy future are microbial platforms that can process aromatic carbons from the largely untapped reservoir of lignin and plastic feedstocks. Comamonas species present promising bacterial candidates for such platforms because they can use a range of natural and xenobiotic aromatic compounds and often possess innate genetic constraints that avoid competition with sugars. However, the metabolic reactions of these species are underexplored, and the regulatory mechanisms are unknown. Here we identify multilevel regulation in the conversion of lignin-related natural aromatic compounds, 4-hydroxybenzoate and vanillate, and the plastics-related xenobiotic aromatic compound, terephthalate, in Comamonas testosteroni KF-1. Transcription-level regulation controls initial catabolism and cleavage, but metabolite-level thermodynamic regulation governs fluxes in central carbon metabolism. Quantitative 13C mapping of tricarboxylic acid cycle and cataplerotic reactions elucidates key carbon routing not evident from enzyme abundance changes. This scheme of transcriptional activation coupled with metabolic fine-tuning challenges outcome predictions during metabolic manipulations.

Haplotype divergence supports long-term asexuality in the oribatid mite Oppiella nova.

Sex strongly impacts genome evolution via recombination and segregation. In the absence of these processes, haplotypes within lineages of diploid organisms are predicted to accumulate mutations independently of each other and diverge over time. This so-called "Meselson effect" is regarded as a strong indicator of the long-term evolution under obligate asexuality. Here, we present genomic and transcriptomic data of three populations of the asexual oribatid mite species Oppiella nova and its sexual relative Oppiella subpectinata We document strikingly different patterns of haplotype divergence between the two species, strongly supporting Meselson effect-like evolution and long-term asexuality in O. nova: I) variation within individuals exceeds variation between populations in O. nova but vice versa in O. subpectinata; II) two O. nova sublineages feature a high proportion of lineage-specific heterozygous single-nucleotide polymorphisms (SNPs), indicating that haplotypes continued to diverge after lineage separation; III) the deepest split in gene trees generally separates the two haplotypes in O. nova, but populations in O. subpectinata; and IV) the topologies of the two haplotype trees match each other. Our findings provide positive evidence for the absence of canonical sex over evolutionary time in O. nova and suggest that asexual oribatid mites can escape the dead-end fate usually associated with asexual lineages.

Evidence for cryptic sex in parthenogenetic stick insects of the genus Timema.

Obligately parthenogenetic species are expected to be short lived since the lack of sex and recombination should translate into a slower adaptation rate and increased accumulation of deleterious alleles. Some, however, are thought to have been reproducing without males for millions of years. It is not clear how these old parthenogens can escape the predicted long-term costs of parthenogenesis, but an obvious explanation is cryptic sex. In this study, we screen for signatures of cryptic sex in eight populations of four parthenogenetic species of Timema stick insects, some estimated to be older than 1 Myr. Low genotype diversity, homozygosity of individuals and high linkage disequilibrium (LD) unaffected by marker distances support exclusively parthenogenetic reproduction in six populations. However, in two populations (namely, of the species Timema douglasi and T. monikensis) we find strong evidence for cryptic sex, most likely mediated by rare males. These populations had comparatively high genotype diversities, lower LD, and a clear LD decay with genetic distance. Rare sex in species that are otherwise largely parthenogenetic could help explain the unusual success of parthenogenesis in the Timema genus and raises the question whether episodes of rare sex are in fact the simplest explanation for the persistence of many old parthenogens in nature.

Spurious regulatory connections dictate the expression-fitness landscape of translation factors.

During steady-state cell growth, individual enzymatic fluxes can be directly inferred from growth rate by mass conservation, but the inverse problem remains unsolved. Perturbing the flux and expression of a single enzyme could have pleiotropic effects that may or may not dominate the impact on cell fitness. Here, we quantitatively dissect the molecular and global responses to varied expression of translation termination factors (peptide release factors, RFs) in the bacterium Bacillus subtilis. While endogenous RF expression maximizes proliferation, deviations in expression lead to unexpected distal regulatory responses that dictate fitness reduction. Molecularly, RF depletion causes expression imbalance at specific operons, which activates master regulators and detrimentally overrides the transcriptome. Through these spurious connections, RF abundances are thus entrenched by focal points within the regulatory network, in one case located at a single stop codon. Such regulatory entrenchment suggests that predictive bottom-up models of expression-fitness landscapes will require near-exhaustive characterization of parts.

X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species.

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

A genome-wide investigation of adaptive signatures in protein-coding genes related to tool behaviour in New Caledonian and Hawaiian crows.

Very few animals habitually manufacture and use tools. It has been suggested that advanced tool behaviour co-evolves with a suite of behavioural, morphological and life history traits. In fact, there are indications for such an adaptive complex in tool-using crows (genus Corvus species). Here, we sequenced the genomes of two habitually tool-using and ten non-tool-using crow species to search for genomic signatures associated with a tool-using lifestyle. Using comparative genomic and population genetic approaches, we screened for signals of selection in protein-coding genes in the tool-using New Caledonian and Hawaiian crows. While we detected signals of recent selection in New Caledonian crows near genes associated with bill morphology, our data indicate that genetic changes in these two lineages are surprisingly subtle, with little evidence at present for convergence. We explore the biological explanations for these findings, such as the relative roles of gene regulation and protein-coding changes, as well as the possibility that statistical power to detect selection in recently diverged lineages may have been insufficient. Our study contributes to a growing body of literature aiming to decipher the genetic basis of recently evolved complex behaviour.

Sex-specific responses to cold in a very cold-tolerant, northern Drosophila species.

Organisms can plastically alter resource allocation in response to changing environmental factors. For example, in harsh conditions, organisms are expected to shift investment from reproduction toward survival; however, the factors and mechanisms that govern the magnitude of such shifts are relatively poorly studied. Here we compared the impact of cold on males and females of the highly cold-tolerant species Drosophila montana at the phenotypic and transcriptomic levels. Although both sexes showed similar changes in cold tolerance and gene expression in response to cold treatment, indicating that the majority of changes are concordant between the sexes, we identified a clear reduction in sexually dimorphic gene expression, suggesting that preparing for the colder season involves reducing investment in sex-specific traits. This reduction was larger in males than females, as expected if male sexual traits are more condition-dependent than female traits, as predicted by theory. Gene expression changes were primarily associated with shifts in metabolic profile, which likely play a role in increasing cold tolerance. Finally, we found that the expression of immune genes was reduced following cold treatment, suggesting that reduced investment in costly immune function may be important in helping flies survive colder periods.

Repeated Evolution of Asexuality Involves Convergent Gene Expression Changes.

Asexual reproduction has evolved repeatedly from sexual ancestors across a wide range of taxa. Whereas the costs and benefits associated with asexuality have received considerable attention, the molecular changes underpinning the evolution of asexual reproduction remain relatively unexplored. In particular, it is completely unknown whether the repeated evolution of asexual phenotypes involves similar molecular changes, as previous studies have focused on changes occurring in single lineages. Here, we investigate the extent of convergent gene expression changes across five independent transitions to asexuality in stick insects. We compared gene expression of asexual females to females of close sexual relatives in whole-bodies, reproductive tracts, and legs. We identified a striking amount of convergent gene expression change (up to 8% of genes), greatly exceeding that expected by chance. Convergent changes were also tissue-specific, and most likely driven by selection for functional changes. Genes showing convergent changes in the reproductive tract were associated with meiotic spindle formation and centrosome organization. These genes are particularly interesting as they can influence the production of unreduced eggs, a key barrier to asexual reproduction. Changes in legs and whole-bodies were likely involved in female sexual trait decay, with enrichment in terms such as sperm-storage and pigmentation. By identifying changes occurring across multiple independent transitions to asexuality, our results provide a rare insight into the molecular basis of asexual phenotypes and suggest that the evolutionary path to asexuality is highly constrained, requiring repeated changes to the same key genes.

Fitness advantages conferred by the L20-interacting RNA cis-regulator of ribosomal protein synthesis in Bacillus subtilis.

In many bacteria, ribosomal proteins autogenously repress their own expression by interacting with RNA structures typically located in the 5'-UTRs of their mRNA transcripts. This regulation is necessary to maintain a balance between ribosomal proteins and rRNA to ensure proper ribosome production. Despite advances in noncoding RNA discovery and validation of RNA-protein regulatory interactions, the selective pressures that govern the formation and maintenance of such RNA cis-regulators in the context of an organism remain largely undetermined. To examine the impact disruptions to this regulation have on bacterial fitness, we introduced point mutations that abolish ribosomal protein binding and regulation into the RNA structure that controls expression of ribosomal proteins L20 and L35 within the Bacillus subtilis genome. Our studies indicate that removing this regulation results in reduced log phase growth, improper rRNA maturation, and the accumulation of a kinetically trapped or misassembled ribosomal particle at low temperatures, suggesting defects in ribosome synthesis. Such work emphasizes the important role regulatory RNAs play in the stoichiometric production of ribosomal components for proper ribosome composition and overall organism viability and reinforces the potential of targeting ribosomal protein production and ribosome assembly with novel antimicrobials.

Rapid Accumulation of Motility-Activating Mutations in Resting Liquid Culture of Escherichia coli.

Expression of motility genes is a potentially beneficial but costly process in bacteria. Interestingly, many isolate strains of Escherichia coli possess motility genes but have lost the ability to activate them under conditions in which motility is advantageous, raising the question of how they respond to these situations. Through transcriptome profiling of strains in the E. coli single-gene knockout Keio collection, we noticed drastic upregulation of motility genes in many of the deletion strains compared to levels in their weakly motile parent strain (BW25113). We show that this switch to a motile phenotype is not a direct consequence of the genes deleted but is instead due to a variety of secondary mutations that increase the expression of the major motility regulator, FlhDC. Importantly, we find that this switch can be reproduced by growing poorly motile E. coli strains in nonshaking liquid medium overnight but not in shaking liquid medium. Individual isolates after the nonshaking overnight incubations acquired distinct mutations upstream of the flhDC operon, including different insertion sequence (IS) elements and, to a lesser extent, point mutations. The rapidity with which genetic changes sweep through the populations grown without shaking shows that poorly motile strains can quickly adapt to a motile lifestyle by genetic rewiring.IMPORTANCE The ability to tune gene expression in times of need outside preordained regulatory networks is an essential evolutionary process that allows organisms to survive and compete. Here, we show that upon overnight incubation in liquid medium without shaking, populations of largely nonmotile Escherichia coli bacteria can rapidly accumulate mutants that have constitutive motility. This effect contributes to widespread secondary mutations in the single-gene knockout library, the Keio collection. As a result, 49/71 (69%) of the Keio strains tested exhibited various degrees of motility, whereas their parental strain is poorly motile. These observations highlight the plasticity of gene expression even in the absence of preexisting regulatory programs and should raise awareness of procedures for handling laboratory strains of E. coli.

Consequences of Asexuality in Natural Populations: Insights from Stick Insects.

Recombination is a fundamental process with significant impacts on genome evolution. Predicted consequences of the loss of recombination include a reduced effectiveness of selection, changes in the amount of neutral polymorphisms segregating in populations, and an arrest of GC-biased gene conversion. Although these consequences are empirically well documented for nonrecombining genome portions, it remains largely unknown if they extend to the whole genome scale in asexual organisms. We identify the consequences of asexuality using de novo transcriptomes of five independently derived, obligately asexual lineages of stick insects, and their sexual sister-species. We find strong evidence for higher rates of deleterious mutation accumulation, lower levels of segregating polymorphisms and arrested GC-biased gene conversion in asexuals as compared with sexuals. Taken together, our study conclusively shows that predicted consequences of genome evolution under asexuality can indeed be found in natural populations.

Validating the Demethylating Effects of 5-aza-2'-deoxycytidine in Insects Requires a Whole-Genome Approach.

We previously demonstrated that treatment with the demethylating agent 5-aza-2'-deoxycytidine (5-aza-dC) alters the offspring sex ratios produced by females of the parasitoid wasp Nasonia vitripennis. Females allocate offspring sex ratio in line with local mate competition theory, producing more or less female-biased sex ratios as the number of other females laying eggs on a patch varies, thereby reducing competition among their sons for mates. Interestingly, treatment with 5-aza-dC did not ablate the facultative sex allocation response. Instead, sex ratios became less female biased, a shift in the direction of the optimum sex ratio for paternally inherited alleles according to genomic conflict theory. This was the first (albeit indirect) experimental evidence for genomic conflict over sex allocation. In their comment, Ellers and colleagues assayed the effects of 5-aza-dC on DNA methylation in 10 Nasonia genes, finding no evidence of demethylation in these 10 genes, from which they conclude that 5-aza-dC has no demethylating capability in N. vitripennis. Quantifying the efficacy of 5-aza-dC in terms of demethylation is indeed crucial to in-depth interpretation of studies using 5-aza-dC to link phenotypes to epigenetic regulation. Here we outline the mode of action of 5-aza-dC and demonstrate that determining the efficacy of 5-aza-dC in insect systems requires a whole-genome approach.

Fundamental and realized feeding niche breadths of sexual and asexual stick insects.

The factors contributing to the maintenance of sex over asexuality in natural populations remain unclear. Ecological divergences between sexual and asexual lineages could help to maintain reproductive polymorphisms, at least transiently, but the consequences of asexuality for the evolution of ecological niches are unknown. Here, we investigated how niche breadths change in transitions from sexual reproduction to asexuality. We used host plant ranges as a proxy to compare the realized feeding niche breadths of five independently derived asexual Timema stick insect species and their sexual relatives at both the species and population levels. Asexual species had systematically narrower realized niches than sexual species, though this pattern was not apparent at the population level. To investigate how the narrower realized niches of asexual species arise, we performed feeding experiments to estimate fundamental niche breadths but found no systematic differences between reproductive modes. The narrow realized niches found in asexual species are therefore probably a consequence of biotic interactions such as predation or competition, that constrain realized niche size in asexuals more strongly than in sexuals.

Codon usage bias and effective population sizes on the X chromosome versus the autosomes in Drosophila melanogaster.

Codon usage bias (CUB) in Drosophila is higher for X-linked genes than for autosomal genes. One possible explanation is that the higher effective recombination rate for genes on the X chromosome compared with the autosomes reduces their susceptibility to Hill-Robertson effects, and thus enhances the efficacy of selection on codon usage. The genome sequence of D. melanogaster was used to test this hypothesis. Contrary to expectation, it was found that, after correcting for the effective recombination rate, CUB remained higher on the X than on the autosomes. In contrast, an analysis of polymorphism data from a Rwandan population showed that mean nucleotide site diversity at 4-fold degenerate sites for genes on the X is approximately three-quarters of the autosomal value after correcting for the effective recombination rate, compared with approximate equality before correction. In addition, these data show that selection for preferred versus unpreferred synonymous variants is stronger on the X than the autosomes, which accounts for the higher CUB of genes on the X chromosome. This difference in the strength of selection does not appear to reflect the effects of dominance of mutations affecting codon usage, differences in gene expression levels between X and autosomes, or differences in mutational bias. Its cause therefore remains unexplained. The stronger selection on CUB on the X chromosome leads to a lower rate of synonymous site divergence compared with the autosomes; this will cause a stronger upward bias for X than A in estimates of the proportion of nonsynonymous mutations fixed by positive selection, for methods based on the McDonald-Kreitman test.

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution.

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

DNA catalysis of a normally disfavored RNA hydrolysis reaction.

We recently used in vitro selection to identify many deoxyribozymes that catalyze DNA phosphodiester bond hydrolysis and create 5'-phosphate and 3'-hydroxyl termini. Alternatively, numerous deoxyribozymes have been identified for catalysis of RNA cleavage by 2'-hydroxyl transesterification, forming 2',3'-cyclic phosphate and 5'-hydroxyl termini. In this study, we investigated the ability of DNA to catalyze RNA cleavage by hydrolysis rather than transesterification, although normally the hydrolysis reaction is substantially disfavored relative to transesterification. Via a series of in vitro selection experiments, we found that reselection of a DNA-hydrolyzing deoxyribozyme leads either to transesterification or hydrolysis, depending on exclusion or inclusion of a stringent selection pressure for hydrolysis. An entirely new selection starting from a random DNA pool, using an all-RNA substrate and imposing the same selection pressure, also leads to RNA hydrolysis. Collectively, these results establish experimentally that small DNA sequences have the catalytic ability to direct a chemical reaction down a disfavored pathway, even when a more favorable mechanism is readily available. Our view of DNA catalysis is therefore expanded beyond merely increasing the rates of reactions that would have occurred more slowly without the catalyst.

BaM-seq and TBaM-seq, highly multiplexed and targeted RNA-seq protocols for rapid, low-cost library generation from bacterial samples.

The ability to profile transcriptomes and characterize global gene expression changes has been greatly enabled by the development of RNA sequencing technologies (RNA-seq). However, the process of generating sequencing-compatible cDNA libraries from RNA samples can be time-consuming and expensive, especially for bacterial mRNAs which lack poly(A)-tails that are often used to streamline this process for eukaryotic samples. Compared to the increasing throughput and decreasing cost of sequencing, library preparation has had limited advances. Here, we describe bacterial-multiplexed-seq (BaM-seq), an approach that enables simple barcoding of many bacterial RNA samples that decreases the time and cost of library preparation. We also present targeted-bacterial-multiplexed-seq (TBaM-seq) that allows for differential expression analysis of specific gene panels with over 100-fold enrichment in read coverage. In addition, we introduce the concept of transcriptome redistribution based on TBaM-seq that dramatically reduces the required sequencing depth while still allowing for quantification of both highly and lowly abundant transcripts. These methods accurately measure gene expression changes with high technical reproducibility and agreement with gold standard, lower throughput approaches. Together, use of these library preparation protocols allows for fast, affordable generation of sequencing libraries.

Behavioural individuality determines infection risk in clonal ant colonies.

In social groups, infection risk is not distributed evenly across individuals. Individual behaviour is a key source of variation in infection risk, yet its effects are difficult to separate from other factors (e.g., age). Here, we combine epidemiological experiments with chemical, transcriptomic, and automated behavioural analyses in clonal ant colonies, where behavioural individuality emerges among identical workers. We find that: (1) Caenorhabditis-related nematodes parasitise ant heads and affect their survival and physiology, (2) differences in infection emerge from behavioural variation alone, and reflect spatially-organised division of labour, (3) infections affect colony social organisation by causing infected workers to stay in the nest. By disproportionately infecting some workers and shifting their spatial distribution, infections reduce division of labour and increase spatial overlap between hosts, which should facilitate parasite transmission. Thus, division of labour, a defining feature of societies, not only shapes infection risk and distribution but is also modulated by parasites.

Convergent consequences of parthenogenesis on stick insect genomes.

The shift from sexual reproduction to parthenogenesis has occurred repeatedly in animals, but how the loss of sex affects genome evolution remains poorly understood. We generated reference genomes for five independently evolved parthenogenetic species in the stick insect genus Timema and their closest sexual relatives. Using these references and population genomic data, we show that parthenogenesis results in an extreme reduction of heterozygosity and often leads to genetically uniform populations. We also find evidence for less effective positive selection in parthenogenetic species, suggesting that sex is ubiquitous in natural populations because it facilitates fast rates of adaptation. Parthenogenetic species did not show increased transposable element (TE) accumulation, likely because there is little TE activity in the genus. By using replicated sexual-parthenogenetic comparisons, our study reveals how the absence of sex affects genome evolution in natural populations, providing empirical support for the negative consequences of parthenogenesis as predicted by theory.