Exploring the role of polymorphic interspecies structural variants in reproductive isolation and adaptive divergence in Eucalyptus.

Structural variations (SVs) play a significant role in speciation and adaptation in many species, yet few studies have explored the prevalence and impact of different categories of SVs. We conducted a comparative analysis of long-read assembled reference genomes of closely related Eucalyptus species to identify candidate SVs potentially influencing speciation and adaptation. Interspecies SVs can be either fixed differences or polymorphic in one or both species. To describe SV patterns, we employed short-read whole-genome sequencing on over 600 individuals of Eucalyptus melliodora and Eucalyptus sideroxylon, along with recent high-quality genome assemblies. We aligned reads and genotyped interspecies SVs predicted between species reference genomes. Our results revealed that 49,756 of 58,025 and 39,536 of 47,064 interspecies SVs could be typed with short reads in E. melliodora and E. sideroxylon, respectively. Focusing on inversions and translocations, symmetric SVs that are readily genotyped within both populations, 24 were found to be structural divergences, 2,623 structural polymorphisms, and 928 shared structural polymorphisms. We assessed the functional significance of fixed interspecies SVs by examining differences in estimated recombination rates and genetic differentiation between species, revealing a complex history of natural selection. Shared structural polymorphisms displayed enrichment of potentially adaptive genes. Understanding how different classes of genetic mutations contribute to genetic diversity and reproductive barriers is essential for understanding how organisms enhance fitness, adapt to changing environments, and diversify. Our findings reveal the prevalence of interspecies SVs and elucidate their role in genetic differentiation, adaptive evolution, and species divergence within and between populations.

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence.

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.

Genome biology and evolution of mating-type loci in four cereal rust fungi.

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

Multiplexed effector screening for recognition by endogenous resistance genes using positive defense reporters in wheat protoplasts.

Plant resistance (R) and pathogen avirulence (Avr) gene interactions play a vital role in pathogen resistance. Efficient molecular screening tools for crops lack far behind their model organism counterparts, yet they are essential to rapidly identify agriculturally important molecular interactions that trigger host resistance. Here, we have developed a novel wheat protoplast assay that enables efficient screening of Avr/R interactions at scale. Our assay allows access to the extensive gene pool of phenotypically described R genes because it does not require the overexpression of cloned R genes. It is suitable for multiplexed Avr screening, with interactions tested in pools of up to 50 Avr candidates. We identified Avr/R-induced defense genes to create a promoter-luciferase reporter. Then, we combined this with a dual-color ratiometric reporter system that normalizes read-outs accounting for experimental variability and Avr/R-induced cell death. Moreover, we introduced a self-replicative plasmid reducing the amount of plasmid used in the assay. Our assay increases the throughput of Avr candidate screening, accelerating the study of cellular defense signaling and resistance gene identification in wheat. We anticipate that our assay will significantly accelerate Avr identification for many wheat pathogens, leading to improved genome-guided pathogen surveillance and breeding of disease-resistant crops.

Whole-mitochondrial genomes of Nannizziopsis provide insights in evolution and detection.

Infectious fungal diseases can have devastating effects on wildlife health and a detailed understanding of the evolution of related emerging fungal pathogen along with the ability to detect them in the wild is considered indispensable for effective management strategies. Several fungi from the genera Nannizziopsis and Paranannizziopsis are emerging pathogens of reptiles and have been observed to cause disease in a wide range of taxa. Nannizziopsis barbatae has become a particularly important pathogen of Australian reptiles with an increasing number of herpetofauna being reported with cases of infection from across the country. Here, we present the mitochondrial genome sequences and phylogenetic analysis for seven species in this group of fungi uncovering new information on the evolutionary relationship of these emerging pathogens. From this analysis, we designed a species-specific qPCR assay for the rapid detection of N. barbatae and demonstrate its application in a wild urban population of a dragon lizard.

Public health interventions successfully mitigated multiple incursions of SARS-CoV-2 Delta variant in the Australian Capital Territory.

The COVID-19 pandemic has presented a unique opportunity to understand how real-time pathogen genomics can be used for large-scale outbreak investigations. On 12 August 2021, the Australian Capital Territory (ACT) detected an incursion of the SARS-CoV-2 Delta (B.1.617.2) variant. Prior to this date, SARS-CoV-2 had been eliminated locally since 7 July 2020. Several public health interventions were rapidly implemented in response to the incursion, including a territory-wide lockdown and comprehensive contact tracing. The ACT has not previously used pathogen genomics at a population level in an outbreak response; therefore, this incursion also presented an opportunity to investigate the utility of genomic sequencing to support contact tracing efforts in the ACT. Sequencing of >75% of the 1793 laboratory-confirmed cases during the 3 months following the initial notification identified at least 13 independent incursions with onwards spread in the community. Stratification of cases by genomic cluster revealed that distinct cohorts were affected by the different incursions. Two incursions resulted in most of the community transmission during the study period, with persistent transmission in vulnerable sections of the community. Ultimately, both major incursions were successfully mitigated through public health interventions, including COVID-19 vaccines. The high rates of SARS-CoV-2 sequencing in the ACT and the relatively small population size facilitated detailed investigations of the patterns of virus transmission, revealing insights beyond those gathered from traditional contact tracing alone. Genomic sequencing was critical to disentangling complex transmission chains to target interventions appropriately.

Cost-conscious generation of multiplexed short-read DNA libraries for whole-genome sequencing.

Massively parallel, second-generation short-read DNA sequencing has become an integral tool in biology for genomic studies. Offering highly accurate base-pair resolution at the most competitive price, the technology has become widespread. However, high-throughput generation of multiplexed DNA libraries can be costly and cumbersome. Here, we present a cost-conscious protocol for generating multiplexed short-read DNA libraries using a bead-linked transposome from Illumina. We prepare libraries in high-throughput with small reaction volumes that use 1/50th the amount of transposome compared to Illumina DNA Prep tagmentation protocols. By reducing transposome usage and optimising the protocol to circumvent magnetic bead-based clean-ups between steps, we reduce costs, labour time and DNA input requirements. Developing our own dual index primers further reduced costs and enables up to nine 96-well microplate combinations. This facilitates efficient usage of large-scale sequencing platforms, such as the Illumina NovaSeq 6000, which offers up to three terabases of sequencing per S4 flow cell. The protocol presented substantially reduces the cost per library by approximately 1/20th compared to conventional Illumina methods.

Interspecies genome divergence is predominantly due to frequent small scale rearrangements in Eucalyptus.

Synteny, the ordering of sequences within homologous chromosomes, must be maintained within the genomes of sexually reproducing species for the sharing of alleles and production of viable, reproducing offspring. However, when the genomes of closely related species are compared, a loss of synteny is often observed. Unequal homologous recombination is the primary mechanism behind synteny loss, occurring more often in transposon rich regions, and resulting in the formation of chromosomal rearrangements. To examine patterns of synteny among three closely related, interbreeding, and wild Eucalyptus species, we assembled their genomes using long-read DNA sequencing and de novo assembly. We identify syntenic and rearranged regions between these genomes and estimate that ~48% of our genomes remain syntenic while ~36% is rearranged. We observed that rearrangements highly fragment microsynteny. Our results suggest that synteny between these species is primarily lost through small-scale rearrangements, not through sequence loss, gain, or sequence divergence. Further examination of identified rearrangements suggests that rearrangements may be altering the phenotypes of Eucalyptus species. Our study also underscores that the use of single reference genomes in genomic variation studies could lead to reference bias, especially given the scale at which we show potentially adaptive loci have highly diverged, deleted, duplicated and/or rearranged. This study provides an unbiased framework to look at potential speciation and adaptive loci among a rapidly radiating foundation species of woodland trees that are free from selective breeding seen in most crop species.

Running a research group in the next generation: combining sustainable and reproducible research with values-driven leadership.

In the summer of 2021, we held a community workshop at the International Congress of Arabidopsis Research (ICAR) aimed at early career researchers and focused on values-based lab leadership. Here, we elaborate on ideas emerging from the workshop that we hope will allow current and future group leaders to reflect on and adjust to the rapidly evolving nature of the academic scientific enterprise.

Species-specific basecallers improve actual accuracy of nanopore sequencing in plants.

BACKGROUND: Long-read sequencing platforms offered by Oxford Nanopore Technologies (ONT) allow native DNA containing epigenetic modifications to be directly sequenced, but can be limited by lower per-base accuracies. A key step post-sequencing is basecalling, the process of converting raw electrical signals produced by the sequencing device into nucleotide sequences. This is challenging as current basecallers are primarily based on mixtures of model species for training. Here we utilise both ONT PromethION and higher accuracy PacBio Sequel II HiFi sequencing on two plants, Phebalium stellatum and Xanthorrhoea johnsonii, to train species-specific basecaller models with the aim of improving per-base accuracy. We investigate sequencing accuracies achieved by ONT basecallers and assess accuracy gains by training single-species and species-specific basecaller models. We also evaluate accuracy gains from ONT's improved flowcells (R10.4, FLO-PRO112) and sequencing kits (SQK-LSK112). For the truth dataset for both model training and accuracy assessment, we developed highly accurate, contiguous diploid reference genomes with PacBio Sequel II HiFi reads. RESULTS: Basecalling with ONT Guppy 5 and 6 super-accurate gave almost identical results, attaining read accuracies of 91.96% and 94.15%. Guppy's plant-specific model gave highly mixed results, attaining read accuracies of 91.47% and 96.18%. Species-specific basecalling models improved read accuracy, attaining 93.24% and 95.16% read accuracies. R10.4 sequencing kits also improve sequencing accuracy, attaining read accuracies of 95.46% (super-accurate) and 96.87% (species-specific). CONCLUSIONS: The use of a single mixed-species basecaller model, such as ONT Guppy super-accurate, may be reducing the accuracy of nanopore sequencing, due to conflicting genome biology within the training dataset and study species. Training of single-species and genome-specific basecaller models improves read accuracy. Studies that aim to do large-scale long-read genotyping would primarily benefit from training their own basecalling models. Such studies could use sequencing accuracy gains and improving bioinformatics tools to improve study outcomes.

Recommendations for empowering early career researchers to improve research culture and practice.

Early career researchers (ECRs) are important stakeholders leading efforts to catalyze systemic change in research culture and practice. Here, we summarize the outputs from a virtual unconventional conference (unconference), which brought together 54 invited experts from 20 countries with extensive experience in ECR initiatives designed to improve the culture and practice of science. Together, we drafted 2 sets of recommendations for (1) ECRs directly involved in initiatives or activities to change research culture and practice; and (2) stakeholders who wish to support ECRs in these efforts. Importantly, these points apply to ECRs working to promote change on a systemic level, not only those improving aspects of their own work. In both sets of recommendations, we underline the importance of incentivizing and providing time and resources for systems-level science improvement activities, including ECRs in organizational decision-making processes, and working to dismantle structural barriers to participation for marginalized groups. We further highlight obstacles that ECRs face when working to promote reform, as well as proposed solutions and examples of current best practices. The abstract and recommendations for stakeholders are available in Dutch, German, Greek (abstract only), Italian, Japanese, Polish, Portuguese, Spanish, and Serbian.

The First Annotated Genome Assembly of Macrophomina tecta Associated with Charcoal Rot of Sorghum.

Charcoal rot is an important soilborne disease caused by a range of Macrophomina species, which affects a broad range of commercially important crops worldwide. Even though Macrophomina species are fungal pathogens of substantial economic importance, their mechanism of pathogenicity and host spectrum are poorly understood. There is an urgent need to better understand the biology, epidemiology, and evolution of Macrophomina species, which, in turn, will aid in improving charcoal rot management strategies. Here, we present the first high-quality genome assembly and annotation of Macrophomina tecta strain BRIP 70781 associated with charcoal rot symptoms on sorghum. Hybrid assembly integrating long reads generated by Oxford Nanopore Technology and short Illumina paired-end reads resulted in 43 contigs with a total assembly size of ∼54 Mb, and an N50 of 3.4 Mb. In total, 12,926 protein-coding genes and 7,036 repeats were predicted. Genome comparisons detected accumulation of DNA transposons in Macrophomina species associated with sorghum. The first reference genome of M. tecta generated in this study will contribute to more comparative and population genomics studies of Macrophomina species.

Inferring Species Compositions of Complex Fungal Communities from Long- and Short-Read Sequence Data.

The kingdom Fungi is highly diverse in morphology and ecosystem function. Yet fungi are challenging to characterize as they can be difficult to culture and morphologically indistinct. Overall, their description and analysis lag far behind other microbes such as bacteria. Classification of species via high-throughput sequencing is increasingly becoming the norm for pathogen detection, microbiome studies, and environmental monitoring. With the rapid development of sequencing technologies, however, standardized procedures for taxonomic assignment of long sequence reads have not yet been well established. Focusing on nanopore sequencing technology, we compared classification and community composition analysis pipelines using shotgun and amplicon sequencing data generated from mock communities comprising 43 fungal species. We show that regardless of the sequencing methodology used, the highest accuracy of species identification was achieved by sequence alignment against a fungal-specific database. During the assessment of classification algorithms, we found that applying cutoffs to the query coverage of each read or contig significantly improved the classification accuracy and community composition analysis without major data loss. We also generated draft genome assemblies for three fungal species from nanopore data which were absent from genome databases. Our study improves sequence-based classification and estimation of relative sequence abundance using real fungal community data and provides a practical guide for the design of metagenomics analyses focusing on fungi. IMPORTANCE Our study is unique in that it provides an in-depth comparative study of a real-life complex fungal community analyzed with multiple long- and short-read sequencing approaches. These technologies and their application are currently of great interest to diverse biologists as they seek to characterize the community compositions of microbiomes. Although great progress has been made on bacterial community compositions, microbial eukaryotes such as fungi clearly lag behind. Our study provides a detailed breakdown of strategies to improve species identification with immediate relevance to real-world studies. We find that real-life data sets do not always behave as expected, distinct from reports based on simulated data sets.

Physical separation of haplotypes in dikaryons allows benchmarking of phasing accuracy in Nanopore and HiFi assemblies with Hi-C data.

BACKGROUND: Most animals and plants have more than one set of chromosomes and package these haplotypes into a single nucleus within each cell. In contrast, many fungal species carry multiple haploid nuclei per cell. Rust fungi are such species with two nuclei (karyons) that contain a full set of haploid chromosomes each. The physical separation of haplotypes in dikaryons means that, unlike in diploids, Hi-C chromatin contacts between haplotypes are false-positive signals. RESULTS: We generate the first chromosome-scale, fully-phased assembly for the dikaryotic leaf rust fungus Puccinia triticina and compare Nanopore MinION and PacBio HiFi sequence-based assemblies. We show that false-positive Hi-C contacts between haplotypes are predominantly caused by phase switches rather than by collapsed regions or Hi-C read mis-mappings. We introduce a method for phasing of dikaryotic genomes into the two haplotypes using Hi-C contact graphs, including a phase switch correction step. In the HiFi assembly, relatively few phase switches occur, and these are predominantly located at haplotig boundaries and can be readily corrected. In contrast, phase switches are widespread throughout the Nanopore assembly. We show that haploid genome read coverage of 30-40 times using HiFi sequencing is required for phasing of the leaf rust genome, with 0.7% heterozygosity, and that HiFi sequencing resolves genomic regions with low heterozygosity that are otherwise collapsed in the Nanopore assembly. CONCLUSIONS: This first Hi-C based phasing pipeline for dikaryons and comparison of long-read sequencing technologies will inform future genome assembly and haplotype phasing projects in other non-haploid organisms.

A high-quality pseudo-phased genome for Melaleuca quinquenervia shows allelic diversity of NLR-type resistance genes.

BACKGROUND: Melaleuca quinquenervia (broad-leaved paperbark) is a coastal wetland tree species that serves as a foundation species in eastern Australia, Indonesia, Papua New Guinea, and New Caledonia. While extensively cultivated for its ornamental value, it has also become invasive in regions like Florida, USA. Long-lived trees face diverse pest and pathogen pressures, and plant stress responses rely on immune receptors encoded by the nucleotide-binding leucine-rich repeat (NLR) gene family. However, the comprehensive annotation of NLR encoding genes has been challenging due to their clustering arrangement on chromosomes and highly repetitive domain structure; expansion of the NLR gene family is driven largely by tandem duplication. Additionally, the allelic diversity of the NLR gene family remains largely unexplored in outcrossing tree species, as many genomes are presented in their haploid, collapsed state. RESULTS: We assembled a chromosome-level pseudo-phased genome for M. quinquenervia and described the allelic diversity of plant NLRs using the novel FindPlantNLRs pipeline. Analysis reveals variation in the number of NLR genes on each haplotype, distinct clustering patterns, and differences in the types and numbers of novel integrated domains. CONCLUSIONS: The high-quality M. quinquenervia genome assembly establishes a new framework for functional and evolutionary studies of this significant tree species. Our findings suggest that maintaining allelic diversity within the NLR gene family is crucial for enabling responses to environmental stress, particularly in long-lived plants.

Activation loop phosphorylation of a non-RD receptor kinase initiates plant innate immune signaling.

Receptor kinases (RKs) are fundamental for extracellular sensing and regulate development and stress responses across kingdoms. In plants, leucine-rich repeat receptor kinases (LRR-RKs) are primarily peptide receptors that regulate responses to myriad internal and external stimuli. Phosphorylation of LRR-RK cytoplasmic domains is among the earliest responses following ligand perception, and reciprocal transphosphorylation between a receptor and its coreceptor is thought to activate the receptor complex. Originally proposed based on characterization of the brassinosteroid receptor, the prevalence of complex activation via reciprocal transphosphorylation across the plant RK family has not been tested. Using the LRR-RK ELONGATION FACTOR TU RECEPTOR (EFR) as a model, we set out to understand the steps critical for activating RK complexes. While the EFR cytoplasmic domain is an active protein kinase in vitro and is phosphorylated in a ligand-dependent manner in vivo, catalytically deficient EFR variants are functional in antibacterial immunity. These results reveal a noncatalytic role for EFR in triggering immune signaling and indicate that reciprocal transphoshorylation is not a ubiquitous requirement for LRR-RK complex activation. Rather, our analysis of EFR along with a detailed survey of the literature suggests a distinction between LRR-RKs with RD- versus non-RD protein kinase domains. Based on newly identified phosphorylation sites that regulate the activation state of the EFR complex in vivo, we propose that LRR-RK complexes containing a non-RD protein kinase may be regulated by phosphorylation-dependent conformational changes of the ligand-binding receptor, which could initiate signaling either allosterically or through driving the dissociation of negative regulators of the complex.

The stem rust fungus Puccinia graminis f. sp. tritici induces centromeric small RNAs during late infection that are associated with genome-wide DNA methylation.

BACKGROUND: Silencing of transposable elements (TEs) is essential for maintaining genome stability. Plants use small RNAs (sRNAs) to direct DNA methylation to TEs (RNA-directed DNA methylation; RdDM). Similar mechanisms of epigenetic silencing in the fungal kingdom have remained elusive. RESULTS: We use sRNA sequencing and methylation data to gain insight into epigenetics in the dikaryotic fungus Puccinia graminis f. sp. tritici (Pgt), which causes the devastating stem rust disease on wheat. We use Hi-C data to define the Pgt centromeres and show that they are repeat-rich regions (~250 kb) that are highly diverse in sequence between haplotypes and, like in plants, are enriched for young TEs. DNA cytosine methylation is particularly active at centromeres but also associated with genome-wide control of young TE insertions. Strikingly, over 90% of Pgt sRNAs and several RNAi genes are differentially expressed during infection. Pgt induces waves of functionally diversified sRNAs during infection. The early wave sRNAs are predominantly 21 nts with a 5' uracil derived from genes. In contrast, the late wave sRNAs are mainly 22-nt sRNAs with a 5' adenine and are strongly induced from centromeric regions. TEs that overlap with late wave sRNAs are more likely to be methylated, both inside and outside the centromeres, and methylated TEs exhibit a silencing effect on nearby genes. CONCLUSIONS: We conclude that rust fungi use an epigenetic silencing pathway that might have similarity with RdDM in plants. The Pgt RNAi machinery and sRNAs are under tight temporal control throughout infection and might ensure genome stability during sporulation.

High-molecular weight DNA extraction, clean-up and size selection for long-read sequencing.

Rapid advancements in long-read sequencing technologies have transformed read lengths from bps to Mbps, which has enabled chromosome-scale genome assemblies. However, read lengths are now becoming limited by the extraction of pure high-molecular weight DNA suitable for long-read sequencing, which is particularly challenging in plants and fungi. To overcome this, we present a protocol collection; high-molecular weight DNA extraction, clean-up and size selection for long-read sequencing. We optimised a gentle magnetic bead based high-molecular weight DNA extraction, which is presented here in detail. The protocol circumvents spin columns and high-centrifugation, to limit DNA fragmentation. The protocol is scalable based on tissue input, which can be used on many species of plants, fungi, reptiles and bacteria. It is also cost effective compared to kit-based protocols and hence applicable at scale in low resource settings. An optional sorbitol wash is listed and is highly recommended for plant and fungal tissues. To further remove any remaining contaminants such as phenols and polysaccharides, optional DNA clean-up and size selection strategies are given. This protocol collection is suitable for all common long-read sequencing platforms, such as technologies offered by PacBio and Oxford Nanopore. Using these protocols, sequencing on the Oxford Nanopore MinION can achieve read length N50 values of 30-50 kb, with reads exceeding 200 kb and outputs ranging from 15-30 Gbp. This has been routinely achieved with various plant, fungi, animal and bacteria samples.

Genomics-informed responses in the elimination of COVID-19 in Victoria, Australia: an observational, genomic epidemiological study.

BACKGROUND: A cornerstone of Australia's ability to control COVID-19 has been effective border control with an extensive supervised quarantine programme. However, a rapid recrudescence of COVID-19 was observed in the state of Victoria in June, 2020. We aim to describe the genomic findings that located the source of this second wave and show the role of genomic epidemiology in the successful elimination of COVID-19 for a second time in Australia. METHODS: In this observational, genomic epidemiological study, we did genomic sequencing of all laboratory-confirmed cases of COVID-19 diagnosed in Victoria, Australia between Jan 25, 2020, and Jan 31, 2021. We did phylogenetic analyses, genomic cluster discovery, and integrated results with epidemiological data (detailed information on demographics, risk factors, and exposure) collected via interview by the Victorian Government Department of Health. Genomic transmission networks were used to group multiple genomic clusters when epidemiological and genomic data suggested they arose from a single importation event and diversified within Victoria. To identify transmission of emergent lineages between Victoria and other states or territories in Australia, all publicly available SARS-CoV-2 sequences uploaded before Feb 11, 2021, were obtained from the national sequence sharing programme AusTrakka, and epidemiological data were obtained from the submitting laboratories. We did phylodynamic analyses to estimate the growth rate, doubling time, and number of days from the first local infection to the collection of the first sequenced genome for the dominant local cluster, and compared our growth estimates to previously published estimates from a similar growth phase of lineage B.1.1.7 (also known as the Alpha variant) in the UK. FINDINGS: Between Jan 25, 2020, and Jan 31, 2021, there were 20 451 laboratory-confirmed cases of COVID-19 in Victoria, Australia, of which 15 431 were submitted for sequencing, and 11 711 met all quality control metrics and were included in our analysis. We identified 595 genomic clusters, with a median of five cases per cluster (IQR 2-11). Overall, samples from 11 503 (98·2%) of 11 711 cases clustered with another sample in Victoria, either within a genomic cluster or transmission network. Genomic analysis revealed that 10 426 cases, including 10 416 (98·4%) of 10 584 locally acquired cases, diagnosed during the second wave (between June and October, 2020) were derived from a single incursion from hotel quarantine, with the outbreak lineage (transmission network G, lineage D.2) rapidly detected in other Australian states and territories. Phylodynamic analyses indicated that the epidemic growth rate of the outbreak lineage in Victoria during the initial growth phase (samples collected between June 4 and July 9, 2020; 47·4 putative transmission events, per branch, per year [1/years; 95% credible interval 26·0-85·0]), was similar to that of other reported variants, such as B.1.1.7 in the UK (mean approximately 71·5 1/years). Strict interventions were implemented, and the outbreak lineage has not been detected in Australia since Oct 29, 2020. Subsequent cases represented independent international or interstate introductions, with limited local spread. INTERPRETATION: Our study highlights how rapid escalation of clonal outbreaks can occur from a single incursion. However, strict quarantine measures and decisive public health responses to emergent cases are effective, even with high epidemic growth rates. Real-time genomic surveillance can alter the way in which public health agencies view and respond to COVID-19 outbreaks. FUNDING: The Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.