A laboratory framework for ongoing optimization of amplification-based genomic surveillance programs.

IMPORTANCE: This study provides a laboratory framework to ensure ongoing relevance and performance of amplification-based whole genome sequencing to strengthen public health surveillance during extended outbreaks or pandemics. The framework integrates regular reviews of the performance of a genomic surveillance system and highlights the importance of ongoing monitoring and the identification and implementation of improvements to whole genome sequencing methods to enhance public health responses to pathogen outbreaks.

Exploring programmatic indicators of tuberculosis control that incorporate routine Mycobacterium tuberculosis sequencing in low incidence settings: a comprehensive (2017-2021) patient cohort analysis.

Background: Routine whole genome sequencing of Mycobacterium tuberculosis has been implemented with increasing frequency. However, its value for tuberculosis (TB) control programs beyond individual case management and enhanced drug resistance detection has not yet been explored. Methods: We analysed routine sequencing data of culture-confirmed TB cases notified between 1st January 2017 and 31st December 2021 in New South Wales (NSW), Australia. Genomic surveillance included evidence of local TB transmission, defined by single nucleotide polymorphism (SNP) clustering over a variable (0-25) SNP threshold, and drug resistance conferring mutations. Findings: M. tuberculosis sequences from 1831 patients were examined, representing 64.8% of all notified TB cases and 96.2% of culture-confirmed cases. Applying a traditional 5-SNP cluster threshold identified 62 transmission clusters with 183 clustered cases; 101/183 (55.2%) had 0 SNP differences. Cluster assessment over a 5-year period, using a 5-SNP threshold, provided a comprehensive overview of likely recent transmission within NSW, Australia, as an indicator of local TB control. Genotypic drug susceptibility testing (DST) was highly concordant with phenotypic DST and provided a 6.8% increase in antimycobacterial resistance detection. Importantly, it detected mutations missed by routine molecular tests. Lineage 2 strains were more likely to be drug resistant (p < 0.0001) and locally transmitted if drug resistant (p < 0.0001). Interpretation: Performing routine prospective WGS in a low incidence country like Australia, provides genomically informed programmatic indicators of local TB control. A rolling 5-year cluster assessment reflects epidemic containment and progress towards 'zero TB transmission'. Genomic DST also provides valuable information for clinical care and drug resistance surveillance. Funding: NHMRC Centre for Research Excellence in Tuberculosis (www.tbcre.org.au) and NSW Health Prevention Research Support Program.

SABRes: in silico detection of drug resistance conferring mutations in subpopulations of SARS-CoV-2 genomes.

The emergence of resistance to antiviral drugs increasingly used to treat SARS-CoV-2 infections has been recognised as a significant threat to COVID-19 control. In addition, some SARS-CoV-2 variants of concern appear to be intrinsically resistant to several classes of these antiviral agents. Therefore, there is a critical need for rapid recognition of clinically relevant polymorphisms in SARS-CoV-2 genomes associated with significant reduction of drug activity in virus neutralisation experiments. Here we present SABRes, a bioinformatic tool, which leverages on expanding public datasets of SARS-CoV-2 genomes and allows detection of drug resistance mutations in consensus genomes as well as in viral subpopulations. We have applied SABRes to detect resistance-conferring mutations in 25,197 genomes generated over the course of the SARS-CoV-2 pandemic in Australia and identified 299 genomes containing resistance conferring mutations to the five antiviral therapeutics that retain effectiveness against currently circulating strains of SARS-CoV-2 - Sotrovimab, Bebtelovimab, Remdesivir, Nirmatrelvir and Molnupiravir. These genomes accounted for a 1.18% prevalence of resistant isolates discovered by SABRes, including 80 genomes with resistance conferring mutations found in viral subpopulations. Timely recognition of these mutations within subpopulations is critical as these mutations can provide an advantage under selective pressure and presents an important step forward in our ability to monitor SARS-CoV-2 drug resistance.

Genome entropy and network centrality contrast exploration and exploitation in evolution of foodborne pathogens.

Modelling evolution of foodborne pathogens is crucial for mitigation and prevention of outbreaks. We apply network-theoretic and information-theoretic methods to trace evolutionary pathways ofSalmonellaTyphimurium in New South Wales, Australia, by studying whole genome sequencing surveillance data over a five-year period which included several outbreaks. The study derives both undirected and directed genotype networks based on genetic proximity, and relates the network's structural property (centrality) to its functional property (prevalence). The centrality-prevalence space derived for the undirected network reveals a salient exploration-exploitation distinction across the pathogens, further quantified by the normalised Shannon entropy and the Fisher information of the corresponding shell genome. This distinction is also analysed by tracing the probability density along evolutionary paths in the centrality-prevalence space. We quantify the evolutionary pathways, and show that pathogens exploring the evolutionary search-space during the considered period begin to exploit their environment (their prevalence increases resulting in outbreaks), but eventually encounter a bottleneck formed by epidemic containment measures.

Emergent Omicron BR.2.1 sublineage of SARS-CoV-2 in New South Wales, Australia: a subvariant with high fitness but without increased disease severity.

OBJECTIVES: To describe the epidemiology and impact of Omicron BR.2.1, an emergent SARS-CoV-2 Omicron BA.2.75 sublineage displaying high fitness compared to other cocirculating subvariants in New South Wales, Australia. METHODS: From September 01 to November 26, 2022, 4971 SARS-CoV-2 consensus genomes from unique patients were generated, and correlated with international travel and reinfection history, and admission to the intensive care unit. RESULTS: BR.2.1 became the predominant variant by late November, and was responsible for a significantly higher proportion of community-acquired cases during the study period (55.1% vs 38.4%, P < 0.001). Reinfections (defined as occurring between 6 and 24 weeks after a prior diagnosis of COVID-19) were significantly higher among BR.2.1 compared to non-BR.2.1 infected persons (17.0% vs 6.0%, P < 0.001). BR.2.1 cases were also significantly younger compared to non-BR.2.1 (median age 48 years (interquartile range [IQR] 32) vs 53 years (IQR 32), P = 0.004). The proportion of patients admitted to the intensive care unit with BR.2.1 was not significantly higher than other subvariants (2.3% vs 2.0%, P = 0.717). CONCLUSION: Having emerged locally within New South Wales, BR.2.1 caused a significant number of SARS-CoV-2 reinfections, but with disease severity comparable with other currently circulating lineages. Given its rapid rise in prevalence, BR.2.1 has the potential to become established internationally.

Genomic Epidemiology and Multilevel Genome Typing of Australian Salmonella enterica Serovar Enteritidis.

Salmonella enterica serovar Enteritidis is one of the leading causes of salmonellosis in Australia. In this study, a total of 568 S. Enteritidis isolates from two Australian states across two consecutive years were analyzed and compared to international strains, using the S. Enteritidis multilevel genome typing (MGT) database, which contained 40,390 publicly available genomes from 99 countries. The Australian S. Enteritidis isolates were divided into three phylogenetic clades (A, B, and C). Clades A and C represented 16.4% and 3.5% of the total isolates, respectively, and were of local origin. Clade B accounted for 80.1% of the isolates which belonged to seven previously defined lineages but was dominated by the global epidemic lineage. At the MGT5 level, three out of five top sequence types (STs) in Australia were also top STs in Asia, suggesting that a fair proportion of Australian S. Enteritidis cases may be epidemiologically linked with Asian strains. In 2018, a large egg-associated local outbreak was caused by a recently defined clade B lineage prevalent in Europe and was closely related, but not directly linked, to three European isolates. Additionally, over half (54.8%) of predicted multidrug resistance (MDR) isolates belonged to 10 MDR-associated MGT-STs, which were also frequent in Asian S. Enteritidis . Overall, this study investigated the genomic epidemiology of S. Enteritidis in Australia, including the first large local outbreak, using MGT. The open MGT platform enables a standardized and sharable nomenclature that can be effectively applied to public health for unified surveillance of S. Enteritidis nationally and globally. IMPORTANCE Salmonella enterica serovar Enteritidis is a leading cause of foodborne infections. We previously developed a genomic typing database (MGTdb) for S. Enteritidis to facilitate global surveillance of this pathogen. In this study, we examined the genomic features of Australian S. Enteritidis using the MGTdb and found that Australian S. Enteritidis is mainly epidemiologically linked with Asian strains (especially strains carrying antimicrobial resistance genes), followed by European strains. The first large-scale egg-associated local outbreak in Australia was caused by a recently defined lineage prevalent in Europe, and three European isolates in the MGTdb were closely related but not directly linked to this outbreak. In summary, the S. Enteritidis MGTdb open platform is shown to be a potentially powerful tool for national and global public health surveillance of this pathogen.

Detection and characterisation of Bordetella hinzii in line-related bacteraemia and respiratory tract infection in Australia.

Bordetella hinzii has emerged as an unusual cause of infection in immunocompromised patients, previously linked to zoonotic transmission. Antimicrobial susceptibility and genetic diversity of B. hinzii are poorly understood. This study reports phenotypic and genomic characteristics of the first four Australian isolates of B. hinzii obtained from elderly immunocompromised patients. Bordetella hinzii isolates were identified by MALDI-TOF and whole genome sequencing (WGS). Antibiotic susceptibility testing was performed using disk diffusion or E-test. Genomes of B. hinzii were analysed in global context. A phylogenetic tree was constructed of all isolates using Roary and a maximum-likelihood tree was generated from the core-snp alignment. Bordetella hinzii minimum inhibitory concentrations (MICs) were largely uniform with high MICs to ampicillin, ceftriaxone and ciprofloxacin and low MICs to meropenem and piperacillin-tazobactam. Genomic analysis of isolate sequences divided strains analysed into two phylogenetically distinct groups, with one Australian B. hinzii isolate (AUS-4) assigned to Group 1, and the remaining isolates (AUS1-AUS3 and AUS-5) to Group 2. Single nucleotide polymorphism (SNP) analysis revealed two isolates, AUS-1 and AUS-2, were closely related with 14 SNP differences between them. All other Australian isolates were unrelated to each and all other isolates from the international dataset. Bordetella hinzii appears to pose a risk to immunocompromised individuals but remains susceptible to extended spectrum ß-lactam and carbapenem antibiotics. Genomic analysis suggested a dissemination of genetically distinct strains.

Guiding the design of SARS-CoV-2 genomic surveillance by estimating the resolution of outbreak detection.

Genomic surveillance of SARS-CoV-2 has been essential to inform public health response to outbreaks. The high incidence of infection has resulted in a smaller proportion of cases undergoing whole genome sequencing due to finite resources. We present a framework for estimating the impact of reduced depths of genomic surveillance on the resolution of outbreaks, based on a clustering approach using pairwise genetic and temporal distances. We apply the framework to simulated outbreak data to show that outbreaks are detected less frequently when fewer cases are subjected to whole genome sequencing. The impact of sequencing fewer cases depends on the size of the outbreaks, and on the genetic and temporal similarity of the index cases of the outbreaks. We also apply the framework to an outbreak of the SARS-CoV-2 Delta variant in New South Wales, Australia. We find that the detection of clusters in the outbreak would have been delayed if fewer cases had been sequenced. Existing recommendations for genomic surveillance estimate the minimum number of cases to sequence in order to detect and monitor new virus variants, assuming representative sampling of cases. Our method instead measures the resolution of clustering, which is important for genomic epidemiology, and accommodates sampling biases.

Whole Genome Sequencing Shows Genetic Diversity, as Well as Clonal Complex and Gene Polymorphisms Associated with Fluconazole Non-Susceptible Isolates of Candida tropicalis.

Resistance to azoles in Candida tropicalis is increasing and may be mediated by genetic characteristics. Using whole genome sequencing (WGS), we examined the genetic diversity of 82 bloodstream C. tropicalis isolates from two countries and one ATCC strain in a global context. Multilocus sequence typing (MLST) and single nucleotide polymorphism (SNP)-based phylogenies were generated. Minimum inhibitory concentrations (MIC) for antifungal agents were determined using Sensititre YeastOne YO10. Eleven (13.2%) isolates were fluconazole-resistant and 17 (20.5%) were classified as fluconazole-non susceptible (FNS). Together with four Canadian isolates, the genomes of 12 fluconazole-resistant (18 FNS) and 69 fluconazole-susceptible strains were examined for gene mutations associated with drug resistance. Fluconazole-resistant isolates contained a mean of 56 non-synonymous SNPs per isolate in contrast to 36 SNPs in fluconazole-susceptible isolates (interquartile range [IQR] 46−59 vs. 31−48 respectively; p < 0.001). Ten of 18 FNS isolates contained missense ERG11 mutations (amino acid substitutions S154F, Y132F, Y257H). Two echinocandin-non susceptible isolates had homozygous FKS1 mutations (S30P). MLST identified high genetic diversity with 61 diploid sequence types (DSTs), including 53 new DSTs. All four isolates in DST 773 were fluconazole-resistant within clonal complex 2. WGS showed high genetic variation in invasive C. tropicalis; azole resistance was distributed across different lineages but with DST 773 associated with in vitro fluconazole resistance.

SARS-CoV-2 Within-Host and in vitro Genomic Variability and Sub-Genomic RNA Levels Indicate Differences in Viral Expression Between Clinical Cohorts and in vitro Culture.

Background: Low frequency intrahost single nucleotide variants (iSNVs) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have been increasingly recognised as predictive indicators of positive selection. Particularly as growing numbers of SARS-CoV-2 variants of interest (VOI) and concern (VOC) emerge. However, the dynamics of subgenomic RNA (sgRNA) expression and its impact on genomic diversity and infection outcome remain poorly understood. This study aims to investigate and quantify iSNVs and sgRNA expression in single and longitudinally sampled cohorts over the course of mild and severe SARS-CoV-2 infection, benchmarked against an in vitro infection model. Methods: Two clinical cohorts of SARS-CoV-2 positive cases in New South Wales, Australia collected between March 2020 and August 2021 were sequenced. Longitudinal samples from cases hospitalised due to SARS-CoV-2 infection (severe) (n = 16) were analysed and compared with cases that presented with SARS-CoV-2 symptoms but were not hospitalised (mild) (n = 23). SARS-CoV-2 genomic diversity profiles were also examined from daily sampling of culture experiments for three SARS-CoV-2 variants (Lineage A, B.1.351, and B.1.617.2) cultured in VeroE6 C1008 cells (n = 33). Results: Intrahost single nucleotide variants were detected in 83% (19/23) of the mild cohort cases and 100% (16/16) of the severe cohort cases. SNP profiles remained relatively fixed over time, with an average of 1.66 SNPs gained or lost, and an average of 4.2 and 5.9 low frequency variants per patient were detected in severe and mild infection, respectively. sgRNA was detected in 100% (25/25) of the mild genomes and 92% (24/26) of the severe genomes. Total sgRNA expressed across all genes in the mild cohort was significantly higher than that of the severe cohort. Significantly higher expression levels were detected in the spike and the nucleocapsid genes. There was significantly less sgRNA detected in the culture dilutions than the clinical cohorts. Discussion and Conclusion: The positions and frequencies of iSNVs in the severe and mild infection cohorts were dynamic overtime, highlighting the importance of continual monitoring, particularly during community outbreaks where multiple SARS-CoV-2 variants may co-circulate. sgRNA levels can vary across patients and the overall level of sgRNA reads compared to genomic RNA can be less than 1%. The relative contribution of sgRNA to the severity of illness warrants further investigation given the level of variation between genomes. Further monitoring of sgRNAs will improve the understanding of SARS-CoV-2 evolution and the effectiveness of therapeutic and public health containment measures during the pandemic.

Co-infection with SARS-CoV-2 Omicron and Delta variants revealed by genomic surveillance.

Co-infections with different variants of SARS-CoV-2 are a key precursor to recombination events that are likely to drive SARS-CoV-2 evolution. Rapid identification of such co-infections is required to determine their frequency in the community, particularly in populations at-risk of severe COVID-19, which have already been identified as incubators for punctuated evolutionary events. However, limited data and tools are currently available to detect and characterise the SARS-CoV-2 co-infections associated with recognised variants of concern. Here we describe co-infection with the SARS-CoV-2 variants of concern Omicron and Delta in two epidemiologically unrelated adult patients with chronic kidney disease requiring maintenance haemodialysis. Both variants were co-circulating in the community at the time of detection. Genomic surveillance based on amplicon- and probe-based sequencing using short- and long-read technologies identified and quantified subpopulations of Delta and Omicron viruses in respiratory samples. These findings highlight the importance of integrated genomic surveillance in vulnerable populations and provide diagnostic pathways to recognise SARS-CoV-2 co-infection using genomic data.

Staying in place and moving in space: Contrasting larval thermal sensitivity explains distributional changes of sympatric sea urchin species to habitat warming.

For marine ectotherms, larval success, planktonic larval duration and dispersal trajectories are strongly influenced by temperature, and therefore, ocean warming and heatwaves have profound impacts on these sensitive stages. Warming, through increased poleward flow in regions with western boundary currents, such as the East Australia Current (EAC), provides opportunities for range extension as propagules track preferred conditions. Two sea urchin species, Centrostephanus rodgersii and Heliocidaris tuberculata, sympatric in the EAC warming hotspot, exhibit contrasting responses to warming. Over half a century, C. rodgersii has undergone marked poleward range extension, but the range of H. tuberculata has not changed. We constructed thermal performance curves (TPC) to determine if contrasting developmental thermal tolerance can explain this difference. The temperatures tested encompassed present-day distribution and forecast ocean warming/heatwave conditions. The broad and narrow thermal optimum (Topt) ranges for C. rodgersii and H. tuberculata larvae (7.2 and 4.7°C range, respectively) matched their realized (adult distribution) thermal niches. The cool and warm temperatures for 50% development to the feeding larva approximated temperatures at adult poleward range limits. Larval cool tolerances with respect to mean local temperature differed, 6.0 and 3.8°C respectively. Larval warm tolerances were similar for both species as are the adult warm range edges. The larvae of both species would be sensitive to heatwaves. Centrostephanus rodgersii has stayed in place and shifted in space, likely due to its broad cold-warm larval thermal tolerance and large thermal safety margins. Phenotypic plasticity of the planktonic stage of C. rodgersii facilitated its range extension. In contrast, larval cold intolerance of H. tuberculata explains its restricted range and will delay poleward extension as the region warms. In a warming ocean, we show that intrinsic thermal biology traits of the pelagic stage provide an integrative tool to explain species-specific variation in range shift patterns.

Documenting elimination of co-circulating COVID-19 clusters using genomics in New South Wales, Australia.

OBJECTIVE: To adapt 'fishplots' to describe real-time evolution of SARS-CoV-2 genomic clusters. RESULTS: This novel analysis adapted the fishplot to depict the size and duration of circulating genomic clusters over time in New South Wales, Australia. It illuminated the effectiveness of interventions on the emergence, spread and eventual elimination of clusters and distilled genomic data into clear information to inform public health action.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2 with grinch.

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Effects of marine heatwave conditions across the metamorphic transition to the juvenile sea urchin (Heliocidaris erythrogramma).

For short development species, like the sea urchin Heliocidaris erythrogramma, the entire planktonic duration can be impacted by marine heatwaves (MHW). Developmental thermal tolerance of this species through metamorphosis was investigated over a broad range (7.6-28.0 °C), including temperatures across its distribution and MHW conditions. In controls (19.5-21.0 °C), 80% of individuals developed to metamorphosis at day 5, doubling to 10 days at 14.0 °C. The thermal range (14.4-21.2 °C) of metamorphosis on day 7 reflected the realised thermal niche with 25.9 °C the upper temperature for success (T40). By day 10, juvenile tolerance narrowed to the local range (16.2-19.0 °C), similar to levels tolerated by adults, indicating negative carryover effects across the metamorphic transition. Without phenotypic adjustment or adaptation, regional warming will be detrimental, although populations may be sustained by thermotolerant offspring. Our results show the importance of the metamorphic transition in understanding the cumulative sensitivity of species to MHW.

Added Value of Genomic Surveillance of Virulence Factors in Shiga Toxin-Producing Escherichia coli in New South Wales, Australia.

The disease caused by Shiga toxin-producing Escherichia coli (STEC) remains a significant public health challenge globally, but the incidence of human STEC infections in Australia remains relatively low. This study examined the virulence characteristics and diversity of STEC isolates in the state of New South Wales between December 2017 and May 2020. Utilisation of both whole and core genome multi-locus sequence typing (MLST) allowed for the inference of genomic diversity and detection of isolates that were likely to be epidemiologically linked. The most common STEC serotype and stx subtype detected in this study were O157:H7 and stx 1a, respectively. A genomic scan of other virulence factors present in STEC suggested interplay between iron uptake system and virulence factors that mediate either iron release or countermeasures against host defence that could result in a reduction of stx 1a expression. This reduced expression of the dominant stx genotype could contribute to the reduced incidence of STEC-related illness in Australia. Genomic surveillance of STEC becomes an important part of public health response and ongoing interrogation of virulence factors in STEC offers additional insights for the public health risk assessment.

Revealing COVID-19 transmission in Australia by SARS-CoV-2 genome sequencing and agent-based modeling.

In January 2020, a novel betacoronavirus (family Coronaviridae), named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was identified as the etiological agent of a cluster of pneumonia cases occurring in Wuhan City, Hubei Province, China1,2. The disease arising from SARS-CoV-2 infection, coronavirus disease 2019 (COVID-19), subsequently spread rapidly causing a worldwide pandemic. Here we examine the added value of near real-time genome sequencing of SARS-CoV-2 in a subpopulation of infected patients during the first 10 weeks of COVID-19 containment in Australia and compare findings from genomic surveillance with predictions of a computational agent-based model (ABM). Using the Australian census data, the ABM generates over 24 million software agents representing the population of Australia, each with demographic attributes of an anonymous individual. It then simulates transmission of the disease over time, spreading from specific infection sources, using contact rates of individuals within different social contexts. We report that the prospective sequencing of SARS-CoV-2 clarified the probable source of infection in cases where epidemiological links could not be determined, significantly decreased the proportion of COVID-19 cases with contentious links, documented genomically similar cases associated with concurrent transmission in several institutions and identified previously unsuspected links. Only a quarter of sequenced cases appeared to be locally acquired and were concordant with predictions from the ABM. These high-resolution genomic data are crucial to track cases with locally acquired COVID-19 and for timely recognition of independent importations once border restrictions are lifted and trade and travel resume.

From pole to pole: the potential for the Arctic seastar Asterias amurensis to invade a warming Southern Ocean.

Due to climatic warming, Asterias amurensis, a keystone boreal predatory seastar that has established extensive invasive populations in southern Australia, is a potential high-risk invader of the sub-Antarctic and Antarctic. To assess the potential range expansion of A. amurensis to the Southern Ocean as it warms, we investigated the bioclimatic envelope of the adult and larval life stages. We analysed the distribution of adult A. amurensis with respect to present-day and future climate scenarios using habitat temperature data to construct species distribution models (SDMs). To integrate the physiological response of the dispersive phase, we determined the thermal envelope of larval development to assess their performance in present-day and future thermal regimes and the potential for success of A. amurensis in poleward latitudes. The SDM indicated that the thermal 'niche' of the adult stage correlates with a 0-17 °C and 1-22.5 °C range, in winter and summer, respectively. As the ocean warms, the range of A. amurensis in Australia will contract, while more southern latitudes will have conditions favourable for range expansion. Successful fertilization occurred from 3 to 23.8 °C. By day 12, development to the early larval stage was successful from 5.5 to 18 °C. Although embryos were able to reach the blastula stage at 2 °C, they had arrested development and high mortality. The optimal thermal range for survival of pelagic stages was 3.5-19.2 °C with a lower and upper critical limit of 2.6 and 20.3 °C, respectively. Our data predict that A. amurensis faces demise in its current invasive range while more favourable conditions at higher latitudes would facilitate invasion of both larval and adult stages to the Southern Ocean. Our results show that vigilance is needed to reduce the risk that this ecologically important Arctic carnivore may invade the Southern Ocean and Antarctica.