Origin of ultradian pulsatility in the hypothalamic–pituitary–adrenal axis

The hypothalamic–pituitary–adrenal (HPA) axis is a neuroendocrine system that regulates the circulating levels of vital glucocorticoid hormones. The activity of the HPA axis is characterized not only by a classic circadian rhythm, but also by an ultradian pattern of discrete pulsatile release of glucocorticoids. A number of psychiatric and metabolic diseases are associated with changes in glucocorticoid pulsatility, and it is now clear that glucocorticoid responsive genes respond to these rapid fluctuations in a biologically meaningful way. Theoretical modelling has enabled us to identify and explore potential mechanisms underlying the ultradian activity in this axis, which to date have not been identified successfully. We demonstrate that the combination of delay with feed-forward and feedback loops in the pituitary–adrenal system is sufficient to give rise to ultradian pulsatility in the absence of an ultradian source from a supra-pituitary site. Moreover, our model enables us to predict the different patterns of glucocorticoid release mediated by changes in hypophysial-portal corticotrophin-releasing hormone levels, with results that parallel our experimental in vivo data

Comparison of Shiga toxin-encoding bacteriophages in highly pathogenic strains of Shiga toxin-producing Escherichia coli O157:H7 in the UK

Over the last 35 years in the UK, the burden of Shiga toxin-producing Escherichia coli (STEC) O157:H7 infection has, during different periods of time, been associated with five different sub-lineages (1983-1995, Ia, I/IIa and I/IIb; 1996-2014, Ic; and 2015-2018, IIb). The acquisition of a stx2a-encoding bacteriophage by these five sub-lineages appears to have coincided with their respective emergences. The Oxford Nanopore Technologies (ONT) system was used to sequence, characterize and compare the stx-encoding prophages harboured by each sub-lineage to investigate the integration of this key virulence factor. The stx2a-encoding prophages from each of the lineages causing clinical disease in the UK were all different, including the two UK sub-lineages (Ia and I/IIa) circulating concurrently and causing severe disease in the early 1980s. Comparisons between the stx2a-encoding prophage in sub-lineages I/IIb and IIb revealed similarity to the prophage commonly found to encode stx2c, and the same site of bacteriophage integration (sbcB) as stx2c-encoding prophage. These data suggest independent acquisition of previously unobserved stx2a-encoding phage is more likely to have contributed to the emergence of STEC O157:H7 sub-lineages in the UK than intra-UK lineage to lineage phage transmission. In contrast, the stx2c-encoding prophage showed a high level of similarity across lineages and time, consistent with the model of stx2c being present in the common ancestor to extant STEC O157:H7 and maintained by vertical inheritance in the majority of the population. Studying the nature of the stx-encoding bacteriophage contributes to our understanding of the emergence of highly pathogenic strains of STEC O157:H7

Analysis of a small outbreak of Shiga toxin-producing Escherichia coli O157:H7 using long-read sequencing

Compared to short-read sequencing data, long-read sequencing facilitates single contiguous de novo assemblies and characterization of the prophage region of the genome. Here, we describe our methodological approach to using Oxford Nanopore Technology (ONT) sequencing data to quantify genetic relatedness and to look for microevolutionary events in the core and accessory genomes to assess the within-outbreak variation of four genetically and epidemiologically linked isolates. Analysis of both Illumina and ONT sequencing data detected one SNP between the four sequences of the outbreak isolates. The variant calling procedure highlighted the importance of masking homologous sequences in the reference genome regardless of the sequencing technology used. Variant calling also highlighted the systemic errors in ONT base-calling and ambiguous mapping of Illumina reads that results in variations in the genetic distance when comparing one technology to the other. The prophage component of the outbreak strain was analysed, and nine of the 16 prophages showed some similarity to the prophage in the Sakai reference genome, including the stx2a-encoding phage. Prophage comparison between the outbreak isolates identified minor genome rearrangements in one of the isolates, including an inversion and a deletion event. The ability to characterize the accessory genome in this way is the first step to understanding the significance of these microevolutionary events and their impact on the evolutionary history, virulence and potentially the likely source and transmission of this zoonotic, foodborne pathogen

Re-analysis of an outbreak of Shiga toxin-producing Escherichia coli O157:H7 associated with raw drinking milk using Nanopore sequencing

The aim of this study was to compare Illumina and Oxford Nanopore Technology (ONT) sequencing data to quantify genetic variation to assess within-outbreak strain relatedness and characterise microevolutionary events in the accessory genomes of a cluster of 23 genetically and epidemiologically linked isolates related to an outbreak of Shiga toxin-producing Escherichia coli O157:H7 caused by the consumption of raw drinking milk. There were seven discrepant variants called between the two technologies, five were false-negative or false-positive variants in the Illumina data and two were false-negative calls in ONT data. After masking horizontally acquired sequences such as prophages, analysis of both short and long-read sequences revealed the 20 isolates linked to the outbreak in 2017 had a maximum SNP distance of one SNP between each other, and a maximum of five SNPs when including three additional strains identified in 2019. Analysis of the ONT data revealed a 47 kbp deletion event in a terminal compound prophage within one sample relative to the remaining samples, and a 0.65 Mbp large chromosomal rearrangement (inversion), within one sample relative to the remaining samples. Furthermore, we detected two bacteriophages encoding the highly pathogenic Shiga toxin (Stx) subtype, Stx2a. One was typical of Stx2a-phage in this sub-lineage (Ic), the other was atypical and inserted into a site usually occupied by Stx2c-encoding phage. Finally, we observed an increase in the size of the pO157 IncFIB plasmid (1.6 kbp) in isolates from 2019 compared to those from 2017, due to the duplication of insertion elements within the plasmids from the more recently isolated strains. The ability to characterize the accessory genome in this way is the first step to understanding the significance of these microevolutionary events and their impact on the genome plasticity and virulence between strains of this zoonotic, foodborne pathogen.</p

Acquisition and loss of CTX-M plasmids in Shigella species associated with MSM transmission in the UK

Shigellosis in men who have sex with men (MSM) is caused by multidrug resistant Shigellae, exhibiting resistance to antimicrobials including azithromycin, ciprofloxacin and more recently the third-generation cephalosporins. We sequenced four bla (CTX-M-27)-positive MSM Shigella isolates (2018–20) using Oxford Nanopore Technologies; three S. sonnei (identified as two MSM clade 2, one MSM clade 5) and one S. flexneri 3a, to explore AMR context. All S. sonnei isolates harboured Tn7/Int2 chromosomal integrons, whereas S. flexneri 3a contained the Shigella Resistance Locus. All strains harboured IncFII pKSR100-like plasmids (67-83kbp); where present bla (CTX-M-27) was located on these plasmids flanked by IS26 and IS903B, however bla (CTX-M-27) was lost in S. flexneri 3a during storage between Illumina and Nanopore sequencing. IncFII AMR regions were mosaic and likely reorganised by IS26; three of the four plasmids contained azithromycin-resistance genes erm(B) and mph(A) and one harboured the pKSR100 integron. Additionally, all S. sonnei isolates possessed a large IncB/O/K/Z plasmid, two of which carried aph(3’)-Ib/aph(6)-Id/sul2 and tet(A). Monitoring the transmission of mobile genetic elements with co-located AMR determinants is necessary to inform empirical treatment guidance and clinical management of MSM-associated shigellosis

Double Averaging Analysis Applied to a Large Eddy Simulation of Coupled Turbulent Overlying and Porewater Flow

Freestream turbulence in rivers is a key contributor to the flux of dissolved nutrients, carbon, and other ecologically important solutes into porewater. To advance understanding of turbulent hyporheic exchange and porewater transport, we investigate flow over and through a rough bed of spheres using large eddy simulation (LES). We apply double averaging (combined space and time averaging) to the LES results to determine the mean velocity distribution, momentum balance, and drag forces. Our simulations show large-scale freestream structures interacting strongly with vortices generated at the surfaces of individual spheres to control turbulent momentum fluxes into the bed. The transition between turbulent flow and Darcy flow occurs over the first row of spheres, where turbulence decays rapidly and turbulent kinetic energy, Reynolds stress, and drag forces peak. Below this region, turbulence is only present in the high-velocity flow in open pore throats. Experimental observations suggest that minimum mean porewater velocity occurs in the first open pore space below the transition region, but our results show that the minimum occurs between the first and second pore spaces. The simulation mean porewater velocities are approximately half those captured in measurements because the model resolves the entire flow continuum while measurements can access high-velocity fluid in open pores. The high-resolution dual time-space averaging of the LES resolves both turbulent and mean flow features that are important to interfacial solute and particle fluxes, providing a means to include turbulent hyporheic exchange in upscaled river models, which has not been achieved to date