Correlated substitutions reveal SARS-like coronaviruses recombine frequently with a diverse set of structured gene pools.

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 analyses suggest recombination was 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 >400,000 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.

Use of Genomics to Track Coronavirus Disease Outbreaks, New Zealand.

Real-time genomic sequencing has played a major role in tracking the global spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), contributing greatly to disease mitigation strategies. In August 2020, after having eliminated the virus, New Zealand experienced a second outbreak. During that outbreak, New Zealand used genomic sequencing in a primary role, leading to a second elimination of the virus. We generated genomes from 78% of the laboratory-confirmed samples of SARS-CoV-2 from the second outbreak and compared them with the available global genomic data. Genomic sequencing rapidly identified that virus causing the second outbreak in New Zealand belonged to a single cluster, thus resulting from a single introduction. However, successful identification of the origin of this outbreak was impeded by substantial biases and gaps in global sequencing data. Access to a broader and more heterogenous sample of global genomic data would strengthen efforts to locate the source of any new outbreaks.

Testing the advantages and disadvantages of short- and long- read eukaryotic metagenomics using simulated reads.

BACKGROUND: The first step in understanding ecological community diversity and dynamics is quantifying community membership. An increasingly common method for doing so is through metagenomics. Because of the rapidly increasing popularity of this approach, a large number of computational tools and pipelines are available for analysing metagenomic data. However, the majority of these tools have been designed and benchmarked using highly accurate short read data (i.e. Illumina), with few studies benchmarking classification accuracy for long error-prone reads (PacBio or Oxford Nanopore). In addition, few tools have been benchmarked for non-microbial communities. RESULTS: Here we compare simulated long reads from Oxford Nanopore and Pacific Biosciences (PacBio) with high accuracy Illumina read sets to systematically investigate the effects of sequence length and taxon type on classification accuracy for metagenomic data from both microbial and non-microbial communities. We show that very generally, classification accuracy is far lower for non-microbial communities, even at low taxonomic resolution (e.g. family rather than genus). We then show that for two popular taxonomic classifiers, long reads can significantly increase classification accuracy, and this is most pronounced for non-microbial communities. CONCLUSIONS: This work provides insight on the expected accuracy for metagenomic analyses for different taxonomic groups, and establishes the point at which read length becomes more important than error rate for assigning the correct taxon.

The predictability of molecular evolution during functional innovation.

Determining the molecular changes that give rise to functional innovations is a major unresolved problem in biology. The paucity of examples has served as a significant hindrance in furthering our understanding of this process. Here we used experimental evolution with the bacterium Escherichia coli to quantify the molecular changes underlying functional innovation in 68 independent instances ranging over 22 different metabolic functions. Using whole-genome sequencing, we show that the relative contribution of regulatory and structural mutations depends on the cellular context of the metabolic function. In addition, we find that regulatory mutations affect genes that act in pathways relevant to the novel function, whereas structural mutations affect genes that act in unrelated pathways. Finally, we use population genetic modeling to show that the relative contributions of regulatory and structural mutations during functional innovation may be affected by population size. These results provide a predictive framework for the molecular basis of evolutionary innovation, which is essential for anticipating future evolutionary trajectories in the face of rapid environmental change.

Combining Shigella Tn-seq data with gold-standard E. coli gene deletion data suggests rare transitions between essential and non-essential gene functionality.

BACKGROUND: Gene essentiality - whether or not a gene is necessary for cell growth - is a fundamental component of gene function. It is not well established how quickly gene essentiality can change, as few studies have compared empirical measures of essentiality between closely related organisms. RESULTS: Here we present the results of a Tn-seq experiment designed to detect essential protein coding genes in the bacterial pathogen Shigella flexneri 2a 2457T on a genome-wide scale. Superficial analysis of this data suggested that 481 protein-coding genes in this Shigella strain are critical for robust cellular growth on rich media. Comparison of this set of genes with a gold-standard data set of essential genes in the closely related Escherichia coli K12 BW25113 revealed that an excessive number of genes appeared essential in Shigella but non-essential in E. coli. Importantly, and in converse to this comparison, we found no genes that were essential in E. coli and non-essential in Shigella, implying that many genes were artefactually inferred as essential in Shigella. Controlling for such artefacts resulted in a much smaller set of discrepant genes. Among these, we identified three sets of functionally related genes, two of which have previously been implicated as critical for Shigella growth, but which are dispensable for E. coli growth. CONCLUSIONS: The data presented here highlight the small number of protein coding genes for which we have strong evidence that their essentiality status differs between the closely related bacterial taxa E. coli and Shigella. A set of genes involved in acetate utilization provides a canonical example. These results leave open the possibility of developing strain-specific antibiotic treatments targeting such differentially essential genes, but suggest that such opportunities may be rare in closely related bacteria.

Automated reconstruction of whole-genome phylogenies from short-sequence reads.

Studies of microbial evolutionary dynamics are being transformed by the availability of affordable high-throughput sequencing technologies, which allow whole-genome sequencing of hundreds of related taxa in a single study. Reconstructing a phylogenetic tree of these taxa is generally a crucial step in any evolutionary analysis. Instead of constructing genome assemblies for all taxa, annotating these assemblies, and aligning orthologous genes, many recent studies 1) directly map raw sequencing reads to a single reference sequence, 2) extract single nucleotide polymorphisms (SNPs), and 3) infer the phylogenetic tree using maximum likelihood methods from the aligned SNP positions. However, here we show that, when using such methods to reconstruct phylogenies from sets of simulated sequences, both the exclusion of nonpolymorphic positions and the alignment to a single reference genome, introduce systematic biases and errors in phylogeny reconstruction. To address these problems, we developed a new method that combines alignments from mappings to multiple reference sequences and show that this successfully removes biases from the reconstructed phylogenies. We implemented this method as a web server named REALPHY (Reference sequence Alignment-based Phylogeny builder), which fully automates phylogenetic reconstruction from raw sequencing reads.

Patterns of evolutionary conservation of essential genes correlate with their compensability.

Essential genes code for fundamental cellular functions required for the viability of an organism. For this reason, essential genes are often highly conserved across organisms. However, this is not always the case: orthologues of genes that are essential in one organism are sometimes not essential in other organisms or are absent from their genomes. This suggests that, in the course of evolution, essential genes can be rendered nonessential. How can a gene become non-essential? Here we used genetic manipulation to deplete the products of 26 different essential genes in Escherichia coli. This depletion results in a lethal phenotype, which could often be rescued by the overexpression of a non-homologous, non-essential gene, most likely through replacement of the essential function. We also show that, in a smaller number of cases, the essential genes can be fully deleted from the genome, suggesting that complete functional replacement is possible. Finally, we show that essential genes whose function can be replaced in the laboratory are more likely to be non-essential or not present in other taxa. These results are consistent with the notion that patterns of evolutionary conservation of essential genes are influenced by their compensability-that is, by how easily they can be functionally replaced, for example through increased expression of other genes.

A genome-wide analysis of promoter-mediated phenotypic noise in Escherichia coli.

Gene expression is subject to random perturbations that lead to fluctuations in the rate of protein production. As a consequence, for any given protein, genetically identical organisms living in a constant environment will contain different amounts of that particular protein, resulting in different phenotypes. This phenomenon is known as "phenotypic noise." In bacterial systems, previous studies have shown that, for specific genes, both transcriptional and translational processes affect phenotypic noise. Here, we focus on how the promoter regions of genes affect noise and ask whether levels of promoter-mediated noise are correlated with genes' functional attributes, using data for over 60% of all promoters in Escherichia coli. We find that essential genes and genes with a high degree of evolutionary conservation have promoters that confer low levels of noise. We also find that the level of noise cannot be attributed to the evolutionary time that different genes have spent in the genome of E. coli. In contrast to previous results in eukaryotes, we find no association between promoter-mediated noise and gene expression plasticity. These results are consistent with the hypothesis that, in bacteria, natural selection can act to reduce gene expression noise and that some of this noise is controlled through the sequence of the promoter region alone.

Opposite effects of KCTD subunit domains on GABA(B) receptor-mediated desensitization.

GABA(B) receptors assemble from principle and auxiliary subunits. The principle subunits GABA(B1) and GABA(B2) form functional heteromeric GABA(B(1,2)) receptors that associate with homotetramers of auxiliary KCTD8, -12, -12b, or -16 (named after their K(+) channel tetramerization domain) subunits. These auxiliary subunits constitute receptor subtypes with distinct functional properties. KCTD12 and -12b generate desensitizing receptor responses while KCTD8 and -16 generate largely non-desensitizing receptor responses. The structural elements of the KCTDs underlying these differences in desensitization are unknown. KCTDs are modular proteins comprising a T1 tetramerization domain, which binds to GABA(B2), and a H1 homology domain. KCTD8 and -16 contain an additional C-terminal H2 homology domain that is not sequence-related to the H1 domains. No functions are known for the H1 and H2 domains. Here we addressed which domains and sequence motifs in KCTD proteins regulate desensitization of the receptor response. We found that the H1 domains in KCTD12 and -12b mediate desensitization through a particular sequence motif, T/NFLEQ, which is not present in the H1 domains of KCTD8 and -16. In addition, the H2 domains in KCTD8 and -16 inhibit desensitization when expressed C-terminal to the H1 domains but not when expressed as a separate protein in trans. Intriguingly, the inhibitory effect of the H2 domain is sequence-independent, suggesting that the H2 domain sterically hinders desensitization by the H1 domain. Evolutionary analysis supports that KCTD12 and -12b evolved desensitizing properties by liberating their H1 domains from antagonistic H2 domains and acquisition of the T/NFLEQ motif.

Quantitative analysis of persister fractions suggests different mechanisms of formation among environmental isolates of E. coli.

BACKGROUND: Bacterial persistence describes a phenomenon wherein a small subpopulation of cells is able to survive a challenge with high doses of an antibiotic (or other stressor) better than the majority of the population. Previous work has shown that cells that are in a dormant or slow-growing state are persistent to antibiotic treatment and that populations with higher fractions of dormant cells exhibit higher levels of persistence. These data suggest that a major determinant of the fraction of persisters within a population is the rate at which cells enter and exit from dormancy. However, it is not known whether there are physiological changes in addition to dormancy that influence persistence. Here, we use quantitative measurements of persister fractions in a set of environmental isolates of E. coli together with a mathematical model of persister formation to test whether a single general physiological change, such as cell dormancy, can explain the differences in persister phenotypes observed in different strains. RESULTS: If a single physiological change (e.g. cell dormancy) underlies most persister phenotypes, then strains should exhibit characteristic fractions of persister cells: some strains will consistently have high fractions of persisters (dormant cells), whereas others will have low fractions. Although we found substantial variation in the fraction of persisters between different environmental isolates of E. coli, these fractions were not correlated across antibiotics. Some strains exhibited high persister fractions in one antibiotic, but low persister fractions in a second antibiotic. Surprisingly, no correlation in persister fractions was observed between any two drugs, even for antibiotics with nearly identical modes of action (ciprofloxacin and nalidixic acid). CONCLUSIONS: These data support the hypothesis that there is no single physiological change that determines the persistence level in a population of cells. Instead, the fraction of cells that survive antibiotic treatment (persist) depends critically on the specific antibiotic that is used, suggesting that physiological changes in addition to dormancy can underlie persister phenotypes.

Growth condition-dependent differences in methylation imply transiently differentiated DNA methylation states in Escherichia coli.

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. However, DNA methylation also affects other cell phenotypes in a heritable manner (i.e. epigenetically). While there are several examples of methylation affecting transcription in an epigenetic manner in highly localized contexts, it is not well-established how frequently methylation serves a more general epigenetic function over larger genomic scales. To address this question, here we use Oxford Nanopore sequencing to profile DNA modification marks in three natural isolates of Escherichia 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 entire 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 differentiate 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 show the utility of using Oxford Nanopore sequencing as an economic means to infer DNA methylation status. They also provide new insights into the dynamics of methylation during bacterial growth and provide evidence of differentiated cell states, a transient analog to what is observed in the differentiation of cell types in multicellular organisms.

Human gene and microbial analyses in rectal cancer complete responses to radiotherapy.

BACKGROUND: The gold standard treatment for locally advanced rectal cancer is total mesorectal excision after preoperative chemoradiotherapy. Response to chemoradiotherapy varies, with some patients completely responding to the treatment and some failing to respond at all. Identifying biomarkers of response to chemoradiotherapy could allow patients to avoid unnecessary treatment-associated morbidity rate. While previous studies have attempted to identify such biomarkers, none have reached clinical utility, which may be due to heterogeneity of the cancer. In this study, potential human gene and microbial biomarkers were explored in a cohort of rectal cancer patients who underwent chemoradiotherapy. METHODS: RNA sequencing was carried out on matched tumour and adjacent normal rectum biopsies from patients with rectal cancer with varying chemoradiotherapy responses treated between 2016 and 2019 at two institutions. Enriched genes and microbes from tumours of complete responders were compared with those from tumours of others with lesser response. RESULTS: In 39 patients analysed, enriched gene sets in complete responders indicate involvement of immune responses, including immunoglobulin production, B cell activation and response to bacteria (adjusted P values <0.050). Bacteria such as Ruminococcaceae bacterium and Bacteroides thetaiotaomicron were documented to be abundant in tumours of complete responders compared with all other patients (adjusted P value <0.100). CONCLUSION: These results identify potential genetic and microbial biomarkers of response to chemoradiotherapy in rectal cancer, as well as suggesting a potential mechanism of complete response to chemoradiotherapy that may benefit further testing in the laboratory.

Gene regulation in Escherichia coli is commonly selected for both high plasticity and low noise.

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. 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.

Long-read sequencing reveals atypical mitochondrial genome structure in a New Zealand marine isopod.

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 suggest 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 an asymmetrical site previously observed across many species within Isopoda is absent in I. armatus, but confirm the presence of two asymmetrical sites recently reported in two other isopod species.

A simple screen to identify promoters conferring high levels of phenotypic noise.

Genetically identical populations of unicellular organisms often show marked variation in some phenotypic traits. To investigate the molecular causes and possible biological functions of this phenotypic noise, it would be useful to have a method to identify genes whose expression varies stochastically on a certain time scale. Here, we developed such a method and used it for identifying genes with high levels of phenotypic noise in Salmonella enterica ssp. I serovar Typhimurium (S. Typhimurium). We created a genomic plasmid library fused to a green fluorescent protein (GFP) reporter and subjected replicate populations harboring this library to fluctuating selection for GFP expression using fluorescent-activated cell sorting (FACS). After seven rounds of fluctuating selection, the populations were strongly enriched for promoters that showed a high amount of noise in gene expression. Our results indicate that the activity of some promoters of S. Typhimurium varies on such a short time scale that these promoters can absorb rapid fluctuations in the direction of selection, as imposed during our experiment. The genomic fragments that conferred the highest levels of phenotypic variation were promoters controlling the synthesis of flagella, which are associated with virulence and host-pathogen interactions. This confirms earlier reports that phenotypic noise may play a role in pathogenesis and indicates that these promoters have among the highest levels of noise in the S. Typhimurium genome. This approach can be applied to many other bacterial and eukaryotic systems as a simple method for identifying genes with noisy expression.

Efficiency of the synthetic self-splicing RiboJ ribozyme is robust to cis- and trans-changes in genetic background.

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 cis- or trans-genetic elements affect the behavior 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.

Understanding the evolutionary fate of finite populations: the dynamics of mutational effects.

The most consistent result in more than two decades of experimental evolution is that the fitness of populations adapting to a constant environment does not increase indefinitely, but reaches a plateau. Using experimental evolution with bacteriophage, we show here that the converse is also true. In populations small enough such that drift overwhelms selection and causes fitness to decrease, fitness declines down to a plateau. We demonstrate theoretically that both of these phenomena must be due either to changes in the ratio of beneficial to deleterious mutations, the size of mutational effects, or both. We use mutation accumulation experiments and molecular data from experimental evolution to show that the most significant change in mutational effects is a drastic increase in the rate of beneficial mutation as fitness decreases. In contrast, the size of mutational effects changes little even as organismal fitness changes over several orders of magnitude. These findings have significant implications for the dynamics of adaptation.

Complete Genome Sequences of 47 Environmental Isolates of Escherichia coli.

Escherichia coli is commonly considered a host-associated bacterium. However, there is evidence that some strains occupy environmental (non-host-associated) niches. Here, we report the complete genomes of 47 Escherichia coli environmental isolates. These will be useful for understanding the dynamics of plasmids, phages, and other repetitive genetic elements.

Rapid and inexpensive whole-genome sequencing of SARS-CoV-2 using 1200 bp tiled amplicons and Oxford Nanopore Rapid Barcoding.

Rapid and cost-efficient whole-genome sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes coronavirus disease 2019, 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 (∼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.

Concordant geographic and genetic structure revealed by genotyping-by-sequencing in a New Zealand marine isopod.

Population genetic structure in the marine environment can be influenced by life-history traits such as developmental mode (biphasic, with distinct adult and larval morphology, and direct development, in which larvae resemble adults) or habitat specificity, as well as geography and selection. Developmental mode is thought to significantly influence dispersal, with direct developers expected to have much lower dispersal potential. However, this prediction can be complicated by the presence of geophysical barriers to dispersal. In this study, we use a panel of 8,020 SNPs to investigate population structure and biogeography over multiple spatial scales for a direct-developing species, the New Zealand endemic marine isopod Isocladus armatus. 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. On a small spatial scale (20 km), gene flow between locations is extremely high, suggestive of an island model of migration. However, over larger spatial scales (600 km), populations exhibit a clear pattern of isolation-by-distance. Our results indicate that I. armatus exhibits significant migration across the hypothesized barriers and suggest that large-scale ocean currents associated with these locations do not present a barrier to dispersal. Interestingly, we find evidence of a north-south population genetic break occurring between Mahia and Wellington. While no known geophysical barrier is apparent in this area, it coincides with the location of a proposed border between bioregions. Analysis of loci under selection revealed that both isolation-by-distance and adaption may be contributing to the degree of population structure we have observed here. We conclude that developmental life history largely predicts dispersal in the intertidal isopod I. armatus. However, localized biogeographic processes can disrupt this expectation, and this may explain the potential meta-population detected in the Auckland region.