A phosphatase gene is linked to nectar dihydroxyacetone accumulation in mānuka (Leptospermum scoparium)

Floral nectar composition beyond common sugars shows great diversity but contributing genetic factors are generally unknown. Mānuka (Leptospermum scoparium) is renowned for the antimicrobial compound methylglyoxal in its derived honey, which originates from the precursor, dihydroxyacetone (DHA), accumulating in the nectar. Although this nectar trait is highly variable, genetic contribution to the trait is unclear. Therefore, we investigated key gene(s) and genomic regions underpinning this trait. We used RNAseq analysis to identify nectary-associated genes differentially expressed between high and low nectar DHA genotypes. We also used a mānuka high-density linkage map and quantitative trait loci (QTL) mapping population, supported by an improved genome assembly, to reveal genetic regions associated with nectar DHA content. Expression and QTL analyses both pointed to the involvement of a phosphatase gene, LsSgpp2. The expression pattern of LsSgpp2 correlated with nectar DHA accumulation, and it co-located with a QTL on chromosome 4. The identification of three QTLs, some of the first reported for a plant nectar trait, indicates polygenic control of DHA content. We have established plant genetics as a key influence on DHA accumulation. The data suggest the hypothesis of LsSGPP2 releasing DHA from DHA-phosphate and variability in LsSgpp2 gene expression contributing to the trait variability.

Blueberry and cranberry pangenomes as a resource for future genetic studies and breeding efforts

Domestication of cranberry and blueberry began in the United States in the early 1800s and 1900s, respectively, and in part owing to their flavors and health-promoting benefits are now cultivated and consumed worldwide. The industry continues to face a wide variety of production challenges (e.g. disease pressures), as well as a demand for higher-yielding cultivars with improved fruit quality characteristics. Unfortunately, molecular tools to help guide breeding efforts for these species have been relatively limited compared with those for other high-value crops. Here, we describe the construction and analysis of the first pangenome for both blueberry and cranberry. Our analysis of these pangenomes revealed both crops exhibit great genetic diversity, including the presence-absence variation of 48.4% genes in highbush blueberry and 47.0% genes in cranberry. Auxiliary genes, those not shared by all cultivars, are significantly enriched with molecular functions associated with disease resistance and the biosynthesis of specialized metabolites, including compounds previously associated with improving fruit quality traits. The discovery of thousands of genes, not present in the previous reference genomes for blueberry and cranberry, will serve as the basis of future research and as potential targets for future breeding efforts. The pangenome, as a multiple-sequence alignment, as well as individual annotated genomes, are publicly available for analysis on the Genome Database for Vaccinium-a curated and integrated web-based relational database. Lastly, the core-gene predictions from the pangenomes will serve useful to develop a community genotyping platform to guide future molecular breeding efforts across the family.

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.

The fate of sex chromosomes during the evolution of monoicy from dioicy in liverworts

Liverworts comprise one of six primary land plant lineages, with the predicted origin of extant liverwort diversity dating to the Silurian. The ancestral liverwort has been inferred to have been dioicous (unisexual) with chromosomal sex determination in which the U chromosome of females and the V chromosome of males were dimorphic with an extensive non-recombining region. In liverworts, sex is determined by a U chromosomal "feminizer" gene that promotes female development, and in its absence, male development ensues. Monoicy (bisexuality) has independently evolved multiple times within liverworts. Here, we explore the evolution of monoicy, focusing on the monoicous species Ricciocarpos natans, and propose that the evolution of monoicy in R. natans involved the appearance of an aneuploid spore that possessed both U and V chromosomes. Chromosomal rearrangements involving the U chromosome resulted in distribution of essential U chromosome genes, including the feminizer, to several autosomal locations. By contrast, we infer that the ancestral V chromosome was inherited largely intact, probably because it carries numerous dispersed "motility" genes distributed across the chromosome. The genetic networks for sex differentiation in R. natans appear largely unchanged except that the feminizer is developmentally regulated, allowing for temporally separated differentiation of female and male reproductive organs on a single plant. A survey of other monoicous liverworts suggests that similar genomic rearrangements may have occurred repeatedly in lineages transitioning to monoicy from dioicy. These data provide a foundation for understanding how genetic networks controlling sex determination can be subtly rewired to produce profound changes in sexual systems.

Blueberry and cranberry pangenomes as a resource for future genetic studies and breeding efforts

Domestication of cranberry and blueberry began in the United States in the early 1800s and 1900s, respectively, and in part owing to their flavors and health-promoting benefits are now cultivated and consumed worldwide. The industry continues to face a wide variety of production challenges (e.g. disease pressures) as well as a demand for higher-yielding cultivars with improved fruit quality characteristics. Unfortunately, molecular tools to help guide breeding efforts for these species have been relatively limited compared with those for other high-value crops. Here, we describe the construction and analysis of the first pangenome for both blueberry and cranberry. Our analysis of these pangenomes revealed both crops exhibit great genetic diversity, including the presence-absence variation of 48.4% genes in highbush blueberry and 47.0% genes in cranberry. Auxiliary genes, those not shared by all cultivars, are significantly enriched with molecular functions associated with disease resistance and the biosynthesis of specialized metabolites, including compounds previously associated with improving fruit quality traits. The discovery of thousands of genes, not present in the previous reference genomes for blueberry and cranberry, will serve as the basis of future research and as potential targets for future breeding efforts. The pangenome, as a multiple-sequence alignment, as well as individual annotated genomes, are publicly available for analysis on the Genome Database for Vaccinium - a curated and integrated web-based relational database. Lastly, the core-gene predictions from the pangenomes will serve useful to develop a community genotyping platform to guide future molecular breeding efforts across the family.

Autopolyploid inheritance and a heterozygous reciprocal translocation shape chromosome genetic behavior in tetraploid blueberry (Vaccinium corymbosum)

Understanding chromosome recombination behavior in polyploidy species is key to advancing genetic discoveries. In blueberry, a tetraploid species, the line of evidences about its genetic behavior still remain poorly understood, owing to the inter-specific, and inter-ploidy admixture of its genome and lack of in depth genome-wide inheritance and comparative structural studies. Here we describe a new high-quality, phased, chromosome-scale genome of a diploid blueberry, clone W85. The genome was integrated with cytogenetics and high-density, genetic maps representing six tetraploid blueberry cultivars, harboring different levels of wild genome admixture, to uncover recombination behavior and structural genome divergence across tetraploid and wild diploid species. Analysis of chromosome inheritance and pairing demonstrated that tetraploid blueberry behaves as an autotetraploid with tetrasomic inheritance. Comparative analysis demonstrated the presence of a reciprocal, heterozygous, translocation spanning one homolog of chr-6 and one of chr-10 in the cultivar Draper. The translocation affects pairing and recombination of chromosomes 6 and 10. Besides the translocation detected in Draper, no other structural genomic divergences were detected across tetraploid cultivars and highly inter-crossable wild diploid species. These findings and resources will facilitate new genetic and comparative genomic studies in Vaccinium and the development of genomic assisted selection strategy for this crop.

Evolution and function of red pigmentation in land plants

BACKGROUND

Land plants commonly produce red pigmentation as a response to environmental stressors, both abiotic and biotic. The type of pigment produced varies among different land plant lineages. In the majority of species they are flavonoids, a large branch of the phenylpropanoid pathway. Flavonoids that can confer red colours include 3-hydroxyanthocyanins, 3-deoxyanthocyanins, sphagnorubins and auronidins, which are the predominant red pigments in flowering plants, ferns, mosses and liverworts, respectively. However, some flowering plants have lost the capacity for anthocyanin biosynthesis and produce nitrogen-containing betalain pigments instead. Some terrestrial algal species also produce red pigmentation as an abiotic stress response, and these include both carotenoid and phenolic pigments.

SCOPE

In this review, we examine: which environmental triggers induce red pigmentation in non-reproductive tissues; theories on the functions of stress-induced pigmentation; the evolution of the biosynthetic pathways; and structure-function aspects of different pigment types. We also compare data on stress-induced pigmentation in land plants with those for terrestrial algae, and discuss possible explanations for the lack of red pigmentation in the hornwort lineage of land plants.

CONCLUSIONS

The evidence suggests that pigment biosynthetic pathways have evolved numerous times in land plants to provide compounds that have red colour to screen damaging photosynthetically active radiation but that also have secondary functions that provide specific benefits to the particular land plant lineage.

Molecular characterization of intergeneric hybrids between Malus and Pyrus

Apple (Malus) and pear (Pyrus) are economically important fruit crops well known for their unique textures, flavours, and nutritional qualities. Both genera are characterised by a distinct pattern of secondary metabolites, which directly affect not only resistance to certain diseases, but also have significant impacts on the flavour and nutritional value of the fruit. The identical chromosome numbers, similar genome size, and their recent divergence date, together with DNA markers have shown that apple and pear genomes are highly co-linear. This study utilized comparative genomic approaches, including simple sequence repeats, high resolution single nucleotide polymorphism melting analysis, and single nucleotide polymorphism chip analysis to identify genetic differences among hybrids of Malus and Pyrus, and F2 offspring. This research has demonstrated and validated that these three marker types, along with metabolomics analysis are very powerful tools to detect and confirm hybridity of progeny derived from crosses between apple and pear in both cross directions. Furthermore, this work analysed the genus-specific metabolite patterns and the resistance to fire blight (Erwinia amylovora) in progeny. The findings of this work will enhance and accelerate the breeding of novel tree fruit crops that benefit producers and consumers, by enabling marker assisted selection of desired traits introgressed between pear and apple.

High-density linkage map construction in an autotetraploid blueberry population and detection of quantitative trait loci for anthocyanin content

Highbush blueberry (Vaccinium corymbosum, 2n = 4x = 48) is the most cultivated type of blueberry, both in New Zealand and overseas. Its perceived nutritional value is conferred by phytonutrients, particularly anthocyanins. Identifying the genetic mechanisms that control the biosynthesis of these metabolites would enable faster development of cultivars with improved fruit qualities. Here, we used recently released tools for genetic mapping in autotetraploids to build a high-density linkage map in highbush blueberry and to detect quantitative trait loci (QTLs) for fruit anthocyanin content. Genotyping was performed by target sequencing, with ∼18,000 single nucleotide polymorphism (SNP) markers being mapped into 12 phased linkage groups (LGs). Fruits were harvested when ripe for two seasons and analyzed with high-performance liquid chromatography-mass spectrometry (HPLC-MS): 25 different anthocyanin compounds were identified and quantified. Two major QTLs that were stable across years were discovered, one on LG2 and one on LG4, and the underlying candidate genes were identified. Interestingly, the presence of anthocyanins containing acylated sugars appeared to be under strong genetic control. Information gained in this study will enable the design of molecular markers for marker-assisted selection and will help build a better understanding of the genetic control of anthocyanin biosynthesis in this crop.

Stress, senescence, and specialized metabolites in bryophytes

Life on land exposes plants to varied abiotic and biotic environmental stresses. These environmental drivers contributed to a large expansion of metabolic capabilities during land plant evolution and species diversification. In this review we summarize knowledge on how the specialized metabolite pathways of bryophytes may contribute to stress tolerance capabilities. Bryophytes are the non-tracheophyte land plant group (comprising the hornworts, liverworts, and mosses) and rapidly diversified following the colonization of land. Mosses and liverworts have as wide a distribution as flowering plants with regard to available environments, able to grow in polar regions through to hot desert landscapes. Yet in contrast to flowering plants, for which the biosynthetic pathways, transcriptional regulation, and compound function of stress tolerance-related metabolite pathways have been extensively characterized, it is only recently that similar data have become available for bryophytes. The bryophyte data are compared with those available for angiosperms, including examining how the differing plant forms of bryophytes and angiosperms may influence specialized metabolite diversity and function. The involvement of stress-induced specialized metabolites in senescence and nutrient response pathways is also discussed.

Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.)

Allele-specific expression (ASE) can lead to phenotypic diversity and evolution. However, the mechanisms regulating ASE are not well understood, particularly in woody perennial plants. In this study, we investigated ASE genes in the apple cultivar 'Royal Gala' (RG). A high quality chromosome-level genome was assembled using a homozygous tetra-haploid RG plant, derived from anther cultures. Using RNA-sequencing (RNA-seq) data from RG flower and fruit tissues, we identified 2091 ASE genes. Compared with the haploid genome of 'Golden Delicious' (GD), a parent of RG, we distinguished the genomic sequences between the two alleles of 817 ASE genes, and further identified allele-specific presence of a transposable element (TE) in the upstream region of 354 ASE genes. These included MYB110a that encodes a transcription factor regulating anthocyanin biosynthesis. Interestingly, another ASE gene, MYB10 also showed an allele-specific TE insertion and was identified using genome data of other apple cultivars. The presence of the TE insertion in both MYB genes was positively associated with ASE and anthocyanin accumulation in apple petals through analysis of 231 apple accessions, and thus underpins apple flower colour evolution. Our study demonstrated the importance of TEs in regulating ASE on a genome-wide scale and presents a novel method for rapid identification of ASE genes and their regulatory elements in plants.

There and back again; historical perspective and future directions for Vaccinium breeding and research studies

The genus Vaccinium L. (Ericaceae) contains a wide diversity of culturally and economically important berry crop species. Consumer demand and scientific research in blueberry (Vaccinium spp.) and cranberry (Vaccinium macrocarpon) have increased worldwide over the crops' relatively short domestication history (~100 years). Other species, including bilberry (Vaccinium myrtillus), lingonberry (Vaccinium vitis-idaea), and ohelo berry (Vaccinium reticulatum) are largely still harvested from the wild but with crop improvement efforts underway. Here, we present a review article on these Vaccinium berry crops on topics that span taxonomy to genetics and genomics to breeding. We highlight the accomplishments made thus far for each of these crops, along their journey from the wild, and propose research areas and questions that will require investments by the community over the coming decades to guide future crop improvement efforts. New tools and resources are needed to underpin the development of superior cultivars that are not only more resilient to various environmental stresses and higher yielding, but also produce fruit that continue to meet a variety of consumer preferences, including fruit quality and health related traits.

Genome-wide patterns of genetic diversity, population structure and demographic history in mānuka (Leptospermum scoparium) growing on indigenous Māori land

Leptospermum scoparium J. R. Forst et G. Forst, known as mānuka by Māori, 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 mānuka 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 mānuka 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, mānuka 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 mānuka by landowners, particularly Māori, and the development of a provenance story for the branding of mānuka 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 Māori 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.

The genome of New Zealand trevally (Carangidae: Pseudocaranx georgianus) uncovers a XY sex determination locus

BACKGROUND

The genetic control of sex determination in teleost species is poorly understood. This is partly because of the diversity of mechanisms that determine sex in this large group of vertebrates, including constitutive genes linked to sex chromosomes, polygenic constitutive mechanisms, environmental factors, hermaphroditism, and unisexuality. Here we use a de novo genome assembly of New Zealand silver trevally (Pseudocaranx georgianus) together with sex-specific whole genome sequencing data to detect sexually divergent genomic regions, identify candidate genes and develop molecular makers.

RESULTS

The de novo assembly of an unsexed trevally (Trevally_v1) resulted in a final assembly of 579.4 Mb in length, with a N50 of 25.2 Mb. Of the assembled scaffolds, 24 were of chromosome scale, ranging from 11 to 31 Mb in length. A total of 28,416 genes were annotated after 12.8 % of the assembly was masked with repetitive elements. Whole genome re-sequencing of 13 wild sexed trevally (seven males and six females) identified two sexually divergent regions located on two scaffolds, including a 6 kb region at the proximal end of chromosome 21. Blast analyses revealed similarity between one region and the aromatase genes cyp19 (a1a/b) (E-value < 1.00E-25, identity > 78.8 %). Males contained higher numbers of heterozygous variants in both regions, while females showed regions of very low read-depth, indicative of male-specificity of this genomic region. Molecular markers were developed and subsequently tested on 96 histologically-sexed fish (42 males and 54 females). Three markers amplified in absolute correspondence with sex (positive in males, negative in females).

CONCLUSIONS

The higher number of heterozygous variants in males combined with the absence of these regions in females support a XY sex-determination model, indicating that the trevally_v1 genome assembly was developed from a male specimen. This sex system contrasts with the ZW sex-determination model documented in closely related carangid species. Our results indicate a sex-determining function of a cyp19a1a-like gene, suggesting the molecular pathway of sex determination is somewhat conserved in this family. The genomic resources developed here will facilitate future comparative work, and enable improved insights into the varied sex determination pathways in teleosts. The sex marker developed in this study will be a valuable resource for aquaculture selective breeding programmes, and for determining sex ratios in wild populations.

The Gillenia trifoliata genome reveals dynamics correlated with growth and reproduction in Rosaceae

The Rosaceae family has striking phenotypic diversity and high syntenic conservation. Gillenia trifoliata is sister species to the Maleae tribe of apple and ~1000 other species. Gillenia has many putative ancestral features, such as herb/sub-shrub habit, dry fruit-bearing and nine base chromosomes. This coalescence of ancestral characters in a phylogenetically important species, positions Gillenia as a 'rosetta stone' for translational science within Rosaceae. We present genomic and phenological resources to facilitate the use of Gillenia for this purpose. The Gillenia genome is the first fully annotated chromosome-level assembly with an ancestral genome complement (x = 9), and with it we developed an improved model of the Rosaceae ancestral genome. MADS and NAC gene family analyses revealed genome dynamics correlated with growth and reproduction and we demonstrate how Gillenia can be a negative control for studying fleshy fruit development in Rosaceae.

Plant biology: Environmental extremes induce a jump in peach fitness

A new study reports that adaptation to climate extremes appears to be driven by replication of a class of transposable elements in peaches and related species. Advanced genomic sequencing techniques may reveal similar events in other plants exposed to extreme stress.

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 Māori 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.

A chromosome-scale assembly of the bilberry genome identifies a complex locus controlling berry anthocyanin composition

Bilberry (Vaccinium myrtillus L.) belongs to the Vaccinium genus, which includes blueberries (Vaccinium spp.) and cranberry (V. macrocarpon). Unlike its cultivated relatives, bilberry remains largely undomesticated, with berry harvesting almost entirely from the wild. As such, it represents an ideal target for genomic analysis, providing comparisons with the domesticated Vaccinium species. Bilberry is prized for its taste and health properties and has provided essential nutrition for Northern European indigenous populations. It contains high concentrations of phytonutrients, with perhaps the most important being the purple colored anthocyanins, found in both skin and flesh. Here, we present the first bilberry genome assembly, comprising 12 pseudochromosomes assembled using Oxford Nanopore (ONT) and Hi-C Technologies. The pseudochromosomes represent 96.6% complete BUSCO genes with an assessed LAI score of 16.3, showing a high conservation of synteny against the blueberry genome. Kmer analysis showed an unusual third peak, indicating the sequenced samples may have been from two individuals. The alternate alleles were purged so that the final assembly represents only one haplotype. A total of 36,404 genes were annotated after nearly 48% of the assembly was masked to remove repeats. To illustrate the genome quality, we describe the complex MYBA locus, and identify the key regulating MYB genes that determine anthocyanin production. The new bilberry genome builds on the genomic resources and knowledge of Vaccinium species, to help understand the genetics underpinning some of the quality attributes that breeding programs aspire to improve. The high conservation of synteny between bilberry and blueberry genomes means that comparative genome mapping can be applied to transfer knowledge about marker-trait association between these two species, as the loci involved in key characters are orthologous.

Peridermal fruit skin formation in Actinidia sp. (kiwifruit) is associated with genetic loci controlling russeting and cuticle formation

BACKGROUND

The skin (exocarp) of fleshy fruit is hugely diverse across species. Most fruit types have a live epidermal skin covered by a layer of cuticle made up of cutin while a few create an outermost layer of dead cells (peridermal layer).

RESULTS

In this study we undertook crosses between epidermal and peridermal skinned kiwifruit, and showed that epidermal skin is a semi-dominant trait. Furthermore, backcrossing these epidermal skinned hybrids to a peridermal skinned fruit created a diverse range of phenotypes ranging from epidermal skinned fruit, through fruit with varying degrees of patches of periderm (russeting), to fruit with a complete periderm. Quantitative trait locus (QTL) analysis of this population suggested that periderm formation was associated with four loci. These QTLs were aligned either to ones associated with russet formation on chromosome 19 and 24, or cuticle integrity and coverage located on chromosomes 3, 11 and 24.

CONCLUSION

From the segregation of skin type and QTL analysis, it appears that skin development in kiwifruit is controlled by two competing factors, cuticle strength and propensity to russet. A strong cuticle will inhibit russeting while a strong propensity to russet can create a continuous dead skinned periderm.

Construction of a high density genetic map for hexaploid kiwifruit (Actinidia chinensis var. deliciosa) using genotyping by sequencing

Commercially grown kiwifruit (genus Actinidia) are generally of two sub-species which have a base haploid genome of 29 chromosomes. The yellow-fleshed Actinidia chinensis var. chinensis, is either diploid (2n = 2x = 58) or tetraploid (2n = 4x = 116) and the green-fleshed cultivar A. chinensis var. deliciosa “Hayward,” is hexaploid (2n = 6x = 174). Advances in breeding green kiwifruit could be greatly sped up by the use of molecular resources for more efficient and faster selection, for example using marker-assisted selection (MAS). The key genetic marker that has been implemented for MAS in hexaploid kiwifruit is for gender testing. The limited marker-trait association has been reported for other polyploid kiwifruit for fruit and production traits. We have constructed a high-density linkage map for hexaploid green kiwifruit using genotyping-by-sequence (GBS). The linkage map obtained consists of 3686 and 3940 markers organized in 183 and 176 linkage groups for the female and male parents, respectively. Both parental linkage maps are co-linear with the A. chinensis “Red5” reference genome of kiwifruit. The linkage map was then used for quantitative trait locus (QTL) mapping, and successfully identified QTLs for king flower number, fruit number and weight, dry matter accumulation, and storage firmness. These are the first QTLs to be reported and discovered for complex traits in hexaploid kiwifruit.

Austropuccinia psidii, causing myrtle rust, has a gigabase-sized genome shaped by transposable elements

Austropuccinia psidii, originating in South America, is a globally invasive fungal plant pathogen that causes rust disease on Myrtaceae. Several biotypes are recognized, with the most widely distributed pandemic biotype spreading throughout the Asia-Pacific and Oceania regions over the last decade. Austropuccinia psidii has a broad host range with more than 480 myrtaceous species. Since first detected in Australia in 2010, the pathogen has caused the near extinction of at least three species and negatively affected commercial production of several Myrtaceae. To enable molecular and evolutionary studies into A. psidii pathogenicity, we assembled a highly contiguous genome for the pandemic biotype. With an estimated haploid genome size of just over 1 Gb (gigabases), it is the largest assembled fungal genome to date. The genome has undergone massive expansion via distinct transposable element (TE) bursts. Over 90% of the genome is covered by TEs predominantly belonging to the Gypsy superfamily. These TE bursts have likely been followed by deamination events of methylated cytosines to silence the repetitive elements. This in turn led to the depletion of CpG sites in TEs and a very low overall GC content of 33.8%. Compared to other Pucciniales, the intergenic distances are increased by an order of magnitude indicating a general insertion of TEs between genes. Overall, we show how TEs shaped the genome evolution of A. psidii and provide a greatly needed resource for strategic approaches to combat disease spread.

QTL Mapping for Resistance to Cankers Induced by Pseudomonas syringae pv. actinidiae (Psa) in a Tetraploid Actinidia chinensis Kiwifruit Population

Polyploidy is a key driver of significant evolutionary changes in plant species. The genus Actinidia (kiwifruit) exhibits multiple ploidy levels, which contribute to novel fruit traits, high yields and resistance to the canker-causing dieback disease incited by Pseudomonas syringae pv. actinidiae (Psa) biovar 3. However, the genetic mechanism for resistance to Psa observed in polyploid kiwifruit is not yet known. In this study we performed detailed genetic analysis of a tetraploid Actinidia chinensis var. chinensis population derived from a cross between a female parent that exhibits weak tolerance to Psa and a highly Psa-resistant male parent. We used the capture-sequencing approach across the whole kiwifruit genome and generated the first ultra-dense maps in a tetraploid kiwifruit population. We located quantitative trait loci (QTLs) for Psa resistance on these maps. Our approach to QTL mapping is based on the use of identity-by-descent trait mapping, which allowed us to relate the contribution of specific alleles from their respective homologues in the male and female parent, to the control of Psa resistance in the progeny. We identified genes in the diploid reference genome whose function is suggested to be involved in plant defense, which underly the QTLs, including receptor-like kinases. Our study is the first to cast light on the genetics of a polyploid kiwifruit and suggest a plausible mechanism for Psa resistance in this species.

Biosynthesis of the Dihydrochalcone Sweetener Trilobatin Requires Phloretin Glycosyltransferase2

Epidemics of obesity and type 2 diabetes drive strong consumer interest in plant-based low-calorie sweeteners. Trilobatin is a sweetener found at high concentrations in the leaves of a range of crabapple (Malus) species, but not in domesticated apple (Malus × domestica) leaves, which contain trilobatin's bitter positional isomer phloridzin. Variation in trilobatin content was mapped to the Trilobatin locus on LG 7 in a segregating population developed from a cross between domesticated apples and crabapples. Phloretin glycosyltransferase2 (PGT2) was identified by activity-directed protein purification and differential gene expression analysis in samples high in trilobatin but low in phloridzin. Markers developed for PGT2 cosegregated strictly with the Trilobatin locus. Biochemical analysis showed PGT2 efficiently catalyzed 4'-o-glycosylation of phloretin to trilobatin as well as 3-hydroxyphloretin to sieboldin. Transient expression of double bond reductase, chalcone synthase, and PGT2 genes reconstituted the apple pathway for trilobatin production in Nicotiana benthamiana Transgenic M. × domestica plants overexpressing PGT2 produced high concentrations of trilobatin in young leaves. Transgenic plants were phenotypically normal, and no differences in disease susceptibility were observed compared to wild-type plants grown under simulated field conditions. Sensory analysis indicated that apple leaf teas from PGT2 transgenics were readily discriminated from control leaf teas and were perceived as significantly sweeter. Identification of PGT2 allows marker-aided selection to be developed to breed apples containing trilobatin, and for high amounts of this natural low-calorie sweetener to be produced via biopharming and metabolic engineering in yeast.

CRISPR-Cas9 enrichment and long read sequencing for fine mapping in plants

BACKGROUND

Genomic methods for identifying causative variants for trait loci applicable to a wide range of germplasm are required for plant biologists and breeders to understand the genetic control of trait variation.

RESULTS

We implemented Cas9-targeted sequencing for fine-mapping in apple, a method combining CRISPR-Cas9 targeted cleavage of a region of interest, followed by enrichment and long-read sequencing using the Oxford Nanopore Technology (ONT). We demonstrated the capability of this methodology to specifically cleave and enrich a plant genomic locus spanning 8 kb. The repeated mini-satellite motif located upstream of the Malus × domestica (apple) MYB10 transcription factor gene, causing red fruit colouration when present in a heterozygous state, was our exemplar to demonstrate the efficiency of this method: it contains a genomic region with a long structural variant normally ignored by short-read sequencing technologiesCleavage specificity of the guide RNAs was demonstrated using polymerase chain reaction products, before using them to specify cleavage of high molecular weight apple DNA. An enriched library was subsequently prepared and sequenced using an ONT MinION flow cell (R.9.4.1). Of the 7,056 ONT reads base-called using both Albacore2 (v2.3.4) and Guppy (v3.2.4), with a median length of 9.78 and 9.89 kb, respectively, 85.35 and 91.38%, aligned to the reference apple genome. Of the aligned reads, 2.98 and 3.04% were on-target with read depths of 180 × and 196 × for Albacore2 and Guppy, respectively, and only five genomic loci were off-target with read depth greater than 25 × , which demonstrated the efficiency of the enrichment method and specificity of the CRISPR-Cas9 cleavage.

CONCLUSIONS

We demonstrated that this method can isolate and resolve single-nucleotide and structural variants at the haplotype level in plant genomic regions. The combination of CRISPR-Cas9 target enrichment and ONT sequencing provides a more efficient technology for fine-mapping loci than genome-walking approaches.

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.

Quantitative trait loci influencing pentacyclic triterpene composition in apple fruit peel

The chemical composition of pentacyclic triterpenes was analysed using a ‘Royal Gala’ x ‘Granny Smith’ segregating population in 2013 and 2015, using apple peels extracted from mature fruit at harvest and after 12 weeks of cold storage. In 2013, 20 compound isoforms from nine unique compound classes were measured for both treatments. In 2015, 20 and 17 compound isoforms from eight unique compound classes were measured at harvest and after cold storage, respectively. In total, 68 quantitative trait loci (QTLs) were detected on 13 linkage groups (LG). Thirty two and 36 QTLs were detected for compounds measured at harvest and after cold storage, respectively. The apple chromosomes with the most QTLs were LG3, LG5, LG9 and LG17. The largest effect QTL was for trihydroxy-urs-12-ene-28-oic acid, located on LG5; this was measured in 2015 after storage, and was inherited from the ‘Royal Gala’ parent (24.9% of the phenotypic variation explained).

Mapping of fire blight resistance in Malus ×robusta 5 flowers following artificial inoculation

BACKGROUND

Although the most common path of infection for fire blight, a severe bacterial disease on apple, is via host plant flowers, quantitative trait loci (QTLs) for fire blight resistance to date have exclusively been mapped following shoot inoculation. It is not known whether the same mechanism underlies flower and shoot resistance.

RESULTS

We report the detection of a fire blight resistance QTL following independent artificial inoculation of flowers and shoots on two F1 segregating populations derived from crossing resistant Malus ×robusta 5 (Mr5) with susceptible 'Idared' and 'Royal Gala' in experimental orchards in Germany and New Zealand, respectively. QTL mapping of phenotypic datasets from artificial flower inoculation of the 'Idared' × Mr5 population with Erwinia amylovora over several years, and of the 'Royal Gala' × Mr5 population in a single year, revealed a single major QTL controlling floral fire blight resistance on linkage group 3 (LG3) of Mr5. This QTL corresponds to the QTL on LG3 reported previously for the 'Idared' × Mr5 and an 'M9' × Mr5 population following shoot inoculation in the glasshouse. Interval mapping of phenotypic data from shoot inoculations of subsets from both flower resistance populations re-confirmed that the resistance QTL is in the same position on LG3 of Mr5 as that for flower inoculation. These results provide strong evidence that fire blight resistance in Mr5 is controlled by a major QTL on LG3, independently of the mode of infection, rootstock and environment.

CONCLUSIONS

This study demonstrates for the first time that resistance to fire blight caused by Erwinia amylovora is independent of the mode of inoculation at least in Malus ×robusta 5.

Genetic control of α-farnesene production in apple fruit and its role in fungal pathogenesis

Terpenes are important compounds in plant trophic interactions. A meta-analysis of GC-MS data from a diverse range of apple (Malus × domestica) genotypes revealed that apple fruit produces a range of terpene volatiles, with the predominant terpene being the acyclic branched sesquiterpene (E,E)-α-farnesene. Four quantitative trait loci (QTLs) for α-farnesene production in ripe fruit were identified in a segregating 'Royal Gala' (RG) × 'Granny Smith' (GS) population with one major QTL on linkage group 10 co-locating with the MdAFS1 (α-farnesene synthase-1) gene. Three of the four QTLs were derived from the GS parent, which was consistent with GC-MS analysis of headspace and solvent-extracted terpenes showing that cold-treated GS apples produced higher levels of (E,E)-α-farnesene than RG. Transgenic RG fruit downregulated for MdAFS1 expression produced significantly lower levels of (E,E)-α-farnesene. To evaluate the role of (E,E)-α-farnesene in fungal pathogenesis, MdAFS1 RNA interference transgenic fruit and RG controls were inoculated with three important apple post-harvest pathogens [Colletotrichum acutatum, Penicillium expansum and Neofabraea alba (synonym Phlyctema vagabunda)]. From results obtained over four seasons, we demonstrate that reduced (E,E)-α-farnesene is associated with decreased disease initiation rates of all three pathogens. In each case, the infection rate was significantly reduced 7 days post-inoculation, although the size of successful lesions was comparable with infections on control fruit. These results indicate that (E,E)-α-farnesene production is likely to be an important factor involved in fungal pathogenesis in apple fruit.

Pseudo-chromosome-length genome assembly of a double haploid "Bartlett" pear (Pyrus communis L.)

BACKGROUND

We report an improved assembly and scaffolding of the European pear (Pyrus communis L.) genome (referred to as BartlettDHv2.0), obtained using a combination of Pacific Biosciences RSII long-read sequencing, Bionano optical mapping, chromatin interaction capture (Hi-C), and genetic mapping. The sample selected for sequencing is a double haploid derived from the same "Bartlett" reference pear that was previously sequenced. Sequencing of di-haploid plants makes assembly more tractable in highly heterozygous species such as P. communis.

FINDINGS

A total of 496.9 Mb corresponding to 97% of the estimated genome size were assembled into 494 scaffolds. Hi-C data and a high-density genetic map allowed us to anchor and orient 87% of the sequence on the 17 pear chromosomes. Approximately 50% (247 Mb) of the genome consists of repetitive sequences. Gene annotation confirmed the presence of 37,445 protein-coding genes, which is 13% fewer than previously predicted.

CONCLUSIONS

We showed that the use of a doubled-haploid plant is an effective solution to the problems presented by high levels of heterozygosity and duplication for the generation of high-quality genome assemblies. We present a high-quality chromosome-scale assembly of the European pear Pyrus communis and demostrate its high degree of synteny with the genomes of Malus x Domestica and Pyrus x bretschneideri.

Two Loci, RiAF3 and RiAF4, Contribute to the Annual-Fruiting Trait in Rubus

Most Rubus species have a biennial cycle of flowering and fruiting with an intervening period of winter dormancy, in common with many perennial fruit crops. Annual-fruiting (AF) varieties of raspberry (Rubus idaeus and Rubus occidentalis L.) and blackberry (Rubus subgenus Rubus) are able to flower and fruit in one growing season, without the intervening dormant period normally required in biennial-fruiting (BF) varieties. We used a red raspberry (R. idaeus) population segregating for AF obtained from a cross between NC493 and 'Chilliwack' to identify genetic factors controlling AF. Genotyping by sequencing (GBS) was used to generate saturated linkage maps in both parents. Trait mapping in this population indicated that AF is controlled by two newly identified loci (RiAF3 and RiAF4) located on Rubus linkage groups (LGs) 3 and 4. The location of these loci was analyzed using single-nucleotide polymorphism (SNP) markers on independent red raspberry and blackberry populations segregating for the AF trait. This confirmed that AF in Rubus is regulated by loci on LG 3 and 4, in addition to a previously reported locus on LG 7. Comparative RNAseq analysis at the time of floral bud differentiation in an AF and a BF variety revealed candidate genes potentially regulating the trait.

Development of a highly efficient Axiom™ 70 K SNP array for Pyrus and evaluation for high-density mapping and germplasm characterization

Background

Both a source of diversity and the development of genomic tools, such as reference genomes and molecular markers, are equally important to enable faster progress in plant breeding. Pear (Pyrus spp.) lags far behind other fruit and nut crops in terms of employment of available genetic resources for new cultivar development. To address this gap, we designed a high-density, high-efficiency and robust single nucleotide polymorphism (SNP) array for pear, with the main objectives of conducting genetic diversity and genome-wide association studies.

Results

By applying a two-step design process, which consisted of the construction of a first ‘draft’ array for the screening of a small subset of samples, we were able to identify the most robust and informative SNPs to include in the Applied Biosystems™ Axiom™ Pear 70 K Genotyping Array, currently the densest SNP array for pear. Preliminary evaluation of this 70 K array in 1416 diverse pear accessions from the USDA National Clonal Germplasm Repository (NCGR) in Corvallis, OR identified 66,616 SNPs (93% of all the tiled SNPs) as high quality and polymorphic (PolyHighResolution). We further used the Axiom Pear 70 K Genotyping Array to construct high-density linkage maps in a bi-parental population, and to make a direct comparison with available genotyping-by-sequencing (GBS) data, which suggested that the SNP array is a more robust method of screening for SNPs than restriction enzyme reduced representation sequence-based genotyping.

Conclusions

The Axiom Pear 70 K Genotyping Array, with its high efficiency in a widely diverse panel of Pyrus species and cultivars, represents a valuable resource for a multitude of molecular studies in pear. The characterization of the USDA-NCGR collection with this array will provide important information for pear geneticists and breeders, as well as for the optimization of conservation strategies for Pyrus.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5712-3) contains supplementary material, which is available to authorized users.

Genetic and cytological analyses reveal the recombination landscape of a partially differentiated plant sex chromosome in kiwifruit

BACKGROUND

Angiosperm sex chromosomes, where present, are generally recently evolved. The key step in initiating the development of sex chromosomes from autosomes is the establishment of a sex-determining locus within a region of non-recombination. To better understand early sex chromosome evolution, it is important to determine the process by which recombination is suppressed around the sex determining genes. We have used the dioecious angiosperm kiwifruit Actinidia chinensis var. chinensis, which has an active-Y sex chromosome system, to study recombination rates around the sex locus, to better understand key events in the development of sex chromosomes.

RESULTS

We have confirmed the sex-determining region (SDR) in A. chinensis var. chinensis, using a combination of high density genetic mapping and fluorescent in situ hybridisation (FISH) of Bacterial Artificial Chromosomes (BACs) linked to the sex markers onto pachytene chromosomes. The SDR is a subtelomeric non-recombining region adjacent to the nucleolar organiser region (NOR). A region of restricted recombination of around 6 Mbp in size in both male and female maps spans the SDR and covers around a third of chromosome 25.

CONCLUSIONS

As recombination is suppressed over a similar region between X chromosomes and between and X and Y chromosomes, we propose that recombination is suppressed in this region because of the proximity of the NOR and the centromere, with both the NOR and centromere suppressing recombination, and this predates suppressed recombination due to differences between X and Y chromosomes. Such regions of suppressed recombination in the genome provide an opportunity for the evolution of sex chromosomes, if a sex-determining locus develops there or translocates into this region.

Indigenous Perspectives and Gene Editing in Aotearoa New Zealand

Gene editing is arguably the most significant recent addition to the modern biotechnology toolbox, bringing both profoundly challenging and enabling opportunities. From a technical point of view the specificity and relative simplicity of these new tools has broadened the potential applications. However, from an ethical point of view it has re-ignited the debates generated by earlier forms of genetic modification. In New Zealand gene editing is currently considered genetic modification and is subject to approval processes under the Environmental Protection Authority (EPA). This process requires decision makers to take into account Māori perspectives. This article outlines previously articulated Māori perspectives on genetic modification and considers the continuing influence of those cultural and ethical arguments within the new context of gene editing. It also explores the range of ways cultural values might be used to analyse the risks and benefits of gene editing in the Aotearoa New Zealand context. Methods used to obtain these perspectives consisted of (a) review of relevant literature regarding lessons learned from the responses of Maori to genetic modification, (b) interviews of selected 'key Maori informants' and (c) surveys of self-selected individuals from groups with interests in either genetics or environmental management. The outcomes of this pilot study identified that while Māori informants were not categorically opposed to new and emerging gene editing technologies a priori, they suggest a dynamic approach to regulation is required where specific uses or types of uses are approved on a case by case basis. This study demonstrates how the cultural cues that Māori referenced in the genetic modification debate continue to be relevant in the context of gene editing but that further work is required to characterize the strength of various positions across the broader community.

Apple whole genome sequences: recent advances and new prospects

In 2010, a major scientific milestone was achieved for tree fruit crops: publication of the first draft whole genome sequence (WGS) for apple (Malus domestica). This WGS, v1.0, was valuable as the initial reference for sequence information, fine mapping, gene discovery, variant discovery, and tool development. A new, high quality apple WGS, GDDH13 v1.1, was released in 2017 and now serves as the reference genome for apple. Over the past decade, these apple WGSs have had an enormous impact on our understanding of apple biological functioning, trait physiology and inheritance, leading to practical applications for improving this highly valued crop. Causal gene identities for phenotypes of fundamental and practical interest can today be discovered much more rapidly. Genome-wide polymorphisms at high genetic resolution are screened efficiently over hundreds to thousands of individuals with new insights into genetic relationships and pedigrees. High-density genetic maps are constructed efficiently and quantitative trait loci for valuable traits are readily associated with positional candidate genes and/or converted into diagnostic tests for breeders. We understand the species, geographical, and genomic origins of domesticated apple more precisely, as well as its relationship to wild relatives. The WGS has turbo-charged application of these classical research steps to crop improvement and drives innovative methods to achieve more durable, environmentally sound, productive, and consumer-desirable apple production. This review includes examples of basic and practical breakthroughs and challenges in using the apple WGSs. Recommendations for "what's next" focus on necessary upgrades to the genome sequence data pool, as well as for use of the data, to reach new frontiers in genomics-based scientific understanding of apple.

Validation of SNP markers for fruit quality and disease resistance loci in apple (Malus × domestica Borkh.) using the OpenArray® platform

Genome mapping has promised much to tree fruit breeding during the last 10 years. Nevertheless, one of the greatest challenges remaining to tree fruit geneticists is the translation of trait loci and whole genome sequences into diagnostic genetic markers that are efficient and cost-effective for use by breeders, who must select genetically optimal parents and subsequently select genetically superior individuals among their progeny. To take this translational step, we designed the apple International RosBREED SNP Consortium OpenArray v1.0 (IRSCOA v1.0) assay using a set of 128 apple single nucleotide polymorphisms (SNPs) linked to fruit quality and pest and disease resistance trait loci. The Thermo Fisher Scientific OpenArray® technology enables multiplexed screening of SNP markers using a real-time PCR instrument with fluorescent probe-based Taqman® assays. We validated the apple IRSCOA v1.0 multi-trait assay by screening 240 phenotyped individuals from the Plant & Food Research apple cultivar breeding programme. This set of individuals comprised commercial and heritage cultivars, elite selections, and families segregating for traits of importance to breeders. In total, 33 SNP markers of the IRSCOA v1.0 were validated for use in marker-assisted selection (MAS) for the scab resistances Rvi2/Vh2, Rvi4/Vh4, Rvi6/Vf, fire blight resistance MR5/RLP1, powdery mildew resistance Pl2, fruit firmness, skin colour, flavour intensity, and acidity. The availability of this set of validated trait-associated SNP markers, which can be used individually on multiple genotyping platforms available to various apple breeding programmes or re-designed using the flanking sequences, represents a large translational genetics step from genomics to crop improvement of apple.

Multiple quantitative trait loci contribute to resistance to bacterial canker incited by Pseudomonas syringae pv. actinidiae in kiwifruit (Actinidia chinensis)

Pseudomonas syringae pv. actinidiae (Psa) biovar 3, a virulent, canker-inducing pathogen is an economic threat to the kiwifruit (Actinidia spp.) industry worldwide. The commercially grown diploid (2×) A. chinensis var. chinensis is more susceptible to Psa than tetraploid and hexaploid kiwifruit. However information on the genetic loci modulating Psa resistance in kiwifruit is not available. Here we report mapping of quantitative trait loci (QTLs) regulating resistance to Psa in a diploid kiwifruit population, derived from a cross between an elite Psa-susceptible 'Hort16A' and a resistant male breeding parent P1. Using high-density genetic maps and intensive phenotyping, we identified a single QTL for Psa resistance on Linkage Group (LG) 27 of 'Hort16A' revealing 16-19% phenotypic variance and candidate alleles for susceptibility and resistance at this loci. In addition, six minor QTLs were identified in P1 on distinct LGs, exerting 4-9% variance. Resistance in the F1 population is improved by additive effects from 'Hort16A' and P1 QTLs providing evidence that divergent genetic pathways interact to combat the virulent Psa strain. Two different bioassays further identified new QTLs for tissue-specific responses to Psa. The genetic marker at LG27 QTL was further verified for association with Psa resistance in diploid Actinidia chinensis populations. Transcriptome analysis of Psa-resistant and susceptible genotypes in field revealed hallmarks of basal defense and provided candidate RNA-biomarkers for screening for Psa resistance in greenhouse conditions.

A near complete, chromosome-scale assembly of the black raspberry (Rubus occidentalis) genome

Background

The fragmented nature of most draft plant genomes has hindered downstream gene discovery, trait mapping for breeding, and other functional genomics applications. There is a pressing need to improve or finish draft plant genome assemblies.

Findings

Here, we present a chromosome-scale assembly of the black raspberry genome using single-molecule real-time Pacific Biosciences sequencing and high-throughput chromatin conformation capture (Hi-C) genome scaffolding. The updated V3 assembly has a contig N50 of 5.1 Mb, representing an ∼200-fold improvement over the previous Illumina-based version. Each of the 235 contigs was anchored and oriented into seven chromosomes, correcting several major misassemblies. Black raspberry V3 contains 47 Mb of new sequences including large pericentromeric regions and thousands of previously unannotated protein-coding genes. Among the new genes are hundreds of expanded tandem gene arrays that were collapsed in the Illumina-based assembly. Detailed comparative genomics with the high-quality V4 woodland strawberry genome (Fragaria vesca) revealed near-perfect 1:1 synteny with dramatic divergence in tandem gene array composition. Lineage-specific tandem gene arrays in black raspberry are related to agronomic traits such as disease resistance and secondary metabolite biosynthesis.

Conclusions

The improved resolution of tandem gene arrays highlights the need to reassemble these highly complex and biologically important regions in draft plant genomes. The updated, high-quality black raspberry reference genome will be useful for comparative genomics across the horticulturally important Rosaceae family and enable the development of marker assisted breeding in Rubus.

Diversification and independent domestication of Asian and European pears

BACKGROUND

Pear (Pyrus) is a globally grown fruit, with thousands of cultivars in five domesticated species and dozens of wild species. However, little is known about the evolutionary history of these pear species and what has contributed to the distinct phenotypic traits between Asian pears and European pears.

RESULTS

We report the genome resequencing of 113 pear accessions from worldwide collections, representing both cultivated and wild pear species. Based on 18,302,883 identified SNPs, we conduct phylogenetics, population structure, gene flow, and selective sweep analyses. Furthermore, we propose a model for the divergence, dissemination, and independent domestication of Asian and European pears in which pear, after originating in southwest China and then being disseminated throughout central Asia, has eventually spread to western Asia, and then on to Europe. We find evidence for rapid evolution and balancing selection for S-RNase genes that have contributed to the maintenance of self-incompatibility, thus promoting outcrossing and accounting for pear genome diversity across the Eurasian continent. In addition, separate selective sweep signatures between Asian pears and European pears, combined with co-localized QTLs and differentially expressed genes, underline distinct phenotypic fruit traits, including flesh texture, sugar, acidity, aroma, and stone cells.

CONCLUSIONS

This study provides further clarification of the evolutionary history of pear along with independent domestication of Asian and European pears. Furthermore, it provides substantive and valuable genomic resources that will significantly advance pear improvement and molecular breeding efforts.

A manually annotated Actinidia chinensis var. chinensis (kiwifruit) genome highlights the challenges associated with draft genomes and gene prediction in plants

Background

Most published genome sequences are drafts, and most are dominated by computational gene prediction. Draft genomes typically incorporate considerable sequence data that are not assigned to chromosomes, and predicted genes without quality confidence measures. The current Actinidia chinensis (kiwifruit) ‘Hongyang’ draft genome has 164 Mb of sequences unassigned to pseudo-chromosomes, and omissions have been identified in the gene models.

Results

A second genome of an A. chinensis (genotype Red5) was fully sequenced. This new sequence resulted in a 554.0 Mb assembly with all but 6 Mb assigned to pseudo-chromosomes. Pseudo-chromosomal comparisons showed a considerable number of translocation events have occurred following a whole genome duplication (WGD) event some consistent with centromeric Robertsonian-like translocations. RNA sequencing data from 12 tissues and ab initio analysis informed a genome-wide manual annotation, using the WebApollo tool. In total, 33,044 gene loci represented by 33,123 isoforms were identified, named and tagged for quality of evidential support. Of these 3114 (9.4%) were identical to a protein within ‘Hongyang’ The Kiwifruit Information Resource (KIR v2). Some proportion of the differences will be varietal polymorphisms. However, as most computationally predicted Red5 models required manual re-annotation this proportion is expected to be small. The quality of the new gene models was tested by fully sequencing 550 cloned ‘Hort16A’ cDNAs and comparing with the predicted protein models for Red5 and both the original ‘Hongyang’ assembly and the revised annotation from KIR v2. Only 48.9% and 63.5% of the cDNAs had a match with 90% identity or better to the original and revised ‘Hongyang’ annotation, respectively, compared with 90.9% to the Red5 models.

Conclusions

Our study highlights the need to take a cautious approach to draft genomes and computationally predicted genes. Our use of the manual annotation tool WebApollo facilitated manual checking and correction of gene models enabling improvement of computational prediction. This utility was especially relevant for certain types of gene families such as the EXPANSIN like genes. Finally, this high quality gene set will supply the kiwifruit and general plant community with a new tool for genomics and other comparative analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4656-3) contains supplementary material, which is available to authorized users.

Chromosome-scale scaffolding of the black raspberry (Rubus occidentalis L.) genome based on chromatin interaction data

Black raspberry (Rubus occidentalis L.) is a niche fruit crop valued for its flavor and potential health benefits. The improvement of fruit and cane characteristics via molecular breeding technologies has been hindered by the lack of a high-quality reference genome. The recently released draft genome for black raspberry (ORUS 4115-3) lacks assembly of scaffolds to chromosome scale. We used high-throughput chromatin conformation capture (Hi-C) and Proximity-Guided Assembly (PGA) to cluster and order 9650 out of 11,936 contigs of this draft genome assembly into seven pseudo-chromosomes. The seven pseudo-chromosomes cover ~97.2% of the total contig length (~223.8 Mb). Locating existing genetic markers on the physical map resolved multiple discrepancies in marker order on the genetic map. Centromeric regions were inferred from recombination frequencies of genetic markers, alignment of 303 bp centromeric sequence with the PGA, and heat map showing the physical contact matrix over the entire genome. We demonstrate a high degree of synteny between each of the seven chromosomes of black raspberry and a high-quality reference genome for strawberry (Fragaria vesca L.) assembled using only PacBio long-read sequences. We conclude that PGA is a cost-effective and rapid method of generating chromosome-scale assemblies from Illumina short-read sequencing data.

Integrated high-density consensus genetic map of Pyrus and anchoring of the 'Bartlett' v1.0 (Pyrus communis) genome

Genetic maps are essential tools for pear genetics and genomics research. In this study, we first constructed an integrated simple sequence repeat (SSR) and single nucleotide polymorphism (SNP)-based consensus genetic map for pear based on common SSR markers between nine published maps. A total of 5,085 markers, including 1,232 SSRs and 3,853 SNPs, were localized on a consensus map spanning 3,266.0 cM in total, with an average marker interval of 0.64 cM, which represents the highest density consensus map of pear to date. Using three sets of high-density SNP-based genetic maps with European pear genetic backgrounds, we anchored a total of 291.5 Mb of the ‘Bartlett’ v1.0 (Pyrus communis L.) genome scaffolds into 17 pseudo-chromosomes. This accounted for 50.5% of the genome assembly, which was a great improvement on the 29.7% achieved originally. Intra-genome and inter-genome synteny analyses of the new ‘Bartlett’ v1.1 genome assembly with the Asian pear ‘Dangshansuli’ (Pyrus bretschneideri Rehd.) and apple (Malus × domestica Borkh.) genomes uncovered four new segmental duplication regions. The integrated high-density SSR and SNP-based consensus genetic map provided new insights into the genetic structure patterns of pear and assisted in the genome assembly of ‘Bartlett’ through further exploration of different pear genetic maps.

Alcohol acyl transferase 1 links two distinct volatile pathways that produce esters and phenylpropenes in apple fruit

Fruit accumulate a diverse set of volatiles including esters and phenylpropenes. Volatile esters are synthesised via fatty acid degradation or from amino acid precursors, with the final step being catalysed by alcohol acyl transferases (AATs). Phenylpropenes are produced as a side branch of the general phenylpropanoid pathway. Major quantitative trait loci (QTLs) on apple (Malus × domestica) linkage group (LG)2 for production of the phenylpropene estragole and volatile esters (including 2-methylbutyl acetate and hexyl acetate) both co-located with the MdAAT1 gene. MdAAT1 has previously been shown to be required for volatile ester production in apple (Plant J., 2014, https://doi.org/10.1111/tpj.12518), and here we show it is also required to produce p-hydroxycinnamyl acetates that serve as substrates for a bifunctional chavicol/eugenol synthase (MdoPhR5) in ripe apple fruit. Fruit from transgenic 'Royal Gala' MdAAT1 knockdown lines produced significantly reduced phenylpropene levels, whilst manipulation of the phenylpropanoid pathway using MdCHS (chalcone synthase) knockout and MdMYB10 over-expression lines increased phenylpropene production. Transient expression of MdAAT1, MdoPhR5 and MdoOMT1 (O-methyltransferase) genes reconstituted the apple pathway to estragole production in tobacco. AATs from ripe strawberry (SAAT1) and tomato (SlAAT1) fruit can also utilise p-coumaryl and coniferyl alcohols, indicating that ripening-related AATs are likely to link volatile ester and phenylpropene production in many different fruit.

Analysis of transcripts differentially expressed between fruited and deflowered 'Gala' adult trees: a contribution to biennial bearing understanding in apple

BACKGROUND

The transition from vegetative to floral state in shoot apical meristems (SAM) is a key event in plant development and is of crucial importance for reproductive success. In perennial plants, this event is recurrent during tree life and subject to both within-tree and between-years heterogeneity. In the present study, our goal was to identify candidate processes involved in the repression or induction of flowering in apical buds of adult apple trees.

RESULTS

Genes differentially expressed (GDE) were examined between trees artificially set in either 'ON' or 'OFF' situation, and in which floral induction (FI) was shown to be inhibited or induced in most buds, respectively, using qRT-PCR and microarray analysis. From the period of FI through to flower differentiation, GDE belonged to four main biological processes (i) response to stimuli, including response to oxidative stress; (ii) cellular processes, (iii) cell wall biogenesis, and (iv) metabolic processes including carbohydrate biosynthesis and lipid metabolic process. Several key regulator genes, especially TEMPRANILLO (TEM), FLORAL TRANSITION AT MERISTEM (FTM1) and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) were found differentially expressed. Moreover, homologs of SPL and Leucine-Rich Repeat proteins were present under QTL zones previously detected for biennial bearing.

CONCLUSIONS

This data set suggests that apical buds of 'ON' and 'OFF' trees were in different physiological states, resulting from different metabolic, hormonal and redox status which are likely to contribute to FI control in adult apple trees. Investigations on carbohydrate and hormonal fluxes from sources to SAM and on cell detoxification process are expected to further contribute to the identification of the underlying physiological mechanisms of FI in adult apple trees.

Scab resistance in 'Geneva' apple is conditioned by a resistance gene cluster with complex genetic control

Apple scab, caused by the fungal pathogen Venturia inaequalis, is one of the most severe diseases of apple worldwide. It is the most studied plant-pathogen interaction involving a woody species using modern genetic, genomic, proteomic and bioinformatic approaches in both species. Although 'Geneva' apple was recognized long ago as a potential source of resistance to scab, this resistance has not been characterized previously. Differential interactions between various monoconidial isolates of V. inaequalis and six segregating F1 and F2 populations indicate the presence of at least five loci governing the resistance in 'Geneva'. The 17 chromosomes of apple were screened using genotyping-by-sequencing, as well as single marker mapping, to position loci controlling the V. inaequalis resistance on linkage group 4. Next, we fine mapped a 5-cM region containing five loci conferring both dominant and recessive scab resistance to the distal end of the linkage group. This region corresponds to 2.2 Mbp (from 20.3 to 22.5 Mbp) on the physical map of 'Golden Delicious' containing nine candidate nucleotide-binding site leucine-rich repeat (NBS-LRR) resistance genes. This study increases our understanding of the complex genetic basis of apple scab resistance conferred by 'Geneva', as well as the gene-for-gene (GfG) relationships between the effector genes in the pathogen and resistance genes in the host.

Genome mapping of postzygotic hybrid necrosis in an interspecific pear population

Deleterious epistatic interactions in plant inter- and intraspecific hybrids can cause a phenomenon known as hybrid necrosis, characterized by a typical seedling phenotype whose main distinguishing features are dwarfism, tissue necrosis and in some cases lethality. Identification of the chromosome regions associated with this type of incompatibility is important not only to increase our understanding of the evolutionary diversification that led to speciation but also for breeding purposes. Development of molecular markers linked to the lethal genes will allow breeders to avoid incompatible inbred combinations that could affect the expression of important agronomic tratis co-segregating with these genes. Although hybrid necrosis has been reported in several plant taxa, including Rosaceae species, this phenomenon has not been described previously in pear. In the interspecific pear population resulting from a cross between PEAR3 (Pyrus bretschneideri × Pyrus communis) and 'Moonglow' (P. communis), we observed two types of hybrid necrosis, expressed at different stages of plant development. Using a combination of previously mapped and newly developed genetic markers, we identified three chromosome regions associated with these two types of lethality, which were genetically independent. One type resulted from a negative epistatic interaction between a locus on linkage group 5 (LG5) of PEAR3 and a locus on LG1 of 'Moonglow', while the second type was due to a gene that maps to LG2 of PEAR3 and which either acts alone or more probably interacts with another gene of unknown location inherited from 'Moonglow'.

Staying green postharvest: how three mutations in the Arabidopsis chlorophyll b reductase gene NYC1 delay degreening by distinct mechanisms

Stresses such as energy deprivation, wounding and water-supply disruption often contribute to rapid deterioration of harvested tissues. To uncover the genetic regulation behind such stresses, a simple assessment system was used to detect senescence mutants in conjunction with two rapid mapping techniques to identify the causal mutations. To demonstrate the power of this approach, immature inflorescences of Arabidopsis plants that contained ethyl methanesulfonate-induced lesions were detached and screened for altered timing of dark-induced senescence. Numerous mutant lines displaying accelerated or delayed timing of senescence relative to wild type were discovered. The underlying mutations in three of these were identified using High Resolution Melting analysis to map to a chromosomal arm followed by a whole-genome sequencing-based mapping method, termed 'Needle in the K-Stack', to identify the causal lesions. All three mutations were single base pair changes and occurred in the same gene, NON-YELLOW COLORING1 (NYC1), a chlorophyll b reductase of the short-chain dehydrogenase/reductase (SDR) superfamily. This was consistent with the mutants preferentially retaining chlorophyll b, although substantial amounts of chlorophyll b were still lost. The single base pair mutations disrupted NYC1 function by three distinct mechanisms, one by producing a termination codon, the second by interfering with correct intron splicing and the third by replacing a highly conserved proline with a non-equivalent serine residue. This non-synonymous amino acid change, which occurred in the NADPH binding domain of NYC1, is the first example of such a mutation in an SDR protein inhibiting a physiological response in plants.

The distribution and impact of common copy-number variation in the genome of the domesticated apple, Malus x domestica Borkh

Background

Copy number variation (CNV) is a common feature of eukaryotic genomes, and a growing body of evidence suggests that genes affected by CNV are enriched in processes that are associated with environmental responses. Here we use next generation sequence (NGS) data to detect copy-number variable regions (CNVRs) within the Malus x domestica genome, as well as to examine their distribution and impact.

Methods

CNVRs were detected using NGS data derived from 30 accessions of M. x domestica analyzed using the read-depth method, as implemented in the CNVrd2 software. To improve the reliability of our results, we developed a quality control and analysis procedure that involved checking for organelle DNA, not repeat masking, and the determination of CNVR identity using a permutation testing procedure.

Results

Overall, we identified 876 CNVRs, which spanned 3.5 % of the apple genome. To verify that detected CNVRs were not artifacts, we analyzed the B- allele-frequencies (BAF) within a single nucleotide polymorphism (SNP) array dataset derived from a screening of 185 individual apple accessions and found the CNVRs were enriched for SNPs having aberrant BAFs (P < 1e-13, Fisher’s Exact test). Putative CNVRs overlapped 845 gene models and were enriched for resistance (R) gene models (P < 1e-22, Fisher’s exact test). Of note was a cluster of resistance gene models on chromosome 2 near a region containing multiple major gene loci conferring resistance to apple scab.

Conclusion

We present the first analysis and catalogue of CNVRs in the M. x domestica genome. The enrichment of the CNVRs with R gene models and their overlap with gene loci of agricultural significance draw attention to a form of unexplored genetic variation in apple. This research will underpin further investigation of the role that CNV plays within the apple genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2096-x) contains supplementary material, which is available to authorized users.