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.

Telomere-to-telomere assembly of a complete human X chromosome.

After two decades of improvements, the current human reference genome (GRCh38) is the most accurate and complete vertebrate genome ever produced. However, no single chromosome has been finished end to end, and hundreds of unresolved gaps persist1,2. Here we present a human genome assembly that surpasses the continuity of GRCh382, along with a gapless, telomere-to-telomere assembly of a human chromosome. This was enabled by high-coverage, ultra-long-read nanopore sequencing of the complete hydatidiform mole CHM13 genome, combined with complementary technologies for quality improvement and validation. Focusing our efforts on the human X chromosome3, we reconstructed the centromeric satellite DNA array (approximately 3.1 Mb) and closed the 29 remaining gaps in the current reference, including new sequences from the human pseudoautosomal regions and from cancer-testis ampliconic gene families (CT-X and GAGE). These sequences will be integrated into future human reference genome releases. In addition, the complete chromosome X, combined with the ultra-long nanopore data, allowed us to map methylation patterns across complex tandem repeats and satellite arrays. Our results demonstrate that finishing the entire human genome is now within reach, and the data presented here will facilitate ongoing efforts to complete the other human chromosomes.

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.

The Impact of Real-Time Whole-Genome Sequencing in Controlling Healthcare-Associated SARS-CoV-2 Outbreaks.

Nosocomial severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections have severely affected bed capacity and patient flow. We utilized whole-genome sequencing (WGS) to identify outbreaks and focus infection control resources and intervention during the United Kingdom's second pandemic wave in late 2020. Phylogenetic analysis of WGS and epidemiological data pinpointed an initial transmission event to an admission ward, with immediate prior community infection linkage documented. High incidence of asymptomatic staff infection with genetically identical viral sequences was also observed, which may have contributed to the propagation of the outbreak. WGS allowed timely nosocomial transmission intervention measures, including admissions ward point-of-care testing and introduction of portable HEPA14 filters. Conversely, WGS excluded nosocomial transmission in 2 instances with temporospatial linkage, conserving time and resources. In summary, WGS significantly enhanced understanding of SARS-CoV-2 clusters in a hospital setting, both identifying high-risk areas and conversely validating existing control measures in other units, maintaining clinical service overall.

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.

Retrospective screening of routine respiratory samples revealed undetected community transmission and missed intervention opportunities for SARS-CoV-2 in the United Kingdom.

In the early phases of the SARS coronavirus type 2 (SARS-CoV-2) pandemic, testing focused on individuals fitting a strict case definition involving a limited set of symptoms together with an identified epidemiological risk, such as contact with an infected individual or travel to a high-risk area. To assess whether this impaired our ability to detect and control early introductions of the virus into the UK, we PCR-tested archival specimens collected on admission to a large UK teaching hospital who retrospectively were identified as having a clinical presentation compatible with COVID-19. In addition, we screened available archival specimens submitted for respiratory virus diagnosis, and dating back to early January 2020, for the presence of SARS-CoV-2 RNA. Our data provides evidence for widespread community circulation of SARS-CoV-2 in early February 2020 and into March that was undetected at the time due to restrictive case definitions informing testing policy. Genome sequence data showed that many of these early cases were infected with a distinct lineage of the virus. Sequences obtained from the first officially recorded case in Nottinghamshire - a traveller returning from Daegu, South Korea - also clustered with these early UK sequences suggesting acquisition of the virus occurred in the UK and not Daegu. Analysis of a larger sample of sequences obtained in the Nottinghamshire area revealed multiple viral introductions, mainly in late February and through March. These data highlight the importance of timely and extensive community testing to prevent future widespread transmission of the virus.

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.

Whole genome sequencing of Salmonella Typhimurium illuminates distinct outbreaks caused by an endemic multi-locus variable number tandem repeat analysis type in Australia, 2014.

BACKGROUND: Salmonella Typhimurium (STM) is an important cause of foodborne outbreaks worldwide. Subtyping of STM remains critical to outbreak investigation, yet current techniques (e.g. multilocus variable number tandem repeat analysis, MLVA) may provide insufficient discrimination. Whole genome sequencing (WGS) offers potentially greater discriminatory power to support infectious disease surveillance. METHODS: We performed WGS on 62 STM isolates of a single, endemic MLVA type associated with two epidemiologically independent, food-borne outbreaks along with sporadic cases in New South Wales, Australia, during 2014. Genomes of case and environmental isolates were sequenced using HiSeq (Illumina) and the genetic distance between them was assessed by single nucleotide polymorphism (SNP) analysis. SNP analysis was compared to the epidemiological context. RESULTS: The WGS analysis supported epidemiological evidence and genomes of within-outbreak isolates were nearly identical. Sporadic cases differed from outbreak cases by a small number of SNPs, although their close relationship to outbreak cases may represent an unidentified common food source that may warrant further public health follow up. Previously unrecognised mini-clusters were detected. CONCLUSIONS: WGS of STM can discriminate foodborne community outbreaks within a single endemic MLVA clone. Our findings support the translation of WGS into public health laboratory surveillance of salmonellosis.

It Is Not All about Single Nucleotide Polymorphisms: Comparison of Mobile Genetic Elements and Deletions in Listeria monocytogenes Genomes Links Cases of Hospital-Acquired Listeriosis to the Environmental Source.

The control of food-borne outbreaks caused by Listeria monocytogenes in humans relies on the timely identification of food or environmental sources and the differentiation of outbreak-related isolates from unrelated ones. This study illustrates the utility of whole-genome sequencing for examining the link between clinical and environmental isolates of L. monocytogenes associated with an outbreak of hospital-acquired listeriosis in Sydney, Australia. Comparative genomic analysis confirmed an epidemiological link between the three clinical and two environmental isolates. Single nucleotide polymorphism (SNP) analysis showed that only two SNPs separated the three human outbreak isolates, which differed by 19 to 20 SNPs from the environmental isolates and 71 to >10,000 SNPs from sporadic L. monocytogenes isolates. The chromosomes of all human outbreak isolates and the two suspected environmental isolates were syntenic. In contrast to the genomes of background sporadic isolates, all epidemiologically linked isolates contained two novel prophages and a previously unreported clustered regularly interspaced short palindromic repeat (CRISPR)-associated (Cas) locus subtype sequence. The mobile genetic element (MGE) profile of these isolates was distinct from that of the other serotype 1/2b reference strains and sporadic isolates. The identification of SNPs and clonally distinctive MGEs strengthened evidence to distinguish outbreak-related isolates of L. monocytogenes from cocirculating endemic strains.

Genome sequencing reveals novel deletions associated with secondary resistance to pyrazinamide in MDR Mycobacterium tuberculosis.

OBJECTIVES: Detection of pyrazinamide resistance in Mycobacterium tuberculosis isolates presents significant challenges in settings with no dominant clonal lineages, such as Australia. We assessed the utility of WGS versus standard PCR amplification assays for the characterization of pyrazinamide resistance in MDR-TB isolates identified in New South Wales, Australia, over an 8 year period. METHODS: PCR amplicon sequencing was used to identify molecular markers associated with antibiotic resistance in pyrazinamide-resistant MDR-TB isolates recovered by the New South Wales Mycobacterium Reference Laboratory between 2007 and 2014. WGS was subsequently performed on two isolates for which pncA amplification failed. RESULTS: WGS identified two novel genomic deletions associated with in vitro resistance to pyrazinamide in MDR-TB. One isolate also carried a second deletion involving the genes dfrA and thyA associated with resistance to para-aminosalicylic acid. CONCLUSIONS: Steadily decreasing sequencing costs are increasing the appeal of WGS as an alternative approach for detecting complex patterns of pyrazinamide resistance in MDR-TB.

A single-cell atlas of pig gastrulation as a resource for comparative embryology.

Cell-fate decisions during mammalian gastrulation are poorly understood outside of rodent embryos. The embryonic disc of pig embryos mirrors humans, making them a useful proxy for studying gastrulation. Here we present a single-cell transcriptomic atlas of pig gastrulation, revealing cell-fate emergence dynamics, as well as conserved and divergent gene programs governing early porcine, primate, and murine development. We highlight heterochronicity in extraembryonic cell-types, despite the broad conservation of cell-type-specific transcriptional programs. We apply these findings in combination with functional investigations, to outline conserved spatial, molecular, and temporal events during definitive endoderm specification. We find early FOXA2 + /TBXT- embryonic disc cells directly form definitive endoderm, contrasting later-emerging FOXA2/TBXT+ node/notochord progenitors. Unlike mesoderm, none of these progenitors undergo epithelial-to-mesenchymal transition. Endoderm/Node fate hinges on balanced WNT and hypoblast-derived NODAL, which is extinguished upon endodermal differentiation. These findings emphasise the interplay between temporal and topological signalling in fate determination during gastrulation.

Sub-Optimal Paternal Diet at the Time of Mating Disrupts Maternal Adaptations to Pregnancy in the Late Gestation Mouse.

Pregnancy represents a stage during which maternal physiology and homeostatic regulation undergo dramatic change and adaptation. The fundamental purpose of these adaptations is to ensure the survival of her offspring through adequate nutrient provision and an environment that is tolerant to the semi-allogenic foetus. While poor maternal diet during pregnancy is associated with perturbed maternal adaptations during pregnancy, the influence of paternal diet on maternal well-being is less clearly defined. We fed C57BL/6 male mice either a control (CD), low protein diet (LPD), a high fat/sugar Western diet (WD) or the LPD or WD supplemented with methyl donors (MD-LPD and MD-WD, respectively) for a minimum of 8 weeks prior to mating with C57BL/6 females. Mated females were culled at day 17 of gestation for the analysis of maternal metabolic, gut, cardiac and bone health. Paternal diet had minimal influences on maternal serum and hepatic metabolite levels or gut microbiota diversity. However, analysis of the maternal hepatic transcriptome revealed distinct profiles of differential gene expression in response to the diet of the father. Paternal LPD and MD-LPD resulted in differential expression of genes associated with lipid metabolism, transcription, ubiquitin conjugation and immunity in dams, while paternal WD and MD-WD modified the expression of genes associated with ubiquitin conjugation and cardiac morphology. Finally, we observed changes in maternal femur length, volume of trabecular bone, trabecular connectivity, volume of the cortical medullar cavity and thickness of the cortical bone in response to the father's diets. Our current study demonstrates that poor paternal diet at the time of mating can influence the patterns of maternal metabolism and gestation-associated adaptations to her physiology.

Genome Reference Assembly for Bottlenecked Southern Australian Koalas.

Koala populations show marked differences in inbreeding levels and in the presence or absence of the endogenous Koala retrovirus (KoRV). These genetic differences among populations may lead to severe disease impacts threatening koala population viability. In addition, the recent colonization of the koala genome by KoRV provides a unique opportunity to study the process of retroviral adaptation to vertebrate genomes and the impact this has on speciation, genome structure, and function. The genome build described here is from an animal from the bottlenecked Southern population free of endogenous and exogenous KoRV. It provides a more contiguous genome build than the previous koala reference derived from an animal from a more outbred Northern population and is the first koala genome from a KoRV polymerase-free animal.

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.

Use of Bulk Segregant Analysis for Determining the Genetic Basis of Azole Resistance in the Opportunistic Pathogen Aspergillus fumigatus.

A sexual cycle was described in 2009 for the opportunistic fungal pathogen Aspergillus fumigatus, opening up for the first time the possibility of using techniques reliant on sexual crossing for genetic analysis. The present study was undertaken to evaluate whether the technique 'bulk segregant analysis' (BSA), which involves detection of differences between pools of progeny varying in a particular trait, could be applied in conjunction with next-generation sequencing to investigate the underlying basis of monogenic traits in A. fumigatus. Resistance to the azole antifungal itraconazole was chosen as a model, with a dedicated bioinformatic pipeline developed to allow identification of SNPs that differed between the resistant progeny pool and resistant parent compared to the sensitive progeny pool and parent. A clinical isolate exhibiting monogenic resistance to itraconazole of unknown basis was crossed to a sensitive parent and F1 progeny used in BSA. In addition, the use of backcrossing and increasing the number in progeny pools was evaluated as ways to enhance the efficiency of BSA. Use of F1 pools of 40 progeny led to the identification of 123 candidate genes with SNPs distributed over several contigs when aligned to an A1163 reference genome. Successive rounds of backcrossing enhanced the ability to identify specific genes and a genomic region, with BSA of progeny (using 40 per pool) from a third backcross identifying 46 genes with SNPs, and BSA of progeny from a sixth backcross identifying 20 genes with SNPs in a single 292 kb region of the genome. The use of an increased number of 80 progeny per pool also increased the resolution of BSA, with 29 genes demonstrating SNPs between the different sensitive and resistant groupings detected using progeny from just the second backcross with the majority of variants located on the same 292 kb region. Further bioinformatic analysis of the 292 kb region identified the presence of a cyp51A gene variant resulting in a methionine to lysine (M220K) change in the CYP51A protein, which was concluded to be the causal basis of the observed resistance to itraconazole. The future use of BSA in genetic analysis of A. fumigatus is discussed.

Paternal low protein diet perturbs inter-generational metabolic homeostasis in a tissue-specific manner in mice.

The underlying mechanisms driving paternally-programmed metabolic disease in offspring remain poorly defined. We fed male C57BL/6 mice either a control normal protein diet (NPD; 18% protein) or an isocaloric low protein diet (LPD; 9% protein) for a minimum of 8 weeks. Using artificial insemination, in combination with vasectomised male mating, we generated offspring using either NPD or LPD sperm but in the presence of NPD or LPD seminal plasma. Offspring from either LPD sperm or seminal fluid display elevated body weight and tissue dyslipidaemia from just 3 weeks of age. These changes become more pronounced in adulthood, occurring in conjunction with altered hepatic metabolic and inflammatory pathway gene expression. Second generation offspring also display differential tissue lipid abundance, with profiles similar to those of first generation adults. These findings demonstrate that offspring metabolic homeostasis is perturbed in response to a suboptimal paternal diet with the effects still evident within a second generation.