Highly pathogenic avian influenza (HPAI) virus subtype H5N1 has never been detected in New Zealand. The potential impact of this virus on New Zealand's wild birds would be catastrophic. To expand our knowledge of avian influenza viruses across New Zealand, we sampled wild aquatic birds from New Zealand, its outer islands and its subantarctic territories. Metatranscriptomic analysis of 700 individuals spanning 33 species revealed no detection of H5N1 during the annual 2023-2024 migration. A single detection of H1N9 in red knots (Calidris canutus) was noted. This study provides a baseline for expanding avian influenza virus monitoring in New Zealand.

By early 2022, the highly transmissible Omicron variant of SARS-CoV-2 had spread across most of the world. For the first time since the pandemic began, New Zealand was experiencing high levels of community transmission of SARS-CoV-2. We enroled a cohort of households to better understand differences in transmission dynamics among subvariants of Omicron. We enroled 71 households, comprising 289 participants, and aimed to use viral genomes to gain a clearer understanding of variant-specific differences in epidemiological parameters affecting transmission dynamics. Approximately 80% of the households enroled experienced transmission of BA.2, while most of the remaining households had infections with BA.1 or BA.5. Using a logistic regression generalised linear mixed model, we found no difference in household secondary infection rate between Omicron subvariants BA.1, BA.2 and BA.5. Of the households recruited, the vast majority (92%) experienced a single chain of transmission with one inferred introduction. Further, we found that in 48% of the households studied, all household participants became infected following an index case. Most household participants tested positive within a week following an introduction, supporting the seven-day isolation requirement for household contacts that was in place in New Zealand at the time. By integrating genomic and epidemiological data, we show that viral transmission dynamics can be investigated with a higher level of granularity than with epidemiological data alone. Overall, households are a high risk setting for viral transmission in New Zealand.

Summary Leaf mould, caused by Fulvia fulva, is a devastating disease of tomato plants. In many commercial tomato cultivars, resistance to this disease is governed by the Cf‐9 locus, which encodes five paralogous receptor‐like proteins. Two of these proteins confer resistance: Cf‐9C recognises the previously identified F. fulva effector Avr9 and provides resistance during all plant growth stages, while Cf‐9B recognises the yet‐unidentified F. fulva effector Avr9B and provides mature plant resistance only. In recent years, F. fulva strains have emerged that can overcome the Cf‐9 locus, with Cf‐9C circumvented through Avr9 deletion. To understand how Cf‐9B is circumvented, we set out to identify Avr9B. Comparative genomics, transient expression assays and gene complementation experiments were used to identify Avr9B, while gene sequencing was used to assess Avr9B allelic variation across a world‐wide strain collection. A strict correlation between Avr9 deletion and resistance‐breaking mutations in Avr9B was observed in strains recently collected from Cf‐9 cultivars, whereas Avr9 deletion but no mutations in Avr9B were observed in older strains. This research showcases how F. fulva has evolved to sequentially break down the Cf‐9 locus and stresses the urgent need for commercial tomato cultivars that carry novel, stacked resistance genes active against this pathogen.

Background Genomic technologies have become routine in the surveillance and monitoring of the coronavirus disease 2019 (COVID-19) pandemic, as evidenced by the millions of SARS-CoV-2 sequences uploaded to international databases. Yet the ways in which these technologies have been applied to manage the pandemic are varied. Main text Aotearoa New Zealand was one of a small number of countries to adopt an elimination strategy for COVID-19, establishing a managed isolation and quarantine system for all international arrivals. To aid our response, we rapidly set up and scaled our use of genomic technologies to help identify community cases of COVID-19, to understand how they had arisen, and to determine the appropriate action to maintain elimination. Once New Zealand pivoted from elimination to suppression in late 2021, our genomic response changed to focusing on identifying new variants arriving at the border, tracking their incidence around the country, and examining any links between specific variants and increased disease severity. Wastewater detection, quantitation and variant detection were also phased into the response. Here, we explore New Zealand’s genomic journey through the pandemic and provide a high-level overview of the lessons learned and potential future capabilities to better prepare for future pandemics. Conclusions Our commentary is aimed at health professionals and decision-makers who might not be familiar with genetic technologies, how they can be used, and why this is an area with great potential to assist in disease detection and tracking now and in the future.

Summary Leaf mould, caused by Fulvia fulva , is a devastating disease of tomato plants. In many commercial tomato cultivars, resistance to this disease is governed by the Cf-9 locus, which comprises five paralogous genes ( Cf-9A–9E ) that encode receptor-like proteins. Two of these proteins contribute to resistance: Cf-9C recognizes the previously identified F. fulva effector Avr9 and provides resistance during all plant growth stages, while Cf-9B recognises the yet-unidentified F. fulva effector Avr9B and provides mature plant resistance only. In recent years, F. fulva strains have emerged that have overcome the Cf-9 locus, with Cf-9C circumvented through Avr9 deletion. To understand how Cf-9B is circumvented, we set out to identify Avr9B . Comparative genomics, in planta transient expression assays and gene complementation experiments were used to identify Avr9B , while gene sequencing was used to assess Avr9B allelic variation across a worldwide strain collection. A strict correlation between Avr9 deletion and resistance-breaking mutations in Avr9B was observed in strains recently collected from Cf-9 cultivars, whereas Avr9 deletion but no mutations in Avr9B were observed in older strains. This research showcases how F. fulva has evolved to sequentially break down the two functional resistance genes of the complex Cf-9 locus and highlights that this locus now has limited value for controlling leaf mould disease in worldwide commercial tomato production.

The kākāpō is a critically endangered, intensively managed, long-lived nocturnal parrot endemic to Aotearoa New Zealand. We generated and analyzed whole-genome sequence data for nearly all individuals living in early 2018 (169 individuals) to generate a high-quality species-wide genetic variant callset. We leverage extensive long-term metadata to quantify genome-wide diversity of the species over time and present new approaches using probabilistic programming, combined with a phenotype dataset spanning five decades, to disentangle phenotypic variance into environmental and genetic effects while quantifying uncertainty in small populations. We find associations for growth, disease susceptibility, clutch size, and egg fertility within genic regions previously shown to influence these traits in other species. Finally, we generate breeding values to predict phenotype and illustrate that active management over the past 45 years has maintained both genome-wide diversity and diversity in breeding values, and hence, evolutionary potential. We provide new pathways for informing future conservation management decisions for kākāpō, including prioritizing individuals for translocation and monitoring individuals with poor growth or high disease risk. Overall, by explicitly addressing the challenge of small sample size, we provide a template for the inclusion of genomic data that will be transformational for species recovery efforts around the globe.

Genome rearrangements in filamentous fungi are prevalent but little is known about the modalities of their evolution, in part because few complete genomes are available within a single genus. To address this, we have generated and compared 15 complete telomere-to-telomere genomes across the phylogeny of a single genus of filamentous fungi, Epichloë. We find that the striking distinction between gene-rich and repeat-rich regions previously reported for isolated species is ubiquitous across the Epichloë genus. We built a species phylogeny from single-copy gene orthologs to provide a comparative framing to study chromosome composition and structural change through evolutionary time. All Epichloë genomes have exactly seven nuclear chromosomes, but despite this conserved ploidy, analyses reveal low synteny and substantial rearrangement of gene content across the genus. These rearrangements are highly lineage-dependent, with most occurring over short evolutionary distances, with long periods of structural stasis. Quantification of chromosomal rearrangements shows they are uncorrelated with numbers of substitutions and evolutionary distances, suggesting that different modes of evolution are acting to create nucleotide and chromosome-scale changes.

Hybridization is a route to speciation that occurs widely across the eukaryote tree of life. The success of allopolyploids (hybrid species with increased ploidy) and homoploid hybrids (with unchanged ploidy) is well documented. However, their formation and establishment is not straightforward, with a suite of near-instantaneous and longer term biological repercussions faced by the new species. Central to these challenges is the rewiring of gene regulatory networks following the merger of distinct genomes inherited from both parental species. Research on the evolution of hybrid gene expression has largely involved studies on a single hybrid species or a few gene families. Here, we present the first standardized transcriptome-wide study exploring the fates of genes following hybridization across three kingdoms: animals, plants and fungi. Within each kingdom, we pair an allopolyploid system with a closely related homoploid hybrid to decouple the influence of increased ploidy from genome merger. Genome merger, not changes in ploidy, has the greatest effect on posthybridization expression patterns across all study systems. Strikingly, we find that differentially expressed genes in parent species preferentially switch to more similar expression in hybrids across all kingdoms, likely as a consequence of regulatory trans-acting cross-talk within the hybrid nucleus. We also highlight the prevalence of gene loss or silencing among extremely differentially expressed genes in hybrid species across all kingdoms. These shared patterns suggest that the evolutionary process of hybridization leads to common high-level expression outcomes, regardless of the particular species or kingdom.

In the second quarter of 2022, there was a global surge of emergent SARS-CoV-2 lineages that had a distinct growth advantage over then-dominant Omicron BA.1 and BA.2 lineages. By generating 10,403 Omicron genomes, we show that Aotearoa New Zealand observed an influx of these immune-evasive variants (BA.2.12.1, BA.4, and BA.5) through the border. This is explained by the return to significant levels of international travel following the border's reopening in March 2022. We estimate one Omicron transmission event from the border to the community for every ~5,000 passenger arrivals at the current levels of travel and restriction. Although most of these introductions did not instigate any detected onward transmission, a small minority triggered large outbreaks. Genomic surveillance at the border provides a lens on the rate at which new variants might gain a foothold and trigger new waves of infection.

New Zealand's COVID-19 elimination strategy heavily relied on the use of genomics to inform contact tracing, linking cases to the border and to clusters during community outbreaks. In August 2021, New Zealand entered its second nationwide lockdown after the detection of a single community case with no immediately apparent epidemiological link to the border. This incursion resulted in the largest outbreak seen in New Zealand caused by the Delta Variant of Concern. Here we generated 3806 high quality SARS-CoV-2 genomes from cases reported in New Zealand between 17 August and 1 December 2021, representing 43% of reported cases. We detected wide geographical spread coupled with undetected community transmission, characterised by the apparent extinction and reappearance of genomically linked clusters. We also identified the emergence, and near replacement, of genomes possessing a 10nucleotide frameshift deletion that caused the likely truncation of accessory protein ORF7a. By early October, New Zealand moved from an elimination strategy to a suppression strategy and the role of genomics changed markedly from being used to track and trace, towards population-level surveillance.

New Zealand’s COVID-19 elimination strategy heavily relied on the use of genomics to inform contact tracing, linking cases to the border and to clusters during community outbreaks. In August 2021, New Zealand entered its second nationwide lockdown after the detection of a single community case with no immediately apparent epidemiological link to the border. This incursion resulted in the largest outbreak seen in New Zealand caused by the Delta Variant of Concern. Here we generated 3806 high quality SARS-CoV-2 genomes from cases reported in New Zealand between 17 August and 1 December 2021, representing 43% of reported cases. We detected wide geographical spread coupled with undetected community transmission, characterised by the apparent extinction and reappearance of genomically linked clusters. We also identified the emergence, and near replacement, of genomes possessing a 10-nucleotide frameshift deletion that caused the likely truncation of accessory protein ORF7a. By early October, New Zealand moved from elimination to suppression and the role of genomics changed markedly from being used to track and trace, towards population-level surveillance.

Abstract Fungi from the genus Epichloë form systemic endobiotic infections of cool season grasses, producing a range of host-protective natural products in return for access to nutrients. These infections are asymptomatic during vegetative host growth, with associations between asexual Epichloë spp. and their hosts considered mutualistic. However, the sexual cycle of Epichloë spp. involves virulent growth, characterized by the envelopment and sterilization of a developing host inflorescence by a dense sheath of mycelia known as a stroma. Microscopic analysis of stromata revealed a dramatic increase in hyphal propagation and host degradation compared with asymptomatic tissues. RNAseq was used to identify differentially expressed genes in asymptomatic vs stromatized tissues from 3 diverse Epichloë–host associations. Comparative analysis identified a core set of 135 differentially expressed genes that exhibited conserved transcriptional changes across all 3 associations. The core differentially expressed genes more strongly expressed during virulent growth encode proteins associated with host suppression, digestion, adaptation to the external environment, a biosynthetic gene cluster, and 5 transcription factors that may regulate Epichloë stroma formation. An additional 5 transcription factor encoding differentially expressed genes were suppressed during virulent growth, suggesting they regulate mutualistic processes. Expression of biosynthetic gene clusters for natural products that suppress herbivory was universally suppressed during virulent growth, and additional biosynthetic gene clusters that may encode production of novel host-protective natural products were identified. A comparative analysis of 26 Epichloë genomes found a general decrease in core differentially expressed gene conservation among asexual species, and a specific decrease in conservation for the biosynthetic gene cluster expressed during virulent growth and an unusual uncharacterized gene.

Legionella longbeachae is an environmental bacterium that is the most clinically significant Legionella species in New Zealand (NZ), causing around two-thirds of all notified cases of Legionnaires’ disease. Here we report the sequencing and analysis of the geo-temporal genetic diversity of 54 L. longbeachae serogroup 1 (sg1) clinical isolates, derived from cases from around NZ over a 22-year period, including one complete genome and its associated methylome. The 54 sg1 isolates belonged to two main clades that last shared a common ancestor between 95 BCE and 1694 CE. There was diversity at the genome-structural level, with large-scale arrangements occurring in some regions of the chromosome and evidence of extensive chromosomal and plasmid recombination. This includes the presence of plasmids derived from recombination and horizontal gene transfer between various Legionella species, indicating there has been both intra- and inter-species gene flow. However, because similar plasmids were found among isolates within each clade, plasmid recombination events may pre-empt the emergence of new L. longbeachae strains. Our complete NZ reference genome consisted of a 4.1 Mb chromosome and a 108 kb plasmid. The genome was highly methylated with two known epigenetic modifications, m 4 C and m 6 A, occurring in particular sequence motifs within the genome.

Kākāpō are a critically endangered species of parrots restricted to a few islands off the coast of New Zealand. Kākāpō are very closely monitored, especially during nesting seasons. In 2019, during a highly successful nesting season, an outbreak of aspergillosis affected 21 individuals and led to the deaths of 9, leaving a population of only 211 kākāpō. In monitoring this outbreak, cultures of aspergillus were grown, and genome sequenced. These sequences demonstrate that, very unusually for an aspergillus outbreak, a single strain of aspergillus caused the outbreak. This strain was found on two islands, but only one had an outbreak of aspergillosis; indicating that the strain was necessary, but not sufficient, to cause disease. Our analysis provides an understanding of the 2019 outbreak and provides potential ways to manage such events in the future.

Phytophthora species are notorious plant pathogens, with some causing devastating tree diseases that threaten the survival of their host species. One such example is Phytophthora agathidicida, the causal agent of kauri dieback - a root and trunk rot disease that kills the ancient, iconic and culturally significant tree species, Agathis australis (New Zealand kauri). A deeper understanding of how Phytophthora pathogens infect their hosts and cause disease is critical for the development of effective treatments. Such an understanding can be gained by interrogating pathogen genomes for effector genes, which are involved in virulence or pathogenicity. Although genome sequencing has become more affordable, the complete assembly of Phytophthora genomes has been problematic, particularly for those with a high abundance of repetitive sequences. Therefore, effector genes located in repetitive regions could be truncated or missed in a fragmented genome assembly. Using a combination of long-read PacBio sequences, chromatin conformation capture (Hi-C) and Illumina short reads, we assembled the P. agathidicida genome into ten complete chromosomes, with a genome size of 57 Mb including 34% repeats. This is the first Phytophthora genome assembled to chromosome level and it reveals a high level of syntenic conservation with the complete genome of Peronospora effusa, the only other completely assembled genome sequence of an oomycete. All P. agathidicida chromosomes have clearly defined centromeres and contain candidate effector genes such as RXLRs and CRNs, but in different proportions, reflecting the presence of gene family clusters. Candidate effector genes are predominantly found in gene-poor, repeat-rich regions of the genome, and in some cases showed a high degree of duplication. Analysis of candidate RXLR effector genes that occur in multicopy gene families indicated half of them were not expressed in planta. Candidate CRN effector gene families showed evidence of transposon-mediated recombination leading to new combinations of protein domains, both within and between chromosomes. Further analysis of this complete genome assembly will help inform new methods of disease control against P. agathidicida and other Phytophthora species, ultimately helping decipher how Phytophthora pathogens have evolved to shape their effector repertoires and how they might adapt in the future.

Epiplakin, a cytolinker protein expressed in the epithelium that interacts with intermediate filaments [ [1] Fujiwara S. Takeo N. Otani Y. Parry D.A.D. Kunimatsu M. Lu R. et al. Epiplakin, a novel member of the plakin family originally identified as a 450-kDa human epidermal autoantigen : Structure and tissue localization. J. Biol. Chem. 2001; 276: 13340-13347 Abstract Full Text Full Text PDF PubMed Scopus (46) Google Scholar , [2] Wang W. Sumiyoshi H. Yoshioka H. Fujiwara S. Interactions between epiplakin and intermediate filaments. J. Dermatol. 2006; 33: 518-527 Crossref PubMed Google Scholar ], contains domains that are homologous to the B domain in the carboxy-terminal region of desmoplakin. The last five B domains in human epiplakin, together with their associated linker regions, are highly conserved. Our previous study showed that there is copy number variation (CNV) in the 3′ terminal half of the epiplakin gene (EPPK1) in human individuals and that CNV reflects a variation of molecular size. Indeed, epiplakin isoforms in humans contain between 11 and 17 B domains, although 13 B domains predominate [ [3] Ishikawa K. Furuhashi M. Sasaki T. Kudoh J. Tsuchisaka A. Hashimoto T. et al. Intragenic copy number variation within human epiplakin 1 (EPPK1) generates variation of molecular size of epiplakin. J. Dermatol. Sci. 2018; 91: 228-231 Abstract Full Text Full Text PDF Google Scholar ].

On the basis of sequence homology with mammalian α-keratins, and on the criteria that the coiled-coil segments and central linker in the rod domain of these molecules must have conserved lengths if they are to assemble into viable intermediate filaments, a total of 28 Type I and Type II keratin intermediate filament chains (KIF) have been identified from the genome of the European common wall lizard (Podarcis muralis). Using the same criteria this number may be compared to 33 found here in the green anole lizard (Anole carolinensis) and 25 in the tuatara (Sphenodon punctatus). The Type I and Type II KIF genes in the wall lizard fall in clusters on chromosomes 13 and 2 respectively. Although some differences occur in the terminal domains in the KIF chains of the two lizards and tuatara, the similarities between key indicator residues - cysteine, glycine and proline - are significant. The terminal domains of the KIF chains in the wall lizard also contain sequence repeats commonly based on glycine and large apolar residues and would permit the fine tuning of physical properties when incorporated within the intermediate filaments. The H1 domain in the Type II chain is conserved across the lizards, tuatara and mammals, and has been related to its role in assembly at the 2-4 molecule level. A KIF-like chain (K80) with an extensive tail domain comprised of multiple tandem repeats has been identified as having a potential filament-crosslinking role.

Epichloe fungi are endophytes of cool season grasses, both wild species and commercial cultivars, where they may exhibit mutualistic or pathogenic lifestyles. The Epichloe-grass symbiosis is of great interest to agricultural research for the fungal bioprotective properties conferred to host grasses but also serves as an ideal system to study the evolution of fungal plant-pathogens in natural environments. Here, we assembled and annotated gapless chromosome-level genomes of two pathogenic Epichloe sibling species. Both genomes have a bipartite genome organization, with blocks of highly syntenic gene-rich regions separated by blocks of AT-rich DNA. The AT-rich regions show an extensive signature of RIP (repeat-induced point mutation) and the expansion of this compartment accounts for the large difference in genome size between the two species. This study reveals how the rapid evolution of repeat structure can drive divergence between closely related taxa and highlights the evolutionary role of dynamic compartments in fungal genomes.

Determination of the sequences of the keratin intermediate filament chains in tuatara has shown that these are closely akin to the a-keratins in human and other vertebrates, especially in the central, coiled-coil rod region. The domain lengths within the rod are preserved exactly, both Type I and Type II chains have been recognised, and sequence identity/homology exists between their respective chains. Nonetheless, there are characteristic differences in amino acid composition and sequence between their respective head (N-terminal) domains and their tail (C-terminal) domains, though some similarities are retained. Further, there is evidence of tandem repeats of a variety of lengths in the tuatara heads and tails indicative of sequence duplication events. These are not evident in human a-keratins and would therefore have implications for the physical attributes of the tissues in the two species. Multiple families of keratin-associated proteins that are ultra-high cysteine-rich or glycine + tyrosinerich in human and other species do not have direct equivalents in the tuatara. Although high-sulphur proteins are present in tuatara the cysteine residue contents are significantly lower than in human. Further, no sequence homologies between the HS proteins in the two species have been found, thereby casting considerable doubt as to whether any matrix-forming high-sulphur proteins exist in tuatara. These observations may be correlated with the numerous cysteine-rich beta-keratins (corneous beta-proteins) that are present in tuatara, but which are conspicuously absent in mammals.