A population genomics analysis of the Aotearoa New Zealand endemic rewarewa tree (Knightia excelsa).

Rewarewa (Knightia excelsa, Proteaceae) is a tree species endemic to Aotearoa New Zealand, with a natural distribution spanning Te Ika-a-Maui (North Island) and the top of Te Waipounamu (South Island). We used the pseudo-chromosome genome assembly of rewarewa as a reference and whole genome pooled sequencing from 35 populations sampled across Aotearoa New Zealand, including trees growing on Maori-owned land, to identify 1,443,255 single nucleotide polymorphisms (SNPs). Four genetic clusters located in the northern North Island (NNI), eastern North Island (NIE), western and southern North Island (NIWS), and the South Island (SI) were identified. Gene flow was revealed between the SI and NIE genetic clusters, plus bottleneck and contraction events within the genetic clusters since the mid-late Pleistocene, with divergence between North and South Island clusters estimated to have occurred ~115,000-230,000 years ago. Genotype environment analysis (GEA) was used to identify loci and genes linked with altitude, soil pH, soil carbon, slope, soil size, annual mean temperature, mean diurnal range, isothermality, annual precipitation, and precipitation seasonality. The location of the SNPs associated with these environmental variables was compared with the position of 52,192 gene-coding sequences that were predicted in the rewarewa genome using RNA sequencing. This new understanding of the genetic variation present in rewarewa and insights into the genetic control of adaptive traits will inform efforts to incorporate the species in restoration plantings and for marketing rewarewa honey based on provenance.

A Review of Viruses Infecting Yam (Dioscorea spp.).

Yam is an important food staple for millions of people globally, particularly those in the developing countries of West Africa and the Pacific Islands. To sustain the growing population, yam production must be increased amidst the many biotic and abiotic stresses. Plant viruses are among the most detrimental of plant pathogens and have caused great losses of crop yield and quality, including those of yam. Knowledge and understanding of virus biology and ecology are important for the development of diagnostic tools and disease management strategies to combat the spread of yam-infecting viruses. This review aims to highlight current knowledge on key yam-infecting viruses by examining their characteristics, genetic diversity, disease symptoms, diagnostics, and elimination to provide a synopsis for consideration in developing diagnostic strategy and disease management for yam.

Genome-wide patterns of genetic diversity, population structure and demographic history in manuka (Leptospermum scoparium) growing on indigenous Maori land.

Leptospermum scoparium J. R. Forst et G. Forst, known as manuka by Maori, the indigenous people of Aotearoa (New Zealand), is a culturally and economically significant shrub species, native to New Zealand and Australia. Chemical, morphological and phylogenetic studies have indicated geographical variation of manuka across its range in New Zealand, and genetic differentiation between New Zealand and Australia. We used pooled whole genome re-sequencing of 76 L. scoparium and outgroup populations from New Zealand and Australia to compile a dataset totalling ~2.5 million SNPs. We explored the genetic structure and relatedness of L. scoparium across New Zealand, and between populations in New Zealand and Australia, as well as the complex demographic history of this species. Our population genomic investigation suggests there are five geographically distinct manuka gene pools within New Zealand, with evidence of gene flow occurring between these pools. Demographic modelling suggests three of these gene pools have undergone expansion events, whilst the evolutionary histories of the remaining two have been subjected to contractions. Furthermore, manuka populations in New Zealand are genetically distinct from populations in Australia, with coalescent modelling suggesting these two clades diverged ~9-12 million years ago. We discuss the evolutionary history of this species and the benefits of using pool-seq for such studies. Our research will support the management and conservation of manuka by landowners, particularly Maori, and the development of a provenance story for the branding of manuka based products.

An exploration of assembly strategies and quality metrics on the accuracy of the rewarewa (Knightia excelsa) genome.

We used long read sequencing data generated from Knightia excelsa, a nectar-producing Proteaceae tree endemic to Aotearoa (New Zealand), to explore how sequencing data type, volume and workflows can impact final assembly accuracy and chromosome reconstruction. Establishing a high-quality genome for this species has specific cultural importance to Maori and commercial importance to honey producers in Aotearoa. Assemblies were produced by five long read assemblers using data subsampled based on read lengths, two polishing strategies and two Hi-C mapping methods. Our results from subsampling the data by read length showed that each assembler tested performed differently depending on the coverage and the read length of the data. Subsampling highlighted that input data with longer read lengths but perhaps lower coverage constructed more contiguous, kmers and gene-complete assemblies than short read length input data with higher coverage. The final genome assembly was constructed into 14 pseudochromosomes using an initial flye long read assembly, a racon/medaka/pilon combined polishing strategy, salsa2 and allhic scaffolding, juicebox curation, and Macadamia linkage map validation. We highlighted the importance of developing assembly workflows based on the volume and read length of sequencing data and established a robust set of quality metrics for generating high-quality assemblies. Scaffolding analyses highlighted that problems found in the initial assemblies could not be resolved accurately by Hi-C data and that assembly scaffolding was more successful when the underlying contig assembly was of higher accuracy. These findings provide insight into how quality assessment tools can be implemented throughout genome assembly pipelines to inform the de novo reconstruction of a high-quality genome assembly for nonmodel organisms.

Resistance of New Zealand Provenance Leptospermum scoparium, Kunzea robusta, Kunzea linearis, and Metrosideros excelsa to Austropuccinia psidii.

Resistance to the pandemic strain of Austropuccinia psidii was identified in New Zealand provenance Leptospermum scoparium, Kunzea robusta, and K. linearis plants. Only 1 Metrosideros excelsa-resistant plant was found (of the 570 tested) and no resistant plants of either Lophomyrtus bullata or L. obcordata were found. Three types of resistance were identified in Leptospermum scoparium. The first two, a putative immune response and a hypersensitive response, are leaf resistance mechanisms found in other myrtaceous species while on the lateral and main stems a putative immune stem resistance was also observed. Both leaf and stem infection were found on K. robusta and K. linearis plants as well as branch tip dieback that developed on almost 50% of the plants. L. scoparium, K. robusta, and K. linearis are the first myrtaceous species where consistent infection of stems has been observed in artificial inoculation trials. This new finding and the first observation of significant branch tip dieback of plants of the two Kunzea spp. resulted in the development of two new myrtle rust disease severity assessment scales. Significant seed family and provenance effects were found in L. scoparium, K. robusta, and K. linearis: some families produced significantly more plants with leaf, stem, and (in Kunzea spp.) branch tip dieback resistance, and provenances provided different percentages of resistant families and plants. The distribution of the disease symptoms on plants from the same seed family, and between plants from different seed families, suggested that the leaf, stem, and branch tip dieback resistances were the result of independent disease resistance mechanisms.

Molecular markers for Nassella trichotoma (Poaceae) to study genetic variation in New Zealand.

PREMISE OF THE STUDY: Simple sequence repeat (SSR) markers were developed for the study of genetic diversity of New Zealand Nassella trichotoma (Poaceae) and to support future studies in its native range. METHODS AND RESULTS: Genomic DNA was extracted from N. trichotoma leaf material and subjected to Roche 454 sequencing. From a total of 745 putative SSRs, 48 with di- to pentanucleotide repeats were screened, 32 primer pairs were designed, and 15 polymorphic markers were optimized for multiplex PCR on 105 N. trichotoma samples from four New Zealand regions. Each locus resulted in two to six alleles per locus, and four of the loci cross-amplified in N. tenuissima. The mean observed and expected heterozygosity ranged from 0.00 to 0.90 and 0.00 to 0.50 per locus, respectively. CONCLUSIONS: The novel SSR markers are valuable for the study of genetic diversity of N. trichotoma and might also be useful for closely related species.

Genetic Variation and Hybridisation among Eight Species of kowhai (Sophora: Fabaceae) from New Zealand Revealed by Microsatellite Markers.

We analysed nine microsatellite markers for 626 individuals representing the geographic range of eight closely related endemic New Zealand species of Sophora. Structure analysis identified the optimal K value as seven, with samples identified as Sophorachathamica, Sophorafulvida, Sophoralongicarinata, and Sophoraprostrata retrieved as well-defined groups. The remaining samples formed less resolved groups referable to Sophoratetraptera and Sophoragodleyi, with Sophoramicrophylla and Sophoramolloyi forming the seventh group. Our data suggest that considerable admixture occurs and this is most likely the result of hybridisation or introgression. S.fulvida shows admixture with the sympatric S.chathamica, and the widespread S.microphylla exhibits admixture with the sympatric S.godleyi, S.molloyi, and S.tetraptera.

A plant growth-promoting symbiosis between Mycena galopus and Vaccinium corymbosum seedlings.

Typically, Mycena species are viewed as saprotrophic fungi. However, numerous detections of Mycena spp. in the roots of green plants suggest that a continuum from saprotrophy to biotrophy could exist. In particular, mycenoid species have repeatedly been found in Ericaceae plant roots. Our study asked whether (1) Mycena species are commonly found in the roots of green Ericaceae plants; (2) Mycena sequences are limited to a single group/lineage within the genus; and (3) a Mycena sp. can behave as a beneficial root associate with a typical ericoid mycorrhizal plant (Vaccinium corymbosum), regardless of how much external labile carbon is available. We detected Mycena sequences in roots of all sampled Ericaceae plants. Our Mycena sequences clustered in four different groups distributed across the Mycena genus. Only one group could be assigned with confidence to a named species (M. galopus). Our Mycena sequences clustered with other Mycena sequences detected in roots of ericoid mycorrhizal plant species collected throughout Europe, America, and Australia. An isolate of M. galopus promoted growth of V. corymbosum seedlings in vitro regardless of external carbon supply in the media. Seedlings inoculated with M. galopus grew as well as those inoculated with the ericoid mycorrhizal fungus Rhizoscyphus ericae. Surprisingly, this M. galopus isolate colonized Vaccinium roots and formed distinctive peg-like structures. Our results suggest that Mycena species might operate along a saprotroph-symbiotic continuum with a range of ericoid mycorrhizal plant species. We discuss our results in terms of fungal partner recruitment by Ericaceae plants.

Simple sequence repeat markers for kanuka (Kunzea spp.; Myrtaceae) present in New Zealand.

PREMISE OF THE STUDY: We developed simple sequence repeat (SSR) markers to facilitate population genetic studies on kanuka (Kunzea spp.; Myrtaceae). METHODS AND RESULTS: A shotgun sequencing library was constructed from leaf material of K. robusta using a Roche 454 Junior sequencer, and a total of 3174 putative SSR regions were identified. Sixteen polymorphic markers were optimized for multiplex PCR on 10 endemic New Zealand Kunzea species. Each of these loci cross-amplified in all tested species. The amplified di-, tri-, and pentanucleotide repeats resulted in eight to 24 alleles per locus for a total of 220 specimens. The mean observed and expected heterozygosity per locus ranged from 0.18 to 0.77 and 0.33 to 0.82, respectively. CONCLUSIONS: The SSR markers we produced are valuable for phylogenetic and population studies on all endemic Kunzea spp. and may also be useful for studies on closely related Kunzea species from Australia.

Simple sequence repeat markers for the endangered species Clianthus puniceus and C. maximus (Fabaceae).

PREMISE OF THE STUDY: Microsatellite markers were developed for Clianthus puniceus using a shotgun sequencing library and tested for cross amplification in the closely related C. maximus to inform population management of these two endangered species. • METHODS AND RESULTS: We constructed a shotgun sequencing library using a Roche 454 sequencer and searched the resulting data set for putative microsatellite regions. We optimized 12 of these regions to produce polymorphic markers for Clianthus. We tested these markers on four populations of C. maximus and on four C. puniceus individuals of known provenance. Alleles per locus ranged from two to nine, while observed and expected heterozygosities per locus ranged from 0.000 to 1.000 and 0.178 to 0.600, respectively. • CONCLUSIONS: These markers will be valuable for ongoing monitoring of the genetic variation in naturally occurring populations of Clianthus and for the selection of individuals for revegetation projects in the species' former range.

Sophora microphylla (Fabaceae) microsatellite markers and their utility across the genus.

PREMISE OF THE STUDY: Genus-specific microsatellite markers were developed for Sophora for population genetic and systematic studies of the group in New Zealand, and potentially elsewhere in the geographic range. • METHODS AND RESULTS: From sequencing a total genomic DNA library (using Roche 454), we identified and developed 29 polymorphic microsatellite markers for S. microphylla and S. chathamica. We tested 12 of these markers on 14 S. chathamica individuals and four S. microphylla populations. All loci amplified in both species and species-specific alleles occurred at seven loci. In S. microphylla populations, the observed and expected heterozygosities ranged from 0.000-0.960 and 0.000-0.908, respectively, with alleles per locus ranging from seven to 23. • CONCLUSIONS: The developed markers will be valuable in studies of phylogenetics, population structure, mating system, and selection of provenances for restoration projects.

Genome size variation and evolution in allotetraploid Arabidopsis kamchatica and its parents, Arabidopsis lyrata and Arabidopsis halleri.

Polyploidization and subsequent changes in genome size are fundamental processes in evolution and diversification. Little is currently known about the extent of genome size variation within taxa and the evolutionary forces acting on this variation. Arabidopsis kamchatica has been reported to contain both diploid and tetraploid individuals. The aim of this study was to determine the genome size of A. kamchatica, whether there is variation in ploidy and/or genome size in A. kamchatica and to study how genome size has evolved. We used propidium iodide flow cytometry to measure 2C DNA content of 73 plants from 25 geographically diverse populations of the putative allotetraploid A. kamchatica and its parents, Arabidopsis lyrata and Arabidopsis halleri. All A. kamchatica plants appear to be tetraploids. The mean 2C DNA content of A. kamchatica was 1.034 pg (1011 Mbp), which is slightly smaller than the sum of its diploid parents (A. lyrata: 0.502 pg; A. halleri: 0.571 pg). Arabidopsis kamchatica appears to have lost ∼37.594 Mbp (3.6 %) of DNA from its 2C genome. Tetraploid A. lyrata from Germany and Austria appears to have lost ∼70.366 Mbp (7.2 %) of DNA from the 2C genome, possibly due to hybridization with A. arenosa, which has a smaller genome than A. lyrata. We did find genome size differences among A. kamchatica populations, which varied up to 7 %. Arabidopsis kamchatica ssp. kawasakiana from Japan appears to have a slightly larger genome than A. kamchatica ssp. kamchatica from North America, perhaps due to multiple allopolyploid origins or hybridization with A. halleri. However, the among-population coefficient of variation in 2C DNA content is lower in A. kamchatica than in other Arabidopsis taxa. Due to its close relationship to A. thaliana, A. kamchatica has the potential to be very useful in the study of polyploidy and genome evolution.

Cost-effective microsatellite markers for Banksia integrifolia (Proteaceae).

PREMISE OF THE STUDY: Microsatellite markers were developed to assess the level of genetic variation and population structure in Banksia integrifolia, a widespread species endemic to eastern Australia. • METHODS AND RESULTS: We used next-generation sequencing approaches to identify and develop 11 polymorphic microsatellite markers with perfect tri- and tetranucleotide repeats. We tested these markers with 71 specimens from three populations. Observed and expected heterozygosities ranged from 0.0 to 0.875 and 0.0 to 0.763, respectively. • CONCLUSIONS: The developed markers will be valuable for studies of the population structure, mating system, and selection of provenances for restoration projects involving B. integrifolia.

Nonneutral evolution of organelle genes in Silene vulgaris.

Knowledge of mitochondrial gene evolution in angiosperms has taken a dramatic shift within the past decade, from universal slow rates of nucleotide change to a growing realization of high variation in rates among lineages. Additionally, evidence of paternal inheritance of plant mitochondria and recombination among mitochondrial genomes within heteroplasmic individuals has led to speculation about the potential for independent evolution of organellar genes. We report intraspecific mitochondrial and chloroplast sequence variation in a cosmopolitan sample of 42 Silene vulgaris individuals. There was remarkably high variation in two mitochondrial genes (atp1, atp9) and additional variation within a third gene (cob). Tests for patterns of nonneutral evolution were significant for atp1 and atp9, indicative of the maintenance of balanced polymorphisms. Two chloroplast genes (matK, ndhF) possessed less, but still high, variation and no divergence from neutral expectations. Phylogenetic patterns of organelle genes in both the chloroplast and mitochondria were incongruent, indicating the potential for independent evolutionary trajectories. Evidence indicated reassociation among cytoplasmic genomes and recombination between mitochondrial genes and within atp1, implying transient heteroplasmy in ancestral lineages. Although the mechanisms for long-term maintenance of mitochondrial polymorphism are currently unknown, frequency-dependent selection on linked cytoplasmic male sterility genes is a potential candidate.

Reproductive strategy and population variability in the facultative apomict Hieracium pilosella (Asteraceae).

Molecular studies of apomictic plant species often detect more genotypic variation than predicted from their assumed reproductive mode. The two most commonly invoked mechanisms to explain these high levels of variation are recombination, via facultative sexuality, and mutation. The potential for sexual reproduction in the facultative apomict Hieracium pilosella (Asteraceae) was determined at three field sites by artificially pollinating with the closely related, but morphologically distinct, H. aurantiacum. The level of genotypic variation at the three sites was recorded using inter-simple sequence repeats (ISSRs). There was a significant, positive relationship between the measured potential for sexual reproduction and population genotypic variability, indicating that sex has played a role in the structuring of these populations; however, a causal relationship cannot be stated because of the use of regression. We also applied the recently developed method of compatibility analysis. Compatibility analysis can determine, using the occurrence of "character incompatibilities," whether patterns of variation observed in populations are most parsimoniously explained by mutation or recombination. Compatibility analysis also indicated that sexual reproduction had played a role in generating genotypic diversity in these populations. Combining these different types of data may give a greater understanding of the potential for the generation of genotypic diversity in facultative apomictic populations.

The influence of antibiotic production and pre-emptive colonization on the population dynamics of Pantoea agglomerans (Erwinia herbicola) Eh1087 and Erwinia amylovora in planta.

Stigma colonization by Erwinia amylovora is the crucial first step in the development of most fire blight infections in apple and pear trees. Suppression at this point of the disease process by antagonists of E. amylovora, such as Pantoea agglomerans (Erwinia herbicola) strain Eh1087, is a rational approach to control fire blight. We tested the hypothesis that the ability of E. amylovora to compete with Eh1087 for colonization of a stigma is reduced by the potential for Eh1087 to produce the phenazine antibiotic, d-alanylgriseoluteic acid (AGA). In competition experiments on the stigmas of apple flowers, E. amylovora was significantly less successful against Eh1087 (AGA+) than against EhDeltaAGA (AGA-). Further experiments to test the importance of pre-emptive colonization of the stigma by either the pathogen or the antagonist suggested that AGA production significantly enhanced the competitiveness of Eh1087 when it was applied at the same time or 24 h before the pathogen. We also found that pre-emptive stigma colonization by either the pathogen or the antagonist resulted in a population that was resilient to subsequent invasion by a second species suggesting that niche exclusion has a dominant influence on the dynamics of bacterial populations on stigmas.