NANOG is required to establish the competence for germ-layer differentiation in the basal tetrapod axolotl.

Pluripotency defines the unlimited potential of individual cells of vertebrate embryos, from which all adult somatic cells and germ cells are derived. Understanding how the programming of pluripotency evolved has been obscured in part by a lack of data from lower vertebrates; in model systems such as frogs and zebrafish, the function of the pluripotency genes NANOG and POU5F1 have diverged. Here, we investigated how the axolotl ortholog of NANOG programs pluripotency during development. Axolotl NANOG is absolutely required for gastrulation and germ-layer commitment. We show that in axolotl primitive ectoderm (animal caps; ACs) NANOG and NODAL activity, as well as the epigenetic modifying enzyme DPY30, are required for the mass deposition of H3K4me3 in pluripotent chromatin. We also demonstrate that all 3 protein activities are required for ACs to establish the competency to differentiate toward mesoderm. Our results suggest the ancient function of NANOG may be establishing the competence for lineage differentiation in early cells. These observations provide insights into embryonic development in the tetrapod ancestor from which terrestrial vertebrates evolved.

Icarust, a real-time simulator for Oxford Nanopore adaptive sampling.

MOTIVATION: Oxford Nanopore Technologies (ONT) sequencers enable real-time generation of sequence data, which allows for concurrent analysis during a run. Adaptive sampling leverages this real-time capability in extremis, rejecting or accepting reads for sequencing based on assessment of the sequence from the start of each read. This functionality is provided by ONT's software, MinKNOW (Oxford Nanopore Technologies). Designing and developing software to take advantage of adaptive sampling can be costly in terms of sequencing consumables, using precious samples and preparing sequencing libraries. MinKNOW addresses this in part by allowing the replay of previously sequenced runs for testing. However, as we show, the sequencing output only partially changes in response to adaptive sampling instructions. Here we present Icarust, a tool enabling more accurate approximations of sequencing runs. Icarust recreates all the required endpoints of MinKNOW to perform adaptive sampling and writes output compatible with current base-callers and analysis pipelines. Icarust serves nanopore signal simulating a MinION or PromethION flow cell experiment from any reference genome using either R9 or R10 pore models. We show that simulating sequencing runs with Icarust provides a realistic testing and development environment for software exploiting the real-time nature of Nanopore sequencing. AVAILABILITY AND IMPLEMENTATION: All code is open source and freely available here-https://github.com/LooseLab/Icarust. Icarust is implemented in Rust, with a docker container also available. The data underlying this article will be shared on reasonable request to the corresponding author.

Author Correction: Nanopore native RNA sequencing of a human poly(A) transcriptome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

minoTour, real-time monitoring and analysis for nanopore sequencers.

SUMMARY: minoTour offers a Laboratory Informations Management System (LIMS) system for Oxford Nanopore Technology sequencers, with real-time metrics and analysis available permanently for review. Integration of unique real-time automated analysis can reduce the time required to answer biological questions, including mapping and classification of sequence while a run is in progress. Real-time sequence data require new methods of analysis which do not wait for the completion of a run and minoTour provides a framework to allow users to exploit these features. AVAILABILITY AND IMPLEMENTATION: Source code and documentation are available at https://github.com/LooseLab/minotourcli and https://github.com/LooseLab/minotourapp. Docker images are available from https://hub.docker.com/r/adoni5/, and can be installed using a preconfigured docker-compose script at https://github.com/LooseLab/minotour-docker. An example server is available at http://137.44.59.170. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Nanopore native RNA sequencing of a human poly(A) transcriptome.

High-throughput complementary DNA sequencing technologies have advanced our understanding of transcriptome complexity and regulation. However, these methods lose information contained in biological RNA because the copied reads are often short and modifications are not retained. We address these limitations using a native poly(A) RNA sequencing strategy developed by Oxford Nanopore Technologies. Our study generated 9.9 million aligned sequence reads for the human cell line GM12878, using thirty MinION flow cells at six institutions. These native RNA reads had a median length of 771 bases, and a maximum aligned length of over 21,000 bases. Mitochondrial poly(A) reads provided an internal measure of read-length quality. We combined these long nanopore reads with higher accuracy short-reads and annotated GM12878 promoter regions to identify 33,984 plausible RNA isoforms. We describe strategies for assessing 3' poly(A) tail length, base modifications and transcript haplotypes.

tartan underlies the evolution of Drosophila male genital morphology.

Male genital structures are among the most rapidly evolving morphological traits and are often the only features that can distinguish closely related species. This process is thought to be driven by sexual selection and may reinforce species separation. However, while the genetic bases of many phenotypic differences have been identified, we still lack knowledge about the genes underlying evolutionary differences in male genital organs and organ size more generally. The claspers (surstyli) are periphallic structures that play an important role in copulation in insects. Here, we show that divergence in clasper size and bristle number between Drosophila mauritiana and Drosophila simulans is caused by evolutionary changes in tartan (trn), which encodes a transmembrane leucine-rich repeat domain protein that mediates cell-cell interactions and affinity. There are no fixed amino acid differences in trn between D. mauritiana and D. simulans, but differences in the expression of this gene in developing genitalia suggest that cis-regulatory changes in trn underlie the evolution of clasper morphology in these species. Finally, analyses of reciprocal hemizygotes that are genetically identical, except for the species from which the functional allele of trn originates, determined that the trn allele of D. mauritiana specifies larger claspers with more bristles than the allele of D. simulans Therefore, we have identified a gene underlying evolutionary change in the size of a male genital organ, which will help to better understand not only the rapid diversification of these structures, but also the regulation and evolution of organ size more broadly.

BulkVis: a graphical viewer for Oxford nanopore bulk FAST5 files.

MOTIVATION: The Oxford Nanopore Technologies (ONT) MinION is used for sequencing a wide variety of sample types with diverse methods of sample extraction. Nanopore sequencers output FAST5 files containing signal data subsequently base called to FASTQ format. Optionally, ONT devices can collect data from all sequencing channels simultaneously in a bulk FAST5 file enabling inspection of signal in any channel at any point. We sought to visualize this signal to inspect challenging or difficult to sequence samples. RESULTS: The BulkVis tool can load a bulk FAST5 file and overlays MinKNOW (the software that controls ONT sequencers) classifications on the signal trace and can show mappings to a reference. Users can navigate to a channel and time or, given a FASTQ header from a read, jump to its specific position. BulkVis can export regions as Nanopore base caller compatible reads. Using BulkVis, we find long reads can be incorrectly divided by MinKNOW resulting in single DNA molecules being split into two or more reads. The longest seen to date is 2 272 580 bases in length and reported in eleven consecutive reads. We provide helper scripts that identify and reconstruct split reads given a sequencing summary file and alignment to a reference. We note that incorrect read splitting appears to vary according to input sample type and is more common in 'ultra-long' read preparations. AVAILABILITY AND IMPLEMENTATION: The software is available freely under an MIT license at https://github.com/LooseLab/bulkvis. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Physiological temperatures reduce dimerization of dengue and Zika virus recombinant envelope proteins.

The spread of dengue (DENV) and Zika virus (ZIKV) is a major public health concern. The primary target of antibodies that neutralize DENV and ZIKV is the envelope (E) glycoprotein, and there is interest in using soluble recombinant E (sRecE) proteins as subunit vaccines. However, the most potent neutralizing antibodies against DENV and ZIKV recognize epitopes on the virion surface that span two or more E proteins. Therefore, to create effective DENV and ZIKV vaccines, presentation of these quaternary epitopes may be necessary. The sRecE proteins from DENV and ZIKV crystallize as native-like dimers, but studies in solution suggest that these dimers are marginally stable. To better understand the challenges associated with creating stable sRecE dimers, we characterized the thermostability of sRecE proteins from ZIKV and three DENV serotypes, DENV2-4. All four proteins irreversibly unfolded at moderate temperatures (46-53 °C). At 23 °C and low micromolar concentrations, DENV2 and ZIKV were primarily dimeric, and DENV3-4 were primarily monomeric, whereas at 37 °C, all four proteins were predominantly monomeric. We further show that the dissociation constant for DENV2 dimerization is very temperature-sensitive, ranging from <1 µm at 25 °C to 50 µm at 41 °C, due to a large exothermic enthalpy of binding of -79 kcal/mol. We also found that quaternary epitope antibody binding to DENV2-4 and ZIKV sRecE is reduced at 37 °C. Our observation of reduced sRecE dimerization at physiological temperature highlights the need for stabilizing the dimer as part of its development as a subunit vaccine.

Eye movements show similar adaptations in temporal coordination to movement planning conditions in both people with and without cerebral palsy.

The effects of cerebral palsy on movement planning for simple reaching tasks are not well understood. Movement planning is complex and entails many processes which could be affected. This study specifically sought to evaluate integrating task information, decoupling movements, and adjusting to altered mapping. For a reaching task, the asynchrony between the eye onset and the hand onset was measured across different movement planning conditions for participants with and without cerebral palsy. Previous research shows people without cerebral palsy vary this temporal coordination for different planning conditions. Our measurements show similar adaptations in temporal coordination for groups with and without cerebral palsy, to three of the four variations in planning condition tested. However, movement durations were still longer for the participants with cerebral palsy. Hence for simple goal-directed reaching, movement execution problems appear to limit activity more than movement planning deficits.

Decoupling speciation and extinction reveals both abiotic and biotic drivers shaped 250 million years of diversity in crocodile-line archosaurs.

Whereas living representatives of Pseudosuchia, crocodylians, number fewer than 30 species, more than 700 pseudosuchian species are known from their 250-million-year fossil record, displaying far greater ecomorphological diversity than their extant counterparts. With a new time-calibrated tree of >500 species, we use a phylogenetic framework to reveal that pseudosuchian evolutionary history and diversification dynamics were directly shaped by the interplay of abiotic and biotic processes over hundreds of millions of years, supported by information theory analyses. Speciation, but not extinction, is correlated with higher temperatures in terrestrial and marine lineages, with high sea level associated with heightened extinction in non-marine taxa. Low lineage diversity and increased speciation in non-marine species is consistent with opportunities for niche-filling, whereas increased competition may have led to elevated extinction rates. In marine lineages, competition via increased lineage diversity appears to have driven both speciation and extinction. Decoupling speciation and extinction, in combination with ecological partitioning, reveals a more complex picture of pseudosuchian evolution than previously understood. As the number of species threatened with extinction by anthropogenic climate change continues to rise, the fossil record provides a unique window into the drivers that led to clade success and those that may ultimately lead to extinction.

Socioeconomic deprivation and illness trajectory in the Scottish population after COVID-19 hospitalization.

BACKGROUND: The associations between deprivation and illness trajectory after hospitalisation for coronavirus disease-19 (COVID-19) are uncertain. METHODS: A prospective, multicentre cohort study was conducted on post-COVID-19 patients, enrolled either in-hospital or shortly post-discharge. Two evaluations were carried out: an initial assessment and a follow-up at 28-60 days post-discharge. The study encompassed research blood tests, patient-reported outcome measures, and multisystem imaging (including chest computed tomography (CT) with pulmonary and coronary angiography, cardiovascular and renal magnetic resonance imaging). Primary and secondary outcomes were analysed in relation to socioeconomic status, using the Scottish Index of Multiple Deprivation (SIMD). The EQ-5D-5L, Brief Illness Perception Questionnaire (BIPQ), Patient Health Questionnaire-4 (PHQ-4) for Anxiety and Depression, and the Duke Activity Status Index (DASI) were used to assess health status. RESULTS: Of the 252 enrolled patients (mean age 55.0 ± 12.0 years; 40% female; 23% with diabetes), deprivation status was linked with increased BMI and diabetes prevalence. 186 (74%) returned for the follow-up. Within this group, findings indicated associations between deprivation and lung abnormalities (p = 0.0085), coronary artery disease (p = 0.0128), and renal inflammation (p = 0.0421). Furthermore, patients with higher deprivation exhibited worse scores in health-related quality of life (EQ-5D-5L, p = 0.0084), illness perception (BIPQ, p = 0.0004), anxiety and depression levels (PHQ-4, p = 0.0038), and diminished physical activity (DASI, p = 0.002). At the 3-month mark, those with greater deprivation showed a higher frequency of referrals to secondary care due to ongoing COVID-19 symptoms (p = 0.0438). However, clinical outcomes were not influenced by deprivation. CONCLUSIONS: In a post-hospital COVID-19 population, socioeconomic deprivation was associated with impaired health status and secondary care episodes. Deprivation influences illness trajectory after COVID-19.

In our study, we aimed to understand how socioeconomic factors impact recovery from COVID-19 following hospitalisation. We followed 252 patients, collecting health data and utilising advanced imaging techniques. We discovered that individuals from deprived areas experienced more severe health complications, reported worse quality of life, and required more specialist care. However, their clinical outcomes were not significantly different. This underscores that socioeconomic deprivation affects health recovery, underlining the need for tailored care for these individuals. Our findings emphasise the importance of considering socioeconomic factors in recovery plans post-COVID-19, potentially improving healthcare for those in deprived areas.

Highly automatic quantification of myocardial oedema in patients with acute myocardial infarction using bright blood T2-weighted CMR.

BACKGROUND: T2-weighted cardiovascular magnetic resonance (CMR) is clinically-useful for imaging the ischemic area-at-risk and amount of salvageable myocardium in patients with acute myocardial infarction (MI). However, to date, quantification of oedema is user-defined and potentially subjective. METHODS: We describe a highly automatic framework for quantifying myocardial oedema from bright blood T2-weighted CMR in patients with acute MI. Our approach retains user input (i.e. clinical judgment) to confirm the presence of oedema on an image which is then subjected to an automatic analysis. The new method was tested on 25 consecutive acute MI patients who had a CMR within 48 hours of hospital admission. Left ventricular wall boundaries were delineated automatically by variational level set methods followed by automatic detection of myocardial oedema by fitting a Rayleigh-Gaussian mixture statistical model. These data were compared with results from manual segmentation of the left ventricular wall and oedema, the current standard approach. RESULTS: The mean perpendicular distances between automatically detected left ventricular boundaries and corresponding manual delineated boundaries were in the range of 1-2 mm. Dice similarity coefficients for agreement (0=no agreement, 1=perfect agreement) between manual delineation and automatic segmentation of the left ventricular wall boundaries and oedema regions were 0.86 and 0.74, respectively. CONCLUSION: Compared to standard manual approaches, the new highly automatic method for estimating myocardial oedema is accurate and straightforward. It has potential as a generic software tool for physicians to use in clinical practice.

A framework for real-time monitoring, analysis and adaptive sampling of viral amplicon nanopore sequencing.

The ongoing SARS-CoV-2 pandemic demonstrates the utility of real-time sequence analysis in monitoring and surveillance of pathogens. However, cost-effective sequencing requires that samples be PCR amplified and multiplexed via barcoding onto a single flow cell, resulting in challenges with maximising and balancing coverage for each sample. To address this, we developed a real-time analysis pipeline to maximise flow cell performance and optimise sequencing time and costs for any amplicon based sequencing. We extended our nanopore analysis platform MinoTour to incorporate ARTIC network bioinformatics analysis pipelines. MinoTour predicts which samples will reach sufficient coverage for downstream analysis and runs the ARTIC networks Medaka pipeline once sufficient coverage has been reached. We show that stopping a viral sequencing run earlier, at the point that sufficient data has become available, has no negative effect on subsequent down-stream analysis. A separate tool, SwordFish, is used to automate adaptive sampling on Nanopore sequencers during the sequencing run. This enables normalisation of coverage both within (amplicons) and between samples (barcodes) on barcoded sequencing runs. We show that this process enriches under-represented samples and amplicons in a library as well as reducing the time taken to obtain complete genomes without affecting the consensus sequence.

Evaluation of variant calling algorithms for wastewater-based epidemiology using mixed populations of SARS-CoV-2 variants in synthetic and wastewater samples.

Wastewater-based epidemiology has been used extensively throughout the COVID-19 (coronavirus disease 19) pandemic to detect and monitor the spread and prevalence of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) and its variants. It has proven an excellent, complementary tool to clinical sequencing, supporting the insights gained and helping to make informed public-health decisions. Consequently, many groups globally have developed bioinformatics pipelines to analyse sequencing data from wastewater. Accurate calling of mutations is critical in this process and in the assignment of circulating variants; yet, to date, the performance of variant-calling algorithms in wastewater samples has not been investigated. To address this, we compared the performance of six variant callers (VarScan, iVar, GATK, FreeBayes, LoFreq and BCFtools), used widely in bioinformatics pipelines, on 19 synthetic samples with known ratios of three different SARS-CoV-2 variants of concern (VOCs) (Alpha, Beta and Delta), as well as 13 wastewater samples collected in London between the 15th and 18th December 2021. We used the fundamental parameters of recall (sensitivity) and precision (specificity) to confirm the presence of mutational profiles defining specific variants across the six variant callers. Our results show that BCFtools, FreeBayes and VarScan found the expected variants with higher precision and recall than GATK or iVar, although the latter identified more expected defining mutations than other callers. LoFreq gave the least reliable results due to the high number of false-positive mutations detected, resulting in lower precision. Similar results were obtained for both the synthetic and wastewater samples.

Utility of wastewater genomic surveillance compared to clinical surveillance to track the spread of the SARS-CoV-2 Omicron variant across England.

The world has moved into a new stage of managing the SARS-CoV-2 pandemic with minimal restrictions and reduced testing in the population, leading to reduced genomic surveillance of virus variants in individuals. Wastewater-based epidemiology (WBE) can provide an alternative means of tracking virus variants in the population but decision-makers require confidence that it can be applied to a national scale and is comparable to individual testing data. We analysed 19,911 samples from 524 wastewater sites across England at least twice a week between November 2021 and February 2022, capturing sewage from >70% of the English population. We used amplicon-based sequencing and the phylogeny based de-mixing tool Freyja to estimate SARS-CoV-2 variant frequencies and compared these to the variant dynamics observed in individual testing data from clinical and community settings. We show that wastewater data can reconstruct the spread of the Omicron variant across England since November 2021 in close detail and aligns closely with epidemiological estimates from individual testing data. We also show the temporal and spatial spread of Omicron within London. Our wastewater data further reliably track the transition between Omicron subvariants BA1 and BA2 in February 2022 at regional and national levels. Our demonstration that WBE can track the fast-paced dynamics of SARS-CoV-2 variant frequencies at a national scale and closely match individual testing data in time shows that WBE can reliably fill the monitoring gap left by reduced individual testing in a more affordable way.

Adjudicated myocarditis and multisystem illness trajectory in healthcare workers post-COVID-19.

BACKGROUND: We investigated the associations of healthcare worker status with multisystem illness trajectory in hospitalised post-COVID-19 individuals. METHODS AND RESULTS: One hundred and sixty-eight patients were evaluated 28-60 days after the last episode of hospital care. Thirty-six (21%) were healthcare workers. Compared with non-healthcare workers, healthcare workers were of similar age (51.3 (8.7) years vs 55.0 (12.4) years; p=0.09) more often women (26 (72%) vs 48 (38%); p<0.01) and had lower 10-year cardiovascular risk (%) (8.1 (7.9) vs 15.0 (11.5); p<0.01) and Coronavirus Clinical Characterisation Consortium in-hospital mortality risk (7.3 (10.2) vs 12.7 (9.8); p<0.01). Healthcare worker status associated with less acute inflammation (peak C reactive protein 48 mg/L (IQR: 14-165) vs 112 mg/L (52-181)), milder illness reflected by WHO clinical severity score distribution (p=0.04) and shorter duration of admission (4 days (IQR: 2-6) vs 6 days (3-12)).In adjusted multivariate logistic regression analysis, healthcare worker status associated with a binary classification (probable/very likely vs not present/unlikely) of adjudicated myocarditis (OR: 2.99; 95% CI (1.01 to 8.89) by 28-60 days postdischarge).After a mean (SD, range) duration of follow-up after hospital discharge of 450 (88) days (range 290, 627 days), fewer healthcare workers died or were rehospitalised (1 (3%) vs 22 (17%); p=0.038) and secondary care referrals for post-COVID-19 syndrome were common (42%) and similar to non-healthcare workers (38%; p=0.934). CONCLUSION: Healthcare worker status was independently associated with the likelihood of adjudicated myocarditis, despite better antecedent health. Two in five healthcare workers had a secondary care referral for post-COVID-19 syndrome. TRIAL REGISTRATION NUMBER: NCT04403607.

Ecological Transitions and the Shape of the Decapod Tree of Life.

Understanding the processes that shaped the distribution of species richness across the Tree of Life is a central macroevolutionary research agenda. Major ecological innovations, including transitions between habitats, may help to explain the striking asymmetries of diversity that are often observed between sister clades. Here, we test the impact of such transitions on speciation rates across decapod crustaceans, modeling diversification dynamics within a phylogenetic framework. Our results show that, while terrestrial lineages have higher speciation rates than either marine or freshwater lineages, there is no difference between mean speciation rates in marine and freshwater lineages across Decapoda. Partitioning our data by infraorder reveals that those clades with habitat heterogeneity have higher speciation rates in freshwater and terrestrial lineages, with freshwater rates up to 1.5 times faster than marine rates, and terrestrial rates approximately four times faster. This averaging out of marine and freshwater speciation rates results from the varying contributions of different clades to average speciation rates. However, with the exception of Caridea, we find no evidence for any causal relationship between habitat and speciation rate. Our results demonstrate that while statistical generalizations about ecological traits and evolutionary rates are valuable, there are many exceptions. Hence, while freshwater and terrestrial lineages typically speciate faster than their marine relatives, there are many atypically slow freshwater lineages and fast marine lineages across Decapoda. Future work on diversification patterns will benefit from the inclusion of fossil data, as well as additional ecological factors.

A multisystem, cardio-renal investigation of post-COVID-19 illness.

The pathophysiology and trajectory of post-Coronavirus Disease 2019 (COVID-19) syndrome is uncertain. To clarify multisystem involvement, we undertook a prospective cohort study including patients who had been hospitalized with COVID-19 (ClinicalTrials.gov ID NCT04403607 ). Serial blood biomarkers, digital electrocardiography and patient-reported outcome measures were obtained in-hospital and at 28-60 days post-discharge when multisystem imaging using chest computed tomography with pulmonary and coronary angiography and cardio-renal magnetic resonance imaging was also obtained. Longer-term clinical outcomes were assessed using electronic health records. Compared to controls (n = 29), at 28-60 days post-discharge, people with COVID-19 (n = 159; mean age, 55 years; 43% female) had persisting evidence of cardio-renal involvement and hemostasis pathway activation. The adjudicated likelihood of myocarditis was 'very likely' in 21 (13%) patients, 'probable' in 65 (41%) patients, 'unlikely' in 56 (35%) patients and 'not present' in 17 (11%) patients. At 28-60 days post-discharge, COVID-19 was associated with worse health-related quality of life (EQ-5D-5L score 0.77 (0.23) versus 0.87 (0.20)), anxiety and depression (PHQ-4 total score 3.59 (3.71) versus 1.28 (2.67)) and aerobic exercise capacity reflected by predicted maximal oxygen utilization (20.0 (7.6) versus 29.5 (8.0) ml/kg/min) (all P < 0.01). During follow-up (mean, 450 days), 24 (15%) patients and two (7%) controls died or were rehospitalized, and 108 (68%) patients and seven (26%) controls received outpatient secondary care (P = 0.017). The illness trajectory of patients after hospitalization with COVID-19 includes persisting multisystem abnormalities and health impairments that could lead to substantial demand on healthcare services in the future.