Search

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

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.

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.