Insights into the effect of human civilization on Malus evolution and domestication.

The evolutionary history of the Malus genus has not been well studied. In the current study, we presented genetic evidence on the origin of the Malus genus based on genome sequencing of 297 Malus accessions, revealing the genetic relationship between wild species and cultivated apples. Our results demonstrated that North American and East Asian wild species are closer to the outgroup (pear) than Central Asian species, and hybrid species including natural (separated before the Pleistocene, about 2.5 Mya) and artificial hybrids (including ornamental trees and rootstocks) are between East and Central Asian wild species. Introgressions from M. sylvestris in cultivated apples appeared to be more extensive than those from M. sieversii, whose genetic background flowed westward across Eurasia and eastward to wild species including M. prunifolia, M. × asiatica, M. × micromalus, and M. × robust. Our results suggested that the loss of ancestral gene flow from M. sieversii in cultivated apples accompanied the movement of European traders around the world since the Age of Discovery. Natural SNP variations showed that cultivated apples had higher nucleotide diversity than wild species and more unique SNPs than other apple groups. An apple ERECTA-like gene that underwent selection during domestication on 15th chromosome was identified as a likely major determinant of fruit length and diameter, and an NB-ARC domain-containing gene was found to strongly affect anthocyanin accumulation using a genome-wide association approach. Our results provide new insights into the origin and domestication of apples and will be useful in new breeding programmes and efforts to increase fruit crop productivity.

An alpine grasshopper radiation older than the mountains, on Ka Tiritiri o te Moana (Southern Alps) of Aotearoa (New Zealand).

In New Zealand, 13 flightless species of endemic grasshopper are associated with alpine habitats and freeze tolerance. We examined the phylogenetic relationships of the New Zealand species and a subset of Australian alpine grasshoppers using DNA sequences from the entire mitochondrial genome, nuclear 45S rRNA and Histone H3 and H4 loci. Within our sampling, the New Zealand alpine taxa are monophyletic and sister to a pair of alpine Tasmanian grasshoppers. We used six Orthopteran fossils to calibrate a molecular clock analysis to infer that the most recent common ancestor of New Zealand and Tasmanian grasshoppers existed about 20 million years ago, before alpine habitat was available in New Zealand. We inferred a radiation of New Zealand grasshoppers ~13-15 Mya, suggesting alpine species diversification occurred in New Zealand well before the Southern Alps were formed by the mountain building events of the Kaikoura Orogeny 2-5 Mya. This would suggest that either the ancestors of today's New Zealand grasshoppers were not dependent on living in the alpine zone, or they diversified outside of New Zealand.

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.

Ngaokeoke Aotearoa: The Peripatoides Onychophora of New Zealand.

(1) Background: Originally described as a single taxon, Peripatoides novaezealandiae (Hutton, 1876) are distributed across both main islands of New Zealand; the existence of multiple distinct lineages of live-bearing Onychophora across this spatial range has gradually emerged. Morphological conservatism obscured the true endemic diversity, and the inclusion of molecular tools has been instrumental in revealing these cryptic taxa. (2) Methods: Here, we review the diversity of the ovoviviparous Onychophora of New Zealand through a re-analysis of allozyme genotype data, mitochondrial DNA cytochrome oxidase subunit I sequences, geographic information and morphology. (3) Results: New analysis of the multilocus biallelic nuclear data using methods that do not require a priori assumptions of population assignment support at least six lineages of ovoviviparous Peripatoides in northern New Zealand, and mtDNA sequence variation is consistent with these divisions. Expansion of mitochondrial DNA sequence data, including representation of all existing taxa and additional populations extends our knowledge of the scale of sympatry among taxa and shows that three other lineages from southern South Island can be added to the Peripatoides list, and names are proposed here. In total, 10 species of Peripatoides can be recognised with current data.

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.

Mwhitiwhiti Aotearoa: Phylogeny and synonymy of the silent alpine grasshopper radiation of New Zealand (Orthoptera: Acrididae).

Aotearoa New Zealand has a fauna of endemic alpine grasshoppers, consisting of thirteen species distributed among four genera. The many re-classifications of species within this group and the presence of species complexes highlight the uncertainty that surrounds relationships within and between these genera. High-throughput Next Generation Sequencing was used to assemble the complete mitochondrial genomes, 45S ribosomal cassettes and histone sequences of New Zealands four endemic alpine genera: Alpinacris, Brachaspis, Paprides and Sigaus. Phylogenetic analysis of these molecular datasets, as individual genes, partitions and combinations returned a consistent topology that is incompatible with the current classification. The genera Sigaus, Alpinacris, and Paprides all exhibit paraphyly. A consideration of the pronotum, epiphallus and terminalia of adult specimens reveals species-specific differences, but fails to provide compelling evidence for species groups justifying distinct genera. In combination with phylogenetic, morphological and spatial evidence we propose a simplified taxonomy consisting of a single genus for the mwhitiwhiti Aotearoa species radiation.

A multiplexed plant-animal SNP array for selective breeding and species conservation applications.

Reliable and high-throughput genotyping platforms are of immense importance for identifying and dissecting genomic regions controlling important phenotypes, supporting selection processes in breeding programs, and managing wild populations and germplasm collections. Amongst available genotyping tools, single nucleotide polymorphism arrays have been shown to be comparatively easy to use and generate highly accurate genotypic data. Single-species arrays are the most commonly used type so far; however, some multi-species arrays have been developed for closely related species that share single nucleotide polymorphism markers, exploiting inter-species cross-amplification. In this study, the suitability of a multiplexed plant-animal single nucleotide polymorphism array, including both closely and distantly related species, was explored. The performance of the single nucleotide polymorphism array across species for diverse applications, ranging from intra-species diversity assessments to parentage analysis, was assessed. Moreover, the value of genotyping pooled DNA of distantly related species on the single nucleotide polymorphism array as a technique to further reduce costs was evaluated. Single nucleotide polymorphism performance was generally high, and species-specific single nucleotide polymorphisms proved suitable for diverse applications. The multi-species single nucleotide polymorphism array approach reported here could be transferred to other species to achieve cost savings resulting from the increased throughput when several projects use the same array, and the pooling technique adds another highly promising advancement to additionally decrease genotyping costs by half.

Species-wide genomics of kakapo provides tools to accelerate recovery.

The kakapo is a critically endangered, intensively managed, long-lived nocturnal parrot endemic to Aotearoa New Zealand. We generated and analysed 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 kakapo, including prioritizing individuals for translocation and monitoring individuals with poor growth or high disease risk. Overall, by explicitly addressing the challenge of the small sample size, we provide a template for the inclusion of genomic data that will be transformational for species recovery efforts around the globe.

Climate change and alpine-adapted insects: modelling environmental envelopes of a grasshopper radiation.

Mountains create steep environmental gradients that are sensitive barometers of climate change. We calibrated 10 statistical models to formulate ensemble ecological niche models for 12 predominantly alpine, flightless grasshopper species in Aotearoa New Zealand, using their current distributions and current conditions. Niche models were then projected for two future global climate scenarios: representative concentration pathway (RCP) 2.6 (1.0°C rise) and RCP8.5 (3.7°C rise). Results were species specific, with two-thirds of our models suggesting a reduction in potential range for nine species by 2070, but surprisingly, for six species, we predict an increase in potential suitable habitat under mild (+1.0°C) or severe global warming (+3.7°C). However, when the limited dispersal ability of these flightless grasshoppers is taken into account, all 12 species studied are predicted to suffer extreme reductions in range, with a quarter likely to go extinct due to a 96-100% reduction in suitable habitat. Habitat loss is associated with habitat fragmentation that is likely to escalate stochastic vulnerability of remaining populations. Here, we present the predicted outcomes for an endemic radiation of alpine taxa as an exemplar of the challenges that alpine species, both in New Zealand and internationally, are subject to by anthropogenic climate change.

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.

Genome-wide analysis reveals the genetic stock structure of hoki (Macruronus novaezelandiae).

The assessment of the genetic structuring of biodiversity is crucial for management and conservation. This is particularly critical for widely distributed and highly mobile deep-water teleosts, such as hoki (Macruronus novaezelandiae). This species is significant to Maori people and supports the largest commercial fishery in New Zealand, but uncertainty about its stock structure presents a challenge for management. Here, we apply a comprehensive genomic analysis to shed light on the demographic structure of this species by (1) assembling the genome, (2) generating a catalogue of genome-wide SNPs to infer the stock structure and (3) identifying regions of the genome under selection. The final genome assembly used short and long reads and is near complete, representing 93.8% of BUSCO genes, and consisting of 566 contigs totalling 501 Mb. Whole-genome re-sequencing of 510 hoki sampled from 14 locations around New Zealand and Australia, at a read depth greater than 10×, produced 227,490 filtered SNPs. Analyses of these SNPs were able to resolve the stock structure of hoki into two genetically and geographically distinct clusters, one including the Australian and the other one all New Zealand locations, indicating genetic exchange between these regions is limited. Location differences within New Zealand samples were much more subtle (global F ST  = 0.0006), and while small and significant differences could be detected, they did not conclusively identify additional substructures. Ten putative adaptive SNPs were detected within the New Zealand samples, but these also did not geographically partition the dataset further. Contemporary and historical N e estimation suggest the current New Zealand population of hoki is large yet declining. Overall, our study provides the first genomic resources for hoki and provides detailed insights into the fine-scale population genetic structure to inform the management of this species.