Summary: Genome-wide expression analysis can result in large numbers of clusters of co-expressed genes. Although there are tools for ab initio discovery of transcription factor-binding sites, most do not provide a quick and easy way to study large numbers of clusters. To address this, we introduce a web tool called MEME-LaB. The tool wraps MEME (an ab initio motif finder), providing an interface for users to input multiple gene clusters, retrieve promoter sequences, run motif finding and then easily browse and condense the results, facilitating better interpretation of the results from large-scale datasets. Availability: MEME-LaB is freely accessible at: http://wsbc.warwick.ac.uk/wsbcToolsWebpage/. Contact: p.e.brown@warwick.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Root cell types possess distinct immunity gene networks that are linked to cell identity networks, as revealed by cell type-specific RNA-seq and a paired motif

Background: Robustness is a central property of living systems, enabling function to be maintained against environmental perturbations. A key challenge is to identify the structures in biological circuits that confer system-level properties such as robustness. Circadian clocks allow organisms to adapt to the predictable changes of the 24-hour day/night cycle by generating endogenous rhythms that can be entrained to the external cycle. In all organisms, the clock circuits typically comprise multiple interlocked feedback loops controlling the rhythmic expression of key genes. Previously, we showed that such architectures increase the flexibility of the clock's rhythmic behaviour. We now test the relationship between flexibility and robustness, using a mathematical model of the circuit controlling conidiation in the fungus Neurospora crassa. Results: The circuit modelled in this work consists of a central negative feedback loop, in which the frequency (frq) gene inhibits its transcriptional activator white collar-1 (wc-1), interlocked with a positive feedback loop in which FRQ protein upregulates WC-1 production. Importantly, our model reproduces the observed entrainment of this circuit under light/dark cycles with varying photoperiod and cycle duration. Our simulations show that whilst the level of frq mRNA is driven directly by the light input, the falling phase of FRQ protein, a molecular correlate of conidiation, maintains a constant phase that is uncoupled from the times of dawn and dusk. The model predicts the behaviour of mutants that uncouple WC-1 production from FRQ's positive feedback, and shows that the positive loop enhances the buffering of conidiation phase against seasonal photoperiod changes. This property is quantified using Kitano's measure for the overall robustness of a regulated system output. Further analysis demonstrates that this functional robustness is a consequence of the greater evolutionary flexibility conferred on the circuit by the interlocking loop structure. Conclusions: Our model shows that the behaviour of the fungal clock in light-dark cycles can be accounted for by a transcription-translation feedback model of the central FRQ-WC oscillator. More generally, we provide an example of a biological circuit in which greater flexibility yields improved robustness, while also introducing novel sensitivity analysis techniques applicable to a broader range of cellular oscillators.

Molecular mimicry is one of the leading mechanisms by which infectious agents can induce autoimmunity. Whether a similar mechanism triggers an antitumor immune response is unexplored, and the role of antiviral T cells infiltrating the tumor has remained anecdotal. To address these questions, we first developed a bioinformatic tool to identify tumor peptides with high similarity to viral epitopes. Using peptides identified by this tool, we demonstrated that, in mice, preexisting immunity toward specific viral epitopes enhanced the efficacy of cancer immunotherapy via molecular mimicry in different settings. To understand whether this mechanism could partly explain immunotherapy responsiveness in humans, we analyzed a cohort of patients with melanoma undergoing anti-PD1 treatment who had a high IgG titer for cytomegalovirus (CMV). In this cohort of patients, we showed that high levels of CMV-specific antibodies were associated with prolonged progression-free survival and found that, in some cases, peripheral blood mononuclear cells (PBMC) could cross-react with both melanoma and CMV homologous peptides. Finally, T-cell receptor sequencing revealed expansion of the same CD8+ T-cell clones when PBMCs were expanded with tumor or homologous viral peptides. In conclusion, we have demonstrated that preexisting immunity and molecular mimicry could influence the response to immunotherapies. In addition, we have developed a free online tool that can identify tumor antigens and neoantigens highly similar to pathogen antigens to exploit molecular mimicry and cross-reactive T cells in cancer vaccine development.

Gambiense human African trypanosomiasis (gHAT) is a virulent disease declining in burden but still endemic in West and Central Africa. Although it is targeted for elimination of transmission by 2030, there remain numerous questions about the drivers of infection and how these vary geographically. In this study we focus on the Democratic Republic of Congo (DRC), which accounted for 84% of the global case burden in 2016, to explore changes in transmission across the country and elucidate factors which may have contributed to the persistence of disease or success of interventions in different regions. We present a Bayesian fitting methodology, applied to 168 endemic health zones (∼100,000 population size), which allows for calibration of a mechanistic gHAT model to case data (from the World Health Organization HAT Atlas) in an adaptive and automated framework. It was found that the model needed to capture improvements in passive detection to match observed trends in the data within former Bandundu and Bas Congo provinces indicating these regions have substantially reduced time to detection. Health zones in these provinces generally had longer burn-in periods during fitting due to additional model parameters. Posterior probability distributions were found for a range of fitted parameters in each health zone; these included the basic reproduction number estimates for pre-1998 (R₀) which was inferred to be between 1 and 1.14, in line with previous gHAT estimates, with higher median values typically in health zones with more case reporting in the 2000s. Previously, it was not clear whether a fall in active case finding in the period contributed to the declining case numbers. The modelling here accounts for variable screening and suggests that underlying transmission has also reduced greatly—on average 96% in former Equateur, 93% in former Bas Congo and 89% in former Bandundu—Equateur and Bandundu having had the highest case burdens in 2000. This analysis also sets out a framework to enable future predictions for the country.

Background: In recent years, a programme of vector control, screening and treatment of gambiense human African trypanosomiasis (gHAT) infections led to a rapid decline in cases in the Mandoul focus of Chad. To represent the biology of transmission between humans and tsetse, we previously developed a mechanistic transmission model, fitted to data between 2000 and 2013 which suggested that transmission was interrupted by 2015. The present study outlines refinements to the model to: (1) Assess whether elimination of transmission has already been achieved despite low-level case reporting; (2) quantify the role of intensified interventions in transmission reduction; and (3) predict the trajectory of gHAT in Mandoul for the next decade under different strategies. Method: Our previous gHAT transmission model for Mandoul was updated using human case data (2000–2019) and a series of model refinements. These include how diagnostic specificity is incorporated into the model and improvements to the fitting method (increased variance in observed case reporting and how underreporting and improvements to passive screening are captured). A side-by-side comparison of fitting to case data was performed between the models. Results: We estimated that passive detection rates have increased due to improvements in diagnostic availability in fixed health facilities since 2015, by 2.1-fold for stage 1 detection, and 1.5-fold for stage 2. We find that whilst the diagnostic algorithm for active screening is estimated to be highly specific (95% credible interval (CI) 99.9–100%, Specificity = 99.9%), the high screening and low infection levels mean that some recently reported cases with no parasitological confirmation might be false positives. We also find that the focus-wide tsetse reduction estimated through model fitting (95% CI 96.1–99.6%, Reduction = 99.1%) is comparable to the reduction previously measured by the decline in tsetse catches from monitoring traps. In line with previous results, the model suggests that transmission was interrupted in 2015 due to intensified interventions. Conclusions: We recommend that additional confirmatory testing is performed in Mandoul to ensure the endgame can be carefully monitored. More specific measurement of cases, would better inform when it is safe to stop active screening and vector control, provided there is a strong passive surveillance system in place. Graphical