Search

Quantifying SARS-like coronavirus (SL-CoV) evolution is critical to understanding the origins of SARS-CoV-2 and the molecular processes that could underlie future epidemic viruses. While genomic evidence implicates recombination as a factor in the emergence of SARS-CoV-2, few studies have quantified recombination rates among SL-CoVs. Here, we infer recombination rates of SL-CoVs from correlated substitutions in sequencing data using a coalescent model with recombination. Our computationally-efficient, non-phylogenetic method infers recombination parameters of both sampled sequences and the unsampled gene pools with which they recombine. We apply this approach to infer recombination parameters for a range of positive-sense RNA viruses. We then analyze a set of 191 SL-CoV sequences (including SARS-CoV-2) and find that ORF1ab and S genes frequently undergo recombination. We identify which SL-CoV sequence clusters have recombined with shared gene pools, and show that these pools have distinct structures and high recombination rates, with multiple recombination events occurring per synonymous substitution. We find that individual genes have recombined with different viral reservoirs. By decoupling contributions from mutation and recombination, we recover the phylogeny of non-recombined portions for many of these SL-CoVs, including the position of SARS-CoV-2 in this clonal phylogeny. Lastly, by analyzing 444,145 SARS-CoV-2 whole genome sequences, we show current diversity levels are insufficient to infer the within-population recombination rate of the virus since the pandemic began. Our work offers new methods for inferring recombination rates in RNA viruses with implications for understanding recombination in SARS-CoV-2 evolution and the structure of clonal relationships and gene pools shaping its origins.Significance StatementQuantifying the population genetics of SARS-like coronavirus (SL-CoV) evolution is vital to deciphering the origins of SARS-CoV-2 and pinpointing viruses with epidemic potential. While some Bayesian approaches can quantify recombination for these pathogens, the required simulations of recombination networks do not scale well with the massive amounts of sequences available in the genomics era. Our approach circumvents this by measuring correlated substitutions in sequences and fitting these data to a coalescent model with recombination. This allows us to analyze hundreds of thousands of sample sequences, and infer recombination rates for unsampled viral reservoirs. Our results provide insights into both the clonal relationships of sampled SL-CoV sequence clusters and the evolutionary dynamics of the gene pools with which they recombine.

DNA methylation in bacteria frequently serves as a simple immune system, allowing recognition of DNA from foreign sources, such as phages or selfish genetic elements. It is not well established whether methylation also frequently serves a more general epigenetic function, modifying bacterial phenotypes in a heritable manner. To address this question, here we use Oxford Nanopore sequencing to profile DNA modification marks in three natural isolates of E. coli. We first identify the DNA sequence motifs targeted by the methyltransferases in each strain. We then quantify the frequency of methylation at each of these motifs across the genome in different growth conditions. We find that motifs in specific regions of the genome consistently exhibit high or low levels of methylation. Furthermore, we show that there are replicable and consistent differences in methylated regions across different growth conditions. This suggests that during growth, E. coli transiently differentiates into distinct methylation states that depend on the growth state, raising the possibility that measuring DNA methylation alone can be used to infer bacterial growth states without additional information such as transcriptome or proteome data. These results provide new insights into the dynamics of methylation during bacterial growth, and provide evidence of differentiated cell states, a transient analogue to what is observed in the differentiation of cell types in multicellular organisms.

Bacterial cells often respond to changes in the environment by modifying protein expression. This can be achieved through changes in transcriptional or translational activity, or both. Recent research has shed some light on how natural selection shapes overall protein expression. Still, little is known about how selection acts on transcription or translation individually. To address part of this question, we implement an experimental system which allows us to measure how genetic changes affect transcription only, excluding the effects on translation. We use this system to quantify changes in three regulatory phenotypes of the lacZ promoter: transcriptional activity, plasticity, and cell-to-cell variability. We compare these phenotypes from segregating variants that have been subject to natural selection, and random variants that have never been subjected to natural selection. We show that natural selection filters out mutations causing large changes in transcriptional levels from the lacZ promoter. Further, we detect directional selection acting on transcriptional plasticity in combinations of glucose, galactose and lactose environments. Focusing on cell-to-cell variability in transcription, we describe both directional and diversifying selection acting on this phenotype depending on the environment used. We also observe a link between the phylogeny of the environmental E. coli strains and high and low transcriptional noise levels in glucose which are mediated by just one or two SNPs. Our results thus provide new insight into how one of the most well-characterized bacterial promoters is shaped in nature by selection.

Long read sequencing technologies now allow routine highly contiguous assembly of bacterial genomes. However, because of the lower accuracy of some long read data, it is often combined with short read data (e.g. Illumina), to improve assembly quality. There are a number of methods available for producing such hybrid assemblies. Here we use Illumina and Oxford Nanopore (ONT) data from 49 natural isolates of Escherichia coli to characterise differences in assembly accuracy for five assembly methods (Canu, Unicycler, Raven, Flye, and Redbean). We evaluate assembly accuracy using five metrics designed to measure structural accuracy and sequence accuracy (indel and substitution frequency). We assess structural accuracy by quantifying (1) the contiguity of chromosomes and plasmids; (2) the fraction of concordantly mapped Illumina reads withheld from the assembly; and (3) whether rRNA operons are correctly oriented. We assess indel and substitution frequency by quantifying (1) the fraction of open reading frames that appear truncated and (2) the number of variants that are called using Illumina reads only. Applying these assembly metrics to a large number of E. coli strains, we find that different assembly methods offer different advantages. In particular, we find that Unicycler assemblies have the highest sequence accuracy in non-repetitive regions, while Flye and Raven tend to be the most structurally accurate. In addition, we find that there are unidentified strain-specific characteristics that affect ONT consensus accuracy, despite individual reads having similar levels of accuracy. The differences in consensus accuracy of the ONT reads can preclude accurate assembly regardless of assembly method. These results provide quantitative insight into the best approaches for hybrid assembly of bacterial genomes and the expected levels of structural and sequence accuracy. They also show that there are intrinsic idiosyncratic strain-level differences that inhibit accurate long read bacterial genome assembly. However, we also show it is possible to diagnose problematic assemblies, even in the absence of ground truth, by comparing long-read first and short-read first assemblies.Author NotesAll supporting data, code and protocols have been provided within the article or through supplementary data files. The supporting code is available from the GitHub repository https://github.com/GeorgiaBreckell/assembly_pipeline. nine supplementary figures and three supplementary tables are available with the online version of this article.Data summarySequence data and genome assemblies for the natural isolates are available at https://www.ebi.ac.uk/ena/browser/view/PRJEB36951. Genome assemblies for additional E. coli strains used here are available from NCBI: (MG1655, SE11, REL606, CFT073, W, IA136, O157:H7-EDL933)

Most animal mitochondrial genomes are small, circular, and structurally conserved. However, recent work indicates that diverse taxa possess unusual mitochondrial genomes. In Isopoda, species in multiple lineages have atypical and rearranged mitochondrial genomes. However, more species of this speciose taxon need to be evaluated to understand the evolutionary origins of atypical mitochondrial genomes in this group. In this study, we report the presence of an atypical mitochondrial structure in the New Zealand endemic marine isopod, Isocladus armatus. Data from long and short read DNA sequencing, suggests that I. armatus has two mitochondrial chromosomes. The first chromosome consists of two mitochondrial genomes that have been inverted and fused together in a circular form, and the second chromosome consists of a single mitochondrial genome in a linearized form. This atypical mitochondrial structure has been detected in other isopod lineages, and our data from an additional divergent isopod lineage (Sphaeromatidae) lends support to the hypothesis that atypical structure evolved early in the evolution of Isopoda. Additionally, we find that a heteroplasmic site previously observed across many species within Isopoda is absent in I. armatus, but confirm the presence of two heteroplasmic sites recently reported in two other isopod species.

Bacteria often respond to dynamically changing environments by regulating gene expression. Despite this regulation being critically important for growth and survival, little is known about how selection shapes gene regulation in natural populations. To better understand the role natural selection plays in shaping bacterial gene regulation, here we compare differences in the regulatory behaviour of naturally segregating promoter variants from Escherichia coli (which have been subject to natural selection) to randomly mutated promoter variants (which have never been exposed to natural selection). We quantify gene expression phenotypes (expression level, plasticity, and noise) for hundreds of promoter variants across multiple environments, and show that segregating promoter variants are enriched for mutations with minimal effects on expression level. In many promoters, we infer that there is strong selection to maintain high levels of plasticity, and direct selection to decrease or increase cell-to-cell variability in expression. Finally, taking an integrated view, we show that across all phenotypes combined, segregating promoter variants are far more phenotypically similar than would be expected given their genetic divergence. This is the consequence of both stabilizing and directional selection acting on individual phenotypes to minimize differences among segregating variants. Taken together, these results expand our knowledge of how gene regulation is affected by natural selection and highlight the power of comparing naturally segregating polymorphisms to de novo random mutations to quantify the action of selection.

The expanding knowledge of the variety of synthetic genetic elements has enabled the construction of new and more efficient genetic circuits and yielded novel insights into molecular mechanisms. However, context dependence, in which interactions between proximal (cis) or distal (trans) elements affect the behaviour of these elements, can reduce their general applicability or predictability. Genetic insulators, which mitigate unintended context-dependent cis-interactions, have been used to address this issue. One of the most commonly used genetic insulators is a self-splicing ribozyme called RiboJ, which can be used to decouple upstream 5’ UTR in mRNA from downstream sequences (e.g., open reading frames). Despite its general use as an insulator, there has been no systematic study quantifying the efficiency of RiboJ splicing or whether this autocatalytic activity is robust to trans- and cis-genetic context. Here, we determine the robustness of RiboJ splicing in the genetic context of six widely divergent E. coli strains. We also check for possible cis-effects by assessing two SNP versions close to the catalytic site of RiboJ. We show that mRNA molecules containing RiboJ are rapidly spliced even during rapid exponential growth and high levels of gene expression, with a mean efficiency of 98%. We also show that neither the cis- nor trans-genetic context has a significant impact on RiboJ activity, suggesting this element is robust to both cis- and trans-genetic changes.

SummaryBackgroundReal-time genomic sequencing has played a major role in tracking the global spread and local transmission of SARS-CoV-2, contributing greatly to disease mitigation strategies. After effectively eliminating the virus, New Zealand experienced a second outbreak of SARS-CoV-2 in August 2020. During this August outbreak, New Zealand utilised genomic sequencing in a primary role to support its track and trace efforts for the first time, leading to a second successful elimination of the virus.MethodsWe generated the genomes of 80% of the laboratory-confirmed samples of SARS-CoV-2 from New Zealand’s August 2020 outbreak and compared these genomes to the available global genomic data.FindingsGenomic sequencing was able to rapidly identify that the new COVID-19 cases in New Zealand belonged to a single cluster and hence resulted from a single introduction. However, successful identification of the origin of this outbreak was impeded by substantial biases and gaps in global sequencing data.InterpretationAccess to a broader and more heterogenous sample of global genomic data would strengthen efforts to locate the source of any new outbreaks.FundingThis work was funded by the Ministry of Health of New Zealand, New Zealand Ministry of Business, Innovation and Employment COVID-19 Innovation Acceleration Fund (CIAF-0470), ESR Strategic Innovation Fund and the New Zealand Health Research Council (20/1018 and 20/1041).

Rapid identification of bacterial pathogens and their antimicrobial resistance (AMR) profiles is critical for minimising patient morbidity and mortality. While many sequencing methods allow deep genomic and metagenomic profiling of samples, widespread use (for example atpoint-of-care settings) is impeded because substantial sequencing and computational infrastructure is required for sequencing and analysis. Here we present Bac-PULCE (Bacterial strain and antimicrobial resistance Profiling Using Long reads via CRISPR Enrichment), which combines CRISPR-cas9 based targeted sequence enrichment with long-read sequencing. We show that this method allows simultaneous bacterial strain-level identification and antimicrobial resistance profiling of single isolates or metagenomic samples with minimal sequencing throughput. In contrast to short read sequencing, long read sequencing used in Bac-PULCE enables strain-level resolution even when targeting and sequencing highly conserved genomic regions, such as 16S rRNA. We show that these long reads allow sequencing of additional AMR genes linked to the targeted region. Additionally, long reads can be used to identify which species in a metagenomic sample harbour specific AMR loci. The ability to massively multiplex crRNAs suggests that this method has the potential to substantially increase the speed and specificity of pathogen strain identification and AMR profiling, while ensuring low computational overhead.ImportanceThere is a critical need for rapid and identification of bacterial strains and antibiotic resistance profiles in clinical settings. However, most current methods require both substantial laboratory infrastructure (e.g. for DNA sequencing), substantial compute infrastructure (e.g. for bioinformatic analyses), or both. Here we present a new method, Bac-PULCE, (Bacterial strain and antimicrobial resistance Profiling Using Long reads via CRISPR Enrichment), which combines CRISPR-cas9 based targeted sequence enrichment with long-read sequencing on the Oxford Nanopore platform. This allows rapid profiling of bacterial strains and antibiotic resistance genes in a sample while requiring very little laboratory or computational infrastructure.

Rapid and cost-efficient whole-genome sequencing of SARS-CoV-2, the virus that causes COVID-19, is critical for understanding viral transmission dynamics. Here we show that using a new multiplexed set of primers in conjunction with the Oxford Nanopore Rapid Barcode library kit allows for faster, simpler, and less expensive SARS-CoV-2 genome sequencing. This primer set results in amplicons that exhibit lower levels of variation in coverage compared to other commonly used primer sets. Using five SARS-CoV-2 patient samples with Cq values between 20 and 31, we show that high-quality genomes can be generated with as few as 10,000 reads (approximately 5 Mbp of sequence data). We also show that mis-classification of barcodes, which may be more likely when using the Oxford Nanopore Rapid Barcode library prep, is unlikely to cause problems in variant calling. This method reduces the time from RNA to genome sequence by more than half compared to the more standard ligation-based Oxford Nanopore library preparation method at considerably lower costs.

Marine organisms generally exhibit one of two developmental modes: biphasic, with distinct adult and larval morphology, and direct development, in which larvae resemble adults. Developmental mode is thought to significantly influence dispersal, with direct developers expected to have much lower dispersal potential. However, in contrast to our relatively good understanding of dispersal and population connectivity for biphasic species, comparatively little is known about direct developers. In this study, we use a panel of 8,020 SNPs to investigate population structure and gene flow for a direct developing species, the New Zealand endemic marine isopodIsocladus armatus. On a small spatial scale (20 kms), gene flow between locations is extremely high and suggests an island model of migration. However, over larger spatial scales (600km), populations exhibit a clear pattern of isolation-by-distance. Because our sampling range is intersected by two well-known biogeographic barriers (the East Cape and the Cook Strait), our study provides an opportunity to understand how such barriers influence dispersal in direct developers. Our results indicate thatI. armatusexhibits significant migration across these barriers, and suggests that ocean currents associated with these locations do not present a barrier to dispersal. Interestingly, we do find evidence of a north-south population genetic break occurring between Māhia and Wellington, two locations where there are no obvious biogeographic barriers between them. We conclude that developmental life history largely predicts dispersal in intertidal marine isopods. However, localised biogeographic processes can disrupt this expectation.

BackgroundUsing metagenomics to determine animal diet offers a new and promising alternative to current methods. Here we show that rapid and inexpensive diet quantification is possible through metagenomic sequencing with the portable Oxford Nanopore Technologies (ONT) MinION. Using an amplification-free approach, we profiled the stomach contents from wild-caught rats.ResultsWe conservatively identified diet items from over 50 taxonomic orders, ranging across nine phyla that include plants, vertebrates, invertebrates, and fungi. This highlights the wide range of taxa that can be identified using this simple approach. We calibrate the accuracy of this method by comparing the characteristics of reads matching the ground-truth host genome (rat) to those matching diet items, and show that at the family-level, false positive taxon assignments are approximately 97.5% accurate. We also suggest a way to mitigate for database biases in metagenomic approaches. Finally, we implement a constrained ordination analysis and show that we can identify the sampling location of an individual rat within tens of kilometres based on diet content alone.ConclusionsThis work establishes proof-of-principle for long-read metagenomic methods in quantitative diet analysis. We show that diet content can be quantified even with limited expertise, using a simple, amplification free workflow and a relatively inexpensive and accessible next generation sequencing method. Continued increases in the accuracy and throughput of ONT sequencing, along with improved genomic databases, suggests that a metagenomic approach to quantification of animal diets will become an important method in the future.

Bacteria adapt to changes in their environment by regulating gene expression, often at the level of transcription. However, since the molecular processes underlying gene regulation are subject to thermodynamic and other stochastic fluctuations, gene expression is inherently noisy, and identical cells in a homogeneous environment can display highly heterogeneous expression levels. To study how stochasticity affects gene regulation at the single-cell level, it is crucial to be able to directly follow gene expression dynamics in single cells under changing environmental conditions. Recently developed microfluidic devices, used in combination with quantitative fluorescence time-lapse microscopy, represent a highly promising experimental approach, allowing tracking of lineages of single cells over long time-scales while simultaneously measuring their growth and gene expression. However, current devices do not allow controlled dynamical changes to the environmental conditions which are needed to study gene regulation. In addition, automated analysis of the imaging data from such devices is still highly challenging and no standard software is currently available. To address these challenges, we here present an integrated experimental and computational setup featuring, on the one hand, a new dual-input microfluidic chip which allows mixing and switching between two growth media and, on the other hand, a novel image analysis software which jointly optimizes segmentation and tracking of the cells and allows interactive user-guided fine-tuning of its results. To demonstrate the power of our approach, we study the lac operon regulation in E. coli cells grown in an environment that switches between glucose and lactose, and quantify stochastic lag times and memory at the single cell level.

Although it is often tacitly assumed that gene regulatory interactions are finely tuned, how accurate gene regulation could evolve from a state without regulation is unclear. Moreover, gene expression noise would seem to impede the evolution of accurate gene regulation, and previous investigations have provided circumstantial evidence that natural selection has acted to lower noise levels. By evolving synthetic Escherichia coli promoters de novo, we here show that, contrary to expectations, promoters exhibit low noise by default. Instead, selection must have acted to increase the noise levels of highly regulated E. coli promoters. We present a general theory of the interplay between gene expression noise and gene regulation that explains these observations. The theory shows that propagation of expression noise from regulators to their targets is not an unwanted side-effect of regulation, but rather acts as a rudimentary form of regulation that facilitates the evolution of more accurate regulation. DOI: http://dx.doi.org/10.7554/eLife.05856.001, eLife digest Genes are stretches of DNA that contain the instructions needed to make proteins and other molecules. By changing how much protein is produced from each gene (i.e., its expression), many organisms—including humans—can produce a wide variety of cell types with very different behaviors. Similarly, single-celled organisms, such as bacteria, can adapt to survive and grow in different environments by changing gene expression levels. It is thus thought that gene expression must be precisely controlled. However, the molecular processes involved in gene expression are subject to random fluctuations, and so gene expression is inherently ‘noisy’. This means that even groups of identical cells in identical environments will show variation in their gene expression patterns. Furthermore, different genes show different levels of noise. The DNA sequence of a part of each gene, called the promoter, has a big effect on these noise levels. Consequently, gene expression noise is a genetically encoded trait, and can therefore be shaped by natural selection. But it remains largely unclear how natural selection has affected gene expression noise. Now, Wolf et al. have carefully measured the gene expression noise of hundreds of synthetic promoters that were evolved in the laboratory from random DNA sequences, and a similar number of natural promoters in a bacterium called E. coli. These experiments revealed that, contrary to expectation, most lab-evolved promoters had low levels of noise. On the other hand, many natural promoters had high levels of noise. Wolf et al. also found that noisy promoters tend to be highly regulated by transcription factors: the proteins that control gene expression by binding to promoter regions. Together, these results imply that unregulated promoters start by having low noise as a default state. Selection pressures must then have caused some E. coli promoters to become regulated by transcription factors and raise their noise levels. But, what might these selection pressures have been? Many genes need to be expressed at different levels in different conditions, and it is generally accepted that regulation by transcription factors evolves to ‘satisfy’ these requirements. However, transcription factors are themselves noisy, and this noise necessarily propagates to their target genes. Wolf et al. have now developed a general theory showing that this noise-propagation can often benefit an organism. This explains why natural selection can favor an increase in noise levels for regulated genes. Importantly, by showing that the main role of a transcription factor can be to increase the noise of its targets, it suddenly becomes very easy to see how new gene regulatory interactions can evolve from scratch. The next steps in understanding of how gene expression noise evolves will involve manipulating the expression noise of a gene, and measuring how selection acts on such changes. DOI: http://dx.doi.org/10.7554/eLife.05856.002