The Genome Database for Rosaceae (GDR): year 10 update

Comparative genomics of N-acetyl-5-methoxytryptamine members in four Prunus species with insights into bud dormancy and abiotic stress responses in Prunus avium

KEY MESSAGE

This study provides novel insights into the evolution, diversification, and functions of melatonin biosynthesis genes in Prunus species, highlighting their potential role in regulating bud dormancy and abiotic stresses. The biosynthesis of melatonin (MEL) in plants is primarily governed by enzymatic reactions involving key enzymes such as serotonin N-acetyltransferase (SNAT), tryptamine 5-hydroxylase (T5H), N-acetylserotonin methyltransferase (ASMT) and tryptophan decarboxylase (TDC). In this study, we analyzed Melatonin genes in four Prunus species such as Prunus avium (Pavi), Prunus pusilliflora (Ppus), Prunus serulata (Pser), and Prunus persica (Pper) based on comparative genomics approach. Among the four Prunus species, a total of 29 TDCs, 998 T5Hs, 16 SNATs, and 115 ASMTs within the genome of four Prunus genomes. A thorough investigation of melatonin-related genes was carried out using systematic biological methods and comparative genomics. Through phylogenetic analysis, orthologous clusters, Go enrichment, syntenic relationship, and gene duplication analysis, we discovered both similarities and variations in Melatonin genes among these Prunus species. Additionally, our study revealed the existence of unique subgroup members in the Melatonin genes of these species, which were distinct from those found in Arabidopsis genes. Furthermore, the transcriptomic expression analysis revealed the potential significance of melatonin genes in bud dormancy regulation and abiotic stresses. Our extensive results offer valuable perspectives on the evolutionary patterns, intricate expansion, and functions of PavMEL genes. Given their promising attributes, PavTDCs, PavT5H, PavNAT, and three PavASMT genes warrant in-depth exploration as prime candidates for manipulating dormancy in sweet cherry. This was done to lay the foundation for future explorations into the structural and functional aspects of these factors in Prunus species. This study offers significant insights into the functions of ASMT, SNAT, T5H, and TDC genes and sheds light on their roles in Prunus avium. Moreover, it established a robust foundation for further exploration functional characterization of melatonin genes in fruit species.

Transcriptomic approach to uncover dynamic events in the development of mid-season sunburn in apple fruit

Apples grown in high heat, high light, and low humidity environments are at risk for sun injury disorders like sunburn and associated crop losses. Understanding the physiological and molecular mechanisms underlying sunburn will support improvement of mitigation strategies and breeding for more resilient varieties. Numerous studies have highlighted key biochemical processes involved in sun injury, such as the phenylpropanoid and reactive oxygen species (ROS) pathways, demonstrating both enzyme activities and expression of related genes in response to sunburn conditions. Most previous studies have focused on at-harvest activity of a small number of genes in response to heat stress. Thus, it remains unclear how stress events earlier in the season affect physiology and gene expression. Here, we applied heat stress to mid-season apples in the field and collected tissue along a time course-24, 48, and 72 h following a heat stimulus-to investigate dynamic gene expression changes using a transcriptomic lens. We found a relatively small number of differentially expressed genes (DEGs) and enriched functional terms in response to heat treatments. Only a few of these belonged to pathways previously described to be involved in sunburn, such as the AsA-GSH pathway, while most DEGs had not yet been implicated in sunburn or heat stress in pome fruit.

Genome-wide identification and characterization of SRLK gene family reveal their roles in self-incompatibility of Erigeron breviscapus

Self-incompatibility (SI) is a reproductive protection mechanism that plants acquired during evolution to prevent self-recession. As the female determinant of SI specificity, SRK has been shown to be the only recognized gene on the stigma and plays important roles in SI response. Asteraceae is the largest family of dicotyledonous plants, many of which exhibit self-incompatibility. However, systematic studies on SRK gene family in Asteraceae are still limited due to lack of high-quality genomic data. In this study, we performed the first systematic genome-wide identification of S-locus receptor like kinases (SRLKs) in the self-incompatible Asteraceae species, Erigeron breviscapus, which is also a widely used perennial medicinal plant endemic to China.52 SRLK genes were identified in the E. breviscapus genome. Structural analysis revealed that the EbSRLK proteins in E. breviscapus are conserved. SRLK proteins from E. breviscapus and other SI plants are clustered into 7 clades, and the majority of the EbSRLK proteins are distributed in Clade I. Chromosomal and duplication analyses indicate that 65% of the EbSRLK genes belong to tandem repeats and could be divided into six tandem gene clusters. Gene expression patterns obtained in E. breviscapus multiple-tissue RNA-Seq data revealed differential temporal and spatial features of EbSRLK genes. Among these, two EbSRLK genes having high expression levels in tongue flowers were cloned. Subcellular localization assay demonstrated that both of their fused proteins are localized on the plasma membrane. All these results indicated that EbSRLK genes possibly involved in SI response in E. breviscapus. This comprehensive genome-wide study of the SRLK gene family in E. breviscapus provides valuable information for understanding the mechanism of SSI in Asteraceae.

Genome of tetraploid sour cherry (Prunus cerasus L.) 'Montmorency' identifies three distinct ancestral Prunus genomes

Sour cherry (Prunus cerasus L.) is a valuable fruit crop in the Rosaceae family and a hybrid between progenitors closely related to extant Prunus fruticosa (ground cherry) and Prunus avium (sweet cherry). Here we report a chromosome-scale genome assembly for sour cherry cultivar Montmorency, the predominant cultivar grown in the USA. We also generated a draft assembly of P. fruticosa to use alongside a published P. avium sequence for syntelog-based subgenome assignments for 'Montmorency' and provide compelling evidence P. fruticosa is also an allotetraploid. Using hierarchal k-mer clustering and phylogenomics, we show 'Montmorency' is trigenomic, containing two distinct subgenomes inherited from a P. fruticosa-like ancestor (A and A') and two copies of the same subgenome inherited from a P. avium-like ancestor (BB). The genome composition of 'Montmorency' is AA'BB and little-to-no recombination has occurred between progenitor subgenomes (A/A' and B). In Prunus, two known classes of genes are important to breeding strategies: the self-incompatibility loci (S-alleles), which determine compatible crosses, successful fertilization, and fruit set, and the Dormancy Associated MADS-box genes (DAMs), which strongly affect dormancy transitions and flowering time. The S-alleles and DAMs in 'Montmorency' and P. fruticosa were manually annotated and support subgenome assignments. Lastly, the hybridization event 'Montmorency' is descended from was estimated to have occurred less than 1.61 million years ago, making sour cherry a relatively recent allotetraploid. The 'Montmorency' genome highlights the evolutionary complexity of the genus Prunus and will inform future breeding strategies for sour cherry, comparative genomics in the Rosaceae, and questions regarding neopolyploidy.

Pedigree-based analysis in multi-parental diploid rose populations reveals QTLs for cercospora leaf spot disease resistance

Cercospora leaf spot (CLS) (Cercospora rosicola) is a major fungal disease of roses (Rosa sp.) in the southeastern U.S. Developing CLS-resistant cultivars offers a potential solution to reduce pesticide use. Yet, no work has been performed on CLS resistance. This study aimed to identify QTLs and to characterize alleles for resistance to CLS. The study used pedigree-based QTL analysis to dissect the genetic basis of CLS resistance using two multi-parental diploid rose populations (TX2WOB and TX2WSE) evaluated across five years in two Texas locations. A total 38 QTLs were identified across both populations and distributed over all linkage groups. Three QTLs on LG3, LG4, and LG6 were consistently mapped over multiple environments. The LG3 QTL was mapped in a region between 18.9 and 27.8 Mbp on the Rosa chinensis genome assembly. This QTL explained 13 to 25% of phenotypic variance. The LG4 QTL detected in the TX2WOB population spanned a 35.2 to 39.7 Mbp region with phenotypic variance explained (PVE) up to 48%. The LG6 QTL detected in the TX2WSE population was localized to 17.9 to 33.6 Mbp interval with PVE up to 36%. Also, this study found multiple degrees of favorable allele effects (q-allele) associated with decreasing CLS at major loci. Ancestors 'OB', 'Violette', and PP-M4-4 were sources of resistance q-alleles. These results will aid breeders in parental selection to develop CLS-resistant rose cultivars. Ultimately, high throughput DNA tests that target major loci for CLS could be developed for routine use in a DNA-informed breeding program.

Comparison of structural variants detected by PacBio-CLR and ONT sequencing in pear

BACKGROUND

Structural variations (SVs) have recently become a topic of great interest in the area of genetic diversity and trait regulation. As genomic sequencing technologies have rapidly advanced, longer reads have been used to identify SVs at high resolution and with increased accuracy. It is important to choose a suitable sequencing platform and appropriate sequencing depth for SV detection in the pear genome.

RESULTS

In this study, two types of long reads from sequencing platforms, continuous long reads from Pacific Biosciences (PB-CLR) and long reads from Oxford Nanopore Technologies (ONT), were used to comprehensively analyze and compare SVs in the pear genome. The mapping rate of long reads was higher when the program Minimap2 rather than the other three mapping tools (NGMLR, LRA and Winnowmap2) was used. Three SV detection programs (Sniffles_v2, CuteSV, and Nanovar) were compared, and Nanovar had the highest sensitivity in detecting SVs at low sequencing depth (10-15×). A sequencing depth of 15× was suitable for SV detection in the pear genome using Nanovar. SVs detected by Sniffles_v2 and CuteSV with ONT reads had the high overlap with presence/absence variations (PAVs) in the pear cultivars 'Bartlett' and 'Dangshansuli', both of them with 38% of insertions and 55% of deletions overlapping with PAVs at sequencing depth of 30×. For the ONT sequencing data, over 37,526 SVs spanning ~ 28 Mb were identified by all three software packages for the 'Bartlett' and 'Dangshansuli' genomes. Those SVs were annotated and combined with transcriptome profiles derived from 'Bartlett' and 'Dangshansuli' fruit flesh at 60 days after cross-pollination. Several genes related to levels of sugars, acid, stone cells, and aromatic compounds were identified among the SVs. Transcription factors were then predicted among those genes, and results included bHLH, ERF, and MYB genes.

CONCLUSION

SV detection is of great significance in exploring phenotypic differences between pear varieties. Our study provides a framework for assessment of different SV software packages and sequencing platforms that can be applied in other plant genome studies. Based on these analyses, ONT sequencing data was determined to be more suitable than PB-CLR for SV detection in the pear genome. This analysis model will facilitate screening of genes related to agronomic traits in other crops.

Transcription factor MdLSD1 negatively regulates α-farnesene biosynthesis in apple-fruit skin tissue

α-Farnesene is a sesquiterpene present in plants. It was first discovered in apples. It plays an important role in the plant defence response and is considered a key factor in the occurrence of superficial scald. The gene encoding α-farnesene synthase, which is the last key enzyme in the biosynthetic pathway of α-farnesene in apple fruit, has become the primary target enzyme for controlling the genetic manipulation of α-farnesene biosynthesis. In this study, the yeast one-hybrid assay and the dual luciferase assay were used to ascertain the relationship between MdLSD1 and MdAFS. Real-time PCR was used to analyse the molecular mechanism underlying the regulation of MdAFS by MdLSD1. Our results revealed that transcription factor MdLSD1, which is closely related to programmed cell death in apple fruit tissues, binds to MdAFS. Transient transformation of apple skin with vectors overexpressing MdLSD1 showed that the gene negatively regulates MdAFS. Overall, we suggest that MdLSD1 negatively regulates MdAFS. Our results are of great significance for future research on the transcriptional regulation of the α-farnesene synthase gene and provide a new direction for exploring the specific mechanism of programmed cell death involved in superficial-scald incidence.

The banana genome hub: a community database for genomics in the Musaceae

The Banana Genome Hub provides centralized access for genome assemblies, annotations, and the extensive related omics resources available for bananas and banana relatives. A series of tools and unique interfaces are implemented to harness the potential of genomics in bananas, leveraging the power of comparative analysis, while recognizing the differences between datasets. Besides effective genomic tools like BLAST and the JBrowse genome browser, additional interfaces enable advanced gene search and gene family analyses including multiple alignments and phylogenies. A synteny viewer enables the comparison of genome structures between chromosome-scale assemblies. Interfaces for differential expression analyses, metabolic pathways and GO enrichment were also added. A catalogue of variants spanning the banana diversity is made available for exploration, filtering, and export to a wide variety of software. Furthermore, we implemented new ways to graphically explore gene presence-absence in pangenomes as well as genome ancestry mosaics for cultivated bananas. Besides, to guide the community in future sequencing efforts, we provide recommendations for nomenclature of locus tags and a curated list of public genomic resources (assemblies, resequencing, high density genotyping) and upcoming resources-planned, ongoing or not yet public. The Banana Genome Hub aims at supporting the banana scientific community for basic, translational, and applied research and can be accessed at https://banana-genome-hub.southgreen.fr.

Evolutionary and Integrative Analysis of Gibberellin-Dioxygenase Gene Family and Their Expression Profile in Three Rosaceae Genomes (F. vesca, P. mume, and P. avium) Under Phytohormone Stress

The gibberellin-dioxygenase (GAox) gene family plays a crucial role in regulating plant growth and development. GAoxs, which are encoded by many gene subfamilies, are extremely critical in regulating bioactive GA levels by catalyzing the subsequent stages in the biosynthesis process. Moreover, GAoxs are important enzymes in the GA synthesis pathway, and the GAox gene family has not yet been identified in Rosaceae species (Prunus avium L., F. vesca, and P. mume), especially in response to gibberellin and PCa (prohexadione calcium; reduce biologically active GAs). In the current investigation, 399 GAox members were identified in sweet cherry, Japanese apricot, and strawberry. Moreover, they were further classified into six (A-F) subgroups based on phylogeny. According to motif analysis and gene structure, the majority of the PavGAox genes have a remarkably well-maintained exon-intron and motif arrangement within the same subgroup, which may lead to functional divergence. In the systematic investigation, PavGAox genes have several duplication events, but segmental duplication occurs frequently. A calculative analysis of orthologous gene pairs in Prunus avium L., F. vesca, and P. mume revealed that GAox genes are subjected to purifying selection during the evolutionary process, resulting in functional divergence. The analysis of cis-regulatory elements in the upstream region of the 140 PavGAox members suggests a possible relationship between genes and specific functions of hormone response-related elements. Moreover, the PavGAox genes display a variety of tissue expression patterns in diverse tissues, with most of the PavGAox genes displaying tissue-specific expression patterns. Furthermore, most of the PavGAox genes express significant expression in buds under phytohormonal stresses. Phytohormones stress analysis demonstrated that some of PavGAox genes are responsible for maintaining the GA level in plant-like Pav co4017001.1 g010.1.br, Pav sc0000024.1 g340.1.br, and Pav sc0000024.1 g270.1.mk. The subcellular localization of PavGAox protein utilizing a tobacco transient transformation system into the tobacco epidermal cells predicted that GFP signals were mostly found in the cytoplasm. These findings will contribute to a better understanding of the GAox gene family's interaction with prohexadione calcium and GA, as well as provide a strong framework for future functional characterization of GAox genes in sweet cherry.

Genome-wide analysis of the rose (Rosa chinensis) NAC family and characterization of RcNAC091

A genome-wide analysis identified 116 NAC genes in rose, including stress-related ones with different expression patterns under drought and salt stress. Silencing of RcNAC091, a member of the ATAF subfamily, decreased dehydration tolerance in rose. The NAC (NAM, ATAF, and CUC) transcription factors (TFs) are plant-specific proteins that regulate various developmental processes and stress responses. However, knowledge of the NAC TFs in rose (Rosa chinensis), one of the most important horticultural crops, is limited. In this study, 116 NAC genes were identified from the rose genome and classified into 16 subfamilies based on protein phylogenetic analysis. Chromosomal mapping revealed that the RcNAC genes were unevenly distributed on the seven chromosomes of rose. Gene structure and motif analysis identified a total of ten conserved motifs, of which motifs 1-7 were highly conserved and present in most rose NACs, while motifs 8-10 were present only in a few subfamilies. Further study of the stress-related RcNACs based on the transcriptome data showed differences in the expression patterns among the organs, at various floral developmental stages, and under drought and salt stress in rose leaves and roots. The stress-related RcNACs possessed cis-regulatory elements (CREs) categorized into three groups corresponding to plant growth and development, phytohormone response, and abiotic and biotic stress response. Reverse transcription-quantitative real-time PCR (RT-qPCR) analysis of 11 representative RcNACs revealed their differential expression in rose leaves and roots under abscisic acid (ABA), polyethylene glycol (PEG), and sodium chloride (NaCl) treatments. Furthermore, the silencing of RcNAC091 verified its role in positively regulating the dehydration stress response. Overall, the present study provides valuable insights into stress-related RcNACs and paves the way for stress tolerance in rose.

MYB transcription factor family in sweet cherry (Prunus avium L.): genome-wide investigation, evolution, structure, characterization and expression patterns

BACK GROUND

MYB Transcription factors (TFs) are most imperative and largest gene family in plants, which participate in development, metabolism, defense, differentiation and stress response. The MYB TFs has been studied in various plant species. However, comprehensive studies of MYB gene family in the sweet cherry (Prunus avium L.) are still unknown.

RESULTS

In the current study, a total of 69 MYB genes were investigated from sweet cherry genome and classified into 28 subfamilies (C1-C28 based on phylogenetic and structural analysis). Microcollinearity analysis revealed that dispersed duplication (DSD) events might play an important role in the MYB genes family expansion. Chromosomal localization, the synonymous (Ks) and nonsynonymous (Ka) analysis, molecular characteristics (pI, weight and length of amino acids) and subcellular localization were accomplished using several bioinformatics tools. Furthermore, the members of distinct subfamilies have diverse cis-acting regions, conserved motifs, and intron-exon architectures, indicating functional heterogeneity in the MYB family. Moreover, the transcriptomic data exposed that MYB genes might play vital role in bud dormancy. The quantitative real-time qRT-PCR was carried out and the expression pattern indicated that MYB genes significantly expressed in floral bud as compared to flower and fruit.

CONCLUSION

Our comprehensive findings provide supportive insights into the evolutions, expansion complexity and functionality of PavMYB genes. These PavMYB genes should be further investigated as they seem to be brilliant candidates for dormancy manipulation in sweet cherry.

Identification and Comprehensive Genome-Wide Analysis of Glutathione S-Transferase Gene Family in Sweet Cherry (Prunus avium) and Their Expression Profiling Reveals a Likely Role in Anthocyanin Accumulation

Glutathione S-transferases (GSTs) in plants are multipurpose enzymes that are involved in growth and development and anthocyanins transportation. However, members of the GST gene family were not identified in sweet cherry (Prunus avium). To identify the GST genes in sweet cherry, a genome-wide analysis was conducted. In this study, we identified 67 GST genes in P. avium genome and nomenclature according to chromosomal distribution. Phylogenetic tree analysis revealed that PavGST genes were classified into seven chief subfamily: TCHQD, Theta, Phi, Zeta, Lambda, DHAR, and Tau. The majority of the PavGST genes had a relatively well-maintained exon-intron and motif arrangement within the same group, according to gene structure and motif analyses. Gene structure (introns-exons) and conserved motif analysis revealed that the majority of the PavGST genes showed a relatively well-maintained motif and exons-introns configuration within the same group. The chromosomal localization, GO enrichment annotation, subcellular localization, syntenic relationship, Ka/Ks analysis, and molecular characteristics were accomplished using various bioinformatics tools. Mode of gene duplication showed that dispersed duplication might play a key role in the expansion of PavGST gene family. Promoter regions of PavGST genes contain numerous cis-regulatory components, which are involved in multiple stress responses, such as abiotic stress and phytohormones responsive factors. Furthermore, the expression profile of sweet cherry PavGSTs showed significant results under LED treatment. Our findings provide the groundwork for future research into induced LED anthocyanin and antioxidants deposition in sweet cherries.

Peduncle Necking in Rosa hybrida Induces Stress-Related Transcription Factors, Upregulates Galactose Metabolism, and Downregulates Phenylpropanoid Biosynthesis Genes

Roses are highly valued as cut flowers worldwide but have limited vase life. Peduncle bending "bent neck" or "necking" is a major cause of reduced vase life, especially in some cultivars. Necking is thought to be caused by either an air embolism or accumulation of microorganisms at or within the stem end, blocking the xylem vessels and preventing water uptake. However, the underlying mechanisms of necking are poorly understood. Here, RNAseq analysis was applied to compare gene expression across three stages of peduncle necking (straight, <90°, and >90°), in the necking-susceptible Rosa hybrida cultivar H30. Most gene expression change was later in bending and there was, overall, more downregulation than upregulation of gene expression during necking. Photosynthetic, starch, and lignin biosynthesis genes were all downregulated, while genes associated with galactose metabolism, producing raffinose and trehalose that are both related to osmoprotection, were upregulated. Genes associated with starch breakdown, autophagy, and senescence were also upregulated, as were most of the NAC and WRKY transcription factors, involved in stress and senescence regulation. Microscopy showed a cellular collapse in the peduncle. These data support a possible mechanism, whereby a reduction in water transport leads to a cellular collapse in the peduncle, accompanied by upregulation of senescence and drought responses.

qtlXplorer: an online systems genetics browser in the Eucalyptus Genome Integrative Explorer (EucGenIE)

BACKGROUND

Affordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking.

RESULTS

We have developed qtlXplorer ( https://eucgenie.org/QTLXplorer ) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases.

CONCLUSIONS

qtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data.

Genome-wide investigation and comparative analysis of MATE gene family in Rosaceae species and their regulatory role in abiotic stress responses in Chinese pear (Pyrus bretschneideri)

The Multidrug and Toxic Compound Extrusion (MATE) protein belongs to a secondary transporter gene family, which plays a primary role in transporting many kinds of substrates such as organic compounds, secondary metabolites, and phytohormones. MATE protein members exist in both prokaryotes and eukaryotes. However, evolution and comprehensive analysis of the MATE genes has not been performed in Rosaceae species. In the present study, a total of 404 MATEs genes were identified from six Rosaceae genomes (Prunus avium, Pyrus bretschneideri, Prunus persica, Fragaria vesca, Prunus mume, and Malus domestica) and classified into eight main subfamilies (I-VII) based on structural and phylogenetic analysis. Microcollinearity analysis showed that whole-genome duplication events might play a vital role in the expansion of the MATE genes family. The Ka/Ks analysis, chromosomal localization, subcellular localization, and molecular characteristics (length, weight, and pI) were performed using various bioinformatics tools. Furthermore, different subfamilies have different introns-exons structures, cis-acting elements, and conserved motifs analysis, indicating functional divergence in the MATE family. Subsequently, RNA-seq analysis and real-time qRT-PCR were conducted during Chinese pear fruit development. Moreover, PbMATE genes were significantly expressed under hormonal treatments of MeJA (methyl jasmonate), SA (salicylic acid), and ABA (abscisic acid). Overall, our results provide helpful insights into the functions, expansion complexity, and evolutions of the MATE genes in Chinese pear and five Rosaceae species.

Tripal and Galaxy: supporting reproducible scientific workflows for community biological databases

Online biological databases housing genomics, genetic and breeding data can be constructed using the Tripal toolkit. Tripal is an open-source, internationally developed framework that implements FAIR data principles and is meant to ease the burden of constructing such websites for research communities. Use of a common, open framework improves the sustainability and manageability of such as site. Site developers can create extensions for their site and in turn share those extensions with others. One challenge that community databases often face is the need to provide tools for their users that analyze increasingly larger datasets using multiple software tools strung together in a scientific workflow on complicated computational resources. The Tripal Galaxy module, a ‘plug-in’ for Tripal, meets this need through integration of Tripal with the Galaxy Project workflow management system. Site developers can create workflows appropriate to the needs of their community using Galaxy and then share those for execution on their Tripal sites via automatically constructed, but configurable, web forms or using an application programming interface to power web-based analytical applications. The Tripal Galaxy module helps reduce duplication of effort by allowing site developers to spend time constructing workflows and building their applications rather than rebuilding infrastructure for job management of multi-step applications.

Genome-wide association mapping identifies novel loci underlying fire blight resistance in apple

Fire blight, caused by epiphytotic gram-negative bacteria Erwinia amylovora, is the most destructive bacterial disease of apple (Malus spp.). Genetic mechanisms of fire blight resistance have mainly been studied using traditional biparental quantitative trait loci (QTL) mapping approaches. Here, we use large-scale historic shoot and blossom fire blight data collected over multiple years and genotyping-by-sequencing (GBS) markers to identify significant marker-trait associations in a diverse set of 566 apple [Malus domestica (Suckow) Borkh.] accessions. There was large variation in fire blight resistance and susceptibility in these accessions. We identified 23 and 38 QTL significantly (p < .001) associated with shoot and blossom blight resistance, respectively. The QTL are distributed across all 17 chromosomes of apple. Four shoot blight and 19 blossom blight QTL identified in this study colocalized with previously identified QTL associated with resistance to fire blight or apple scab. Using transcriptomics data of two apple cultivars with contrasting fire blight responses, we also identified candidate genes for fire blight resistance that are differentially expressed between resistant and susceptible cultivars and located within QTL intervals for fire blight resistance. However, further experiments are needed to confirm and validate these marker-trait associations and develop diagnostic markers before use in marker-assisted breeding to develop apple cultivars with decreased fire blight susceptibility.

Mapping and characterization QTLs for phenological traits in seven pedigree-connected peach families

Background

Environmental adaptation and expanding harvest seasons are primary goals of most peach [Prunus persica (L.) Batsch] breeding programs. Breeding perennial crops is a challenging task due to their long breeding cycles and large tree size. Pedigree-based analysis using pedigreed families followed by haplotype construction creates a platform for QTL and marker identification, validation, and the use of marker-assisted selection in breeding programs.

Results

Phenotypic data of seven F₁ low to medium chill full-sib families were collected over 2 years at two locations and genotyped using the 9 K SNP Illumina array. Three QTLs were discovered for bloom date (BD) and mapped on linkage group 1 (LG1) (172–182 cM), LG4 (48–54 cM), and LG7 (62–70 cM), explaining 17–54%, 11–55%, and 11–18% of the phenotypic variance, respectively. The QTL for ripening date (RD) and fruit development period (FDP) on LG4 was co-localized at the central part of LG4 (40–46 cM) and explained between 40 and 75% of the phenotypic variance. Haplotype analyses revealed SNP haplotypes and predictive SNP marker(s) associated with desired QTL alleles and the presence of multiple functional alleles with different effects for a single locus for RD and FDP.

Conclusions

A multiple pedigree-linked families approach validated major QTLs for the three key phenological traits which were reported in previous studies across diverse materials, geographical distributions, and QTL mapping methods. Haplotype characterization of these genomic regions differentiates this study from the previous QTL studies. Our results will provide the peach breeder with the haplotypes for three BD QTLs and one RD/FDP QTL to create predictive DNA-based molecular marker tests to select parents and/or seedlings that have desired QTL alleles and cull unwanted genotypes in early seedling stages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07483-8.

The NanDeSyn database for Nannochloropsis systems and synthetic biology

Nannochloropsis species, unicellular industrial oleaginous microalgae, are model organisms for microalgal systems and synthetic biology. To facilitate community-based annotation and mining of the rapidly accumulating functional genomics resources, we have initiated an international consortium and present a comprehensive multi-omics resource database named Nannochloropsis Design and Synthesis (NanDeSyn; http://nandesyn.single-cell.cn). Via the Tripal toolkit, it features user-friendly interfaces hosting genomic resources with gene annotations and transcriptomic and proteomic data for six Nannochloropsis species, including two updated genomes of Nannochloropsis oceanica IMET1 and Nannochloropsis salina CCMP1776. Toolboxes for search, Blast, synteny view, enrichment analysis, metabolic pathway analysis, a genome browser, etc. are also included. In addition, functional validation of genes is indicated based on phenotypes of mutants and relevant bibliography. Furthermore, epigenomic resources are also incorporated, especially for sequencing of small RNAs including microRNAs and circular RNAs. Such comprehensive and integrated landscapes of Nannochloropsis genomics and epigenomics will promote and accelerate community efforts in systems and synthetic biology of these industrially important microalgae.

Genetic determinism of prickles in rose

KEY MESSAGE

The genetic determinism of prickle in rose is complex, with a major locus on LG3 that controls the absence/presence of prickles on the rose stem. Rose is one of the major ornamental plants. The selection of glabrous cultivars is an important breeding target but remains a difficult task due to our limited genetic knowledge. Our objective was to understand the genetic and molecular determinism of prickles. Using a segregating diploid rose F1 population, we detected two types of prickles (glandular and non-glandular) in the progeny. We scored the number of non-glandular prickles on the floral and main stems for three years. We performed QTL analysis and detected four prickle loci on LG1, 3, 4 and 6. We determined the credible interval on the reference genome. The QTL on LG3 is a major locus that controls the presence of prickles, and three QTLs (LG3, 4 and 1) may be responsible for prickle density. We further revealed that glabrous hybrids are caused by the combination of the two recessive alleles from both parents. In order to test whether rose prickles could originate from a 'trichome-like structure,' we used a candidate approach to characterize rose gene homologues known in Arabidopsis, involved in trichome initiation. Four of these homologues were located within the overlapping credible interval of the detected QTLs. Transcript accumulation analysis weakly supports the involvement of trichome homologous genes, in the molecular control of prickle initiation. Our studies provide strong evidence for a complex genetic determinism of stem prickle and could help to establish guidelines for glabrous rose breeding. New insights into the relationship between prickles and trichomes constitute valuable information for reverse genetic research on prickles.

Strawberry FaNAC2 Enhances Tolerance to Abiotic Stress by Regulating Proline Metabolism

The quality and yields of strawberry plants are seriously affected by abiotic stress every year. NAC (NAM, ATAF, CUC) transcription factors are plant-specific, having various functions in plant development and response to stress. In our study, FaNAC2 from strawberry (Fragaria × ananassa, cultivar "Benihoppe") was isolated and found to be a member of the ATAF sub-family, belonging to the NAC family of transcription factors. FaNAC2 was strongly expressed in the shoot apical meristem and older leaves of strawberries, and was induced by cold, high salinity, and drought stress. To investigate how FaNAC2 functions in plant responses to abiotic stress, transgenic Nicotiana benthamiana plants ectopically overexpressing FaNAC2 were generated. The transgenic plants grew better under salt and cold stress, and, during simulated drought treatment, these transgenic lines not only grew better, but also showed higher seed germination rates than wild-type plants. Gene expression analysis revealed that key genes in proline biosynthesis pathways were up-regulated in FaNAC2 overexpression lines, while its catabolic pathway genes were down-regulated and proline was accumulated more with the overexpression of FaNAC2 after stress treatments. Furthermore, the gene expression of abscisic acid biosynthesis was also promoted. Our results demonstrate that FaNAC2 plays an important positive role in response to different abiotic stresses and may be further utilized to improve the stress tolerance of strawberry plants.

Identification and characterization of QTLs for fruit quality traits in peach through a multi-family approach

BACKGROUND

Fruit quality traits have a significant effect on consumer acceptance and subsequently on peach (Prunus persica (L.) Batsch) consumption. Determining the genetic bases of key fruit quality traits is essential for the industry to improve fruit quality and increase consumption. Pedigree-based analysis across multiple peach pedigrees can identify the genomic basis of complex traits for direct implementation in marker-assisted selection. This strategy provides breeders with better-informed decisions and improves selection efficiency and, subsequently, saves resources and time.

RESULTS

Phenotypic data of seven F1 low to medium chill full-sib families were collected over 2 years at two locations and genotyped using the 9 K SNP Illumina array. One major QTL for fruit blush was found on linkage group 4 (LG4) at 40-46 cM that explained from 20 to 32% of the total phenotypic variance and showed three QTL alleles of different effects. For soluble solids concentration (SSC), one QTL was mapped on LG5 at 60-72 cM and explained from 17 to 39% of the phenotypic variance. A major QTL for titratable acidity (TA) co-localized with the major locus for low-acid fruit (D-locus). It was mapped at the proximal end of LG5 and explained 35 to 80% of the phenotypic variance. The new QTL for TA on the distal end of LG5 explained 14 to 22% of the phenotypic variance. This QTL co-localized with the QTL for SSC and affected TA only when the first QTL is homozygous for high acidity (epistasis). Haplotype analyses revealed SNP haplotypes and predictive SNP marker(s) associated with desired QTL alleles.

CONCLUSIONS

A multi-family-based QTL discovery approach enhanced the ability to discover a new TA QTL at the distal end of LG5 and validated other QTLs which were reported in previous studies. Haplotype characterization of the mapped QTLs distinguishes this work from the previous QTL studies. Identified predictive SNPs and their original sources will facilitate the selection of parents and/or seedlings that have desired QTL alleles. Our findings will help peach breeders develop new predictive, DNA-based molecular marker tests for routine use in marker-assisted breeding.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a ‘SyntenyViewer’ to view genome synteny between different cucurbit species and an ‘RNA-Seq’ module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

Investigation of the transcriptomic and metabolic changes associated with superficial scald physiology impaired by lovastatin and 1-methylcyclopropene in pear fruit (cv. "Blanquilla")

To elucidate the physiology underlying the development of superficial scald in pears, susceptible "Blanquilla" fruit was treated with different compounds that either promoted (ethylene) or repressed (1-methylcyclopropene and lovastatin) the incidence of this disorder after 4 months of cold storage. Our data show that scald was negligible for the fruit treated with 1-methylcyclopropene or lovastatin, but highly manifested in untreated (78% incidence) or ethylene-treated fruit (97% incidence). The comparison between the fruit metabolomic profile and transcriptome evidenced a distinct reprogramming associated with each treatment. In all treated samples, cold storage led to an activation of a cold-acclimation-resistance mechanism, including the biosynthesis of very-long-chain fatty acids, which was especially evident in 1-methylcyclopropane-treated fruit. Among the treatments applied, only 1-methylcyclopropene inhibited ethylene production, hence supporting the involvement of this hormone in the development of scald. However, a common repression effect on the PPO gene combined with higher sorbitol content was found for both lovastatin and 1-methylcyclopropene-treated samples, suggesting also a non-ethylene-mediated process preventing the development of this disorder. The results presented in this work represent a step forward to better understand the physiological mechanisms governing the etiology of superficial scald in pears.

SpinachBase: a central portal for spinach genomics

Spinach (Spinacia oleracea L.) is a nutritious vegetable enriched with many essential minerals and vitamins. A reference spinach genome has been recently released, and additional spinach genomic resources are being rapidly developed. Therefore, there is an urgent need of a central database to store, query, analyze and integrate various resources of spinach genomic data. To this end, we developed SpinachBase (http://spinachbase.org), which provides centralized public accesses to genomic data as well as analytical tools to assist research and breeding in spinach. The database currently stores the spinach reference genome sequence, and sequences and comprehensive functional annotations of protein-coding genes predicted from the genome. The database also contains gene expression profiles derived from RNA-Seq experiments as well as highly co-expressed genes and genetic variants called from transcriptome sequences of 120 cultivated and wild Spinacia accessions. Biochemical pathways have been predicted from spinach protein-coding genes and are available through a pathway database (SpinachCyc) within SpinachBase. SpinachBase provides a suite of analysis and visualization tools including a genome browser, sequence similarity searches with BLAST, functional enrichment and functional classification analyses and functions to query and retrieve gene sequences and annotations.

Comparative genomic analysis of the PAL genes in five Rosaceae species and functional identification of Chinese white pear

Phenylalanine ammonia lyase (PAL) plays an important role in the biosynthesis of secondary metabolites regulating plant growth response. To date, the evolutionary history of the PAL family in Rosaceae plants remains unclear. In this study, we identified 16 PAL homologous genes in five Rosaceae plants (Pyrus bretschneideri, Fragaria vesca, Prunus mume, Prunus persica, and Malus × domestica). We classified these PALs into three categories based on phylogenetic analysis, and all PALs were distributed on 13 chromosomes. We tracked gene duplication events and performed sliding window analysis. These results revealed the evolution of PALs in five Rosaceae plants. We predicted the promoter of the PbPALs by PLANT CARE online software, and found that the promoter region of both PbPAL1 and PbPAL3 have at least one AC element. The results of qRT-PCR analysis found that PbPAL1 and PbPAL2 were highly expressed in the stems and roots, while expression level of PbPAL3 was relatively low in different tissues. The expression of PbPAL1 and PbPAL2 increased firstly and then decreased at different developmental periods of pear fruit. Among them, the expression of PbPAL1 reached the highest level 55 days after flowering. Three PbPALs were induced by abiotic stress to varying degrees. We transfected PbPAL1 and PbPAL2 into Arabidopsis thaliana, which resulted in an increase in lignin content and thickening of the cell walls of intervascular fibres and xylem cells. In summary, this research laid a foundation for better understanding the molecular evolution of PALs in five Rosaceae plants. Furthermore, the present study revealed the role of PbPALs in lignin synthesis, and provided basic data for regulating lignin synthesis and stone cells development in pear plants.

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.

Association of the phenylpropanoid pathway with dormancy and adaptive trait variation in apricot (Prunus armeniaca)

Trees use many mechanisms to adapt and respond to stressful conditions. The phenylpropanoid pathway in particular is known to be associated with a diverse suite of plant stress responses. In this study, we explored the relationship between the phenylpropanoid pathway metabolite production, gene expression and adaptive trait variation associated with floral bud reactivation during and following dormancy in Prunus armeniaca L. (apricot). Concentrations of eight phenylpropanoid metabolites were measured during chill accumulation and at developmental stages corresponding to the emergence of sepals and petals in floral buds of varieties that differ phenotypically in bloom date (BD). A significant interaction effect of chill hours and BD phenotype on the concentration of each of the compounds was observed (mixed analysis of variance, P < 0.05), with the concentration of most phenylpropanoid metabolites dropping precipitously when sepals and petals emerged. While phenylpropanoid biosynthetic gene expression patterns were more variable in general, expression changed over time and was impacted, although to a lesser degree, by BD phenotype. Furthermore, separation of BD phenotypic groups was most pronounced when early and late BD varieties were at different developmental stages, i.e., 800 chill hours. Taken together, these results suggest that the phenylpropanoid pathway is associated with floral bud reactivation in apricot. Furthermore, we show that the phenylpropanoid pathway is also impacted by phenological trait variation associated with dormancy. A better understanding of how apricot and other perennial tree species respond and adapt to environmental perturbations will be critical for improvement programs aimed at identifying and breeding trees more suitable for rapidly changing environments.

Cyberinfrastructure to Improve Forest Health and Productivity: The Role of Tree Databases in Connecting Genomes, Phenomes, and the Environment

Despite tremendous advancements in high throughput sequencing, the vast majority of tree genomes, and in particular, forest trees, remain elusive. Although primary databases store genetic resources for just over 2,000 forest tree species, these are largely focused on sequence storage, basic genome assemblies, and functional assignment through existing pipelines. The tree databases reviewed here serve as secondary repositories for community data. They vary in their focal species, the data they curate, and the analytics provided, but they are united in moving toward a goal of centralizing both data access and analysis. They provide frameworks to view and update annotations for complex genomes, interrogate systems level expression profiles, curate data for comparative genomics, and perform real-time analysis with genotype and phenotype data. The organism databases of today are no longer simply catalogs or containers of genetic information. These repositories represent integrated cyberinfrastructure that support cross-site queries and analysis in web-based environments. These resources are striving to integrate across diverse experimental designs, sequence types, and related measures through ontologies, community standards, and web services. Efficient, simple, and robust platforms that enhance the data generated by the research community, contribute to improving forest health and productivity.

Dryas as a Model for Studying the Root Symbioses of the Rosaceae

The nitrogen-fixing root nodule symbiosis is restricted to four plant orders: Fabales (legumes), Fagales, Cucurbitales and Rosales (Elaeagnaceae, Rhamnaceae, and Rosaceae). Interestingly all of the Rosaceae genera confirmed to contain nodulating species (i.e., Cercocarpus, Chamaebatia, Dryas, and Purshia) belong to a single subfamily, the Dryadoideae. The Dryas genus is particularly interesting from an evolutionary perspective because it contains closely related nodulating (Dryas drummondii) and non-nodulating species (Dryas octopetala). The close phylogenetic relationship between these two species makes Dryas an ideal model genus to study the genetic basis of nodulation by whole genome comparison and classical genetics. Therefore, we established methods for plant cultivation, transformation and DNA extraction for these species. We optimized seed surface sterilization and germination methods and tested growth protocols ranging from pots and Petri dishes to a hydroponic system. Transgenic hairy roots were obtained by adapting Agrobacterium rhizogenes-based transformation protocols for Dryas species. We compared several DNA extraction protocols for their suitability for subsequent molecular biological analysis. Using CTAB extraction, reproducible PCRs could be performed, but CsCl gradient purification was essential to obtain DNA in sufficient purity for high quality de novo genome sequencing of both Dryas species. Altogether, we established a basic toolkit for the culture, transient transformation and genetic analysis of Dryas sp.

Exploring sources of resistance to brown rot in an interspecific almond × peach population

BACKGROUND

Monilinia spp. are responsible for brown rot, one of the most significant stone fruit diseases. Planting resistant cultivars seems a promising alternative, although most commercial cultivars are susceptible to brown rot. The aim of this study was to explore resistance to Monilinia fructicola over two seasons in a backcross one interspecific population between almond 'Texas' and peach 'Earlygold' (named T1E).

RESULTS

'Texas' almond was resistant to brown rot inoculation, whereas peach was highly susceptible. Phenotypic data from the T1E population indicated wide differences in response to M. fructicola. Additionally, several non-wounded individuals exhibited resistance to brown rot. Quantitative trait loci (QTLs) were identified in several linkage groups, but only two proximal QTLs in G4 were detected over both seasons and accounted for 11.3-16.2% of the phenotypic variation.

CONCLUSION

Analysis of the progeny allowed the identification of resistant genotypes that could serve as a source of resistance in peach breeding programs. The finding of loci associated with brown rot resistance would shed light on implementing a strategy based on marker-assisted selection (MAS) for introgression of this trait into elite peach materials. New peach cultivars resistant to brown rot may contribute to the implementation of more sustainable crop protection strategies. © 2019 Society of Chemical Industry.

Genome-Wide Association Studies in Apple Reveal Loci for Aroma Volatiles, Sugar Composition, and Harvest Date

Understanding the genetic architecture of fruit quality traits is crucial to target breeding of apple ( L.) cultivars. We linked genotype and phenotype information by combining genotyping-by-sequencing (GBS) generated single nucleotide polymorphism (SNP) markers with fruit flavor volatile data, sugar and acid content, and historical trait data from a gene bank collection. Using gas chromatography-mass spectrometry (GC-MS) analysis of apple juice samples, we identified 49 fruit volatile organic compounds (VOCs). We found a very variable content of VOCs, especially for the esters, among 149 apple cultivars. We identified convincing associations for the acetate esters especially butyl acetate and hexyl acetate on chromosome 2 in a region of several alcohol acyl-transferases including AAT1. For sucrose content and for fructose and sucrose in percentage of total sugars, we revealed significant SNP associations. Here, we suggest a vacuolar invertase close to significant SNPs for this association as candidate gene. Harvest date was in strong SNP association with a NAC transcription factor gene and sequencing identified two haplotypes associated with harvest date. The study shows that SNP marker characterization of a gene bank collection can be successfully combined with new and historical trait data for association studies. Suggested candidate genes may contribute to an improved understanding of the genetic basis for important traits and simultaneously provide tools for targeted breeding using marker-assisted selection (MAS).

Transcriptional regulation of MdPIN3 and MdPIN10 by MdFLP during apple self-rooted stock adventitious root gravitropism

BACKGROUND

The close planting of dwarfing self-rooted rootstocks is currently a widely used method for apple production; however, self-rooted rootstocks are weak with shallow roots and poor grounding. Therefore, understanding the molecular mechanisms that establish the gravitropic set-point angles (GSAs) of the adventitious roots of self-rooted apple stocks is important for developing self-rooted apple rootstock cultivars with deep roots.

RESULTS

We report that the apple FOUR LIPS (MdFLP), an R2R3-MYB transcription factor (TF), functions in establishing the GSA of the adventitious roots of self-rooted apple stocks in response to gravity. Biochemical analyses demonstrate that MdFLP directly binds to the promoters of two auxin efflux carriers, MdPIN3 and MdPIN10, that are involved in auxin transport, activates their transcriptional expression, and thereby promotes the development of adventitious roots in self-rooted apple stocks. Additionally, the apple auxin response factor MdARF19 influences the expression of those auxin efflux carriers and the establishment of the GSA of adventitious roots of apple in response to gravity by directly activating the expression of MdFLP.

CONCLUSION

Our findings provide new insights into the transcriptional regulation of MdFLP by the auxin response factor MdARF19 in the regulation of the GSA of adventitious roots of self-rooted apple stocks in response to gravity.

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.

Edit at will: Genotype independent plant transformation in the era of advanced genomics and genome editing

The combination of advanced genomics, genome editing and plant transformation biology presents a powerful platform for basic plant research and crop improvement. Together these advances provide the tools to identify genes as targets for direct editing as single base pair changes, deletions, insertions and site specific homologous recombination. Recent breakthrough technologies using morphogenic regulators in plant transformation creates the ability to introduce reagents specific toward their identified targets and recover stably transformed and/or edited plants which are genotype independent. These technologies enable the possibility to alter a trait in any variety, without genetic disruption which would require subsequent extensive breeding, but rather to deliver the same variety with one trait changed. Regulatory issues regarding this technology will predicate how broadly these technologies will be implemented. In addition, education will play a crucial role for positive public acceptance. Taken together these technologies comprise a platform for advanced breeding which is an imperative for future world food security.

Indole-3-acetylaspartate and indole-3-acetylglutamate, the IAA-amide conjugates in the diploid strawberry achene, are hydrolyzed in growing seedlings

Indole-3-acetylaspartate and indole-3-acetylglutamate are the stored auxin amino acid conjugates of the achene of the diploid strawberry and serve as sources of auxin during seedling growth. The edible part of the strawberry, a pseudocarp, has long been known to enlarge in response to auxin produced by the developing achenes, the botanical true fruit. Auxin homeostasis involves a complex interaction between biosynthesis, conjugate formation and hydrolysis, catabolism and transport. Strawberry tissues are capable of synthesizing auxin conjugates, and transcriptome data support the expression of genes involved in IAA conjugate formation and hydrolysis throughout embryo development and subsequent seedling growth. Using a highly sensitive and selective mass spectrometric method, we identified all the low molecular weight indole-auxin amino acid conjugates in achenes of F. vesca as consisting of indole-3-acetylaspartate (IAasp) and indole-3-acetylglutamate (IAglu). In contrast to what has been proposed to occur in Arabidopsis, we determined that IAasp and IAglu are hydrolyzed by seedlings to provide a source of free IAA for growth.

The complete sequences of two divergent variants of the rhabdovirus raspberry vein chlorosis virus and the design of improved primers for virus detection

The complete sequence was obtained for two variants of raspberry vein chlorosis virus (RVCV), confirming that this virus is a rhabdovirus most closely related to the cytorhabdoviruses alfalfa dwarf virus and strawberry crinkle virus. The two RVCV variants share only a 68% nucleotide sequence identity so that the previously published RT-PCR diagnostic test for this virus was not able to efficiently detect both variants. Using the new, complete sequence information several new primer sets have been designed that allow a much improved RVCV detection.

Wide transcriptional investigation unravel novel insights of the on-tree maturation and postharvest ripening of 'Abate Fetel' pear fruit

To decipher the transcriptomic regulation of the on-tree fruit maturation in pear cv. 'Abate Fetel', a RNA-seq transcription analysis identified 8939 genes differentially expressed across four harvesting stages. These genes were grouped into 11 SOTA clusters based on their transcriptional pattern, of which three included genes upregulated while the other four were represented by downregulated genes. Fruit ripening was furthermore investigated after 1 month of postharvest cold storage. The most important variation in fruit firmness, production of ethylene and volatile organic compounds were observed after 5 days of shelf-life at room temperature following cold storage. The role of ethylene in controlling the ripening of 'Abate Fetel' pears was furthermore investigated through the application of 1-methylcyclopropene, which efficiently delayed the progression of ripening by reducing fruit softening and repressing both ethylene and volatile production. The physiological response of the interference at the ethylene receptor level was moreover unraveled investigating the expression pattern of 12 candidate genes, initially selected to validate the RNA-seq profile. This analysis confirmed the effective role of the ethylene competitor in downregulating the expression of cell wall (PG) and ethylene-related genes (ACS, ACO, ERS1, and ERS2), as well as inducing one element involved in the auxin signaling pathway (Aux/IAA), highlighting a possible cross-talk between these two hormones. The expression patterns of these six elements suggest their use as molecular toolkit to monitor at molecular level the progression of the fruit on-tree maturation and postharvest ripening.

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

BACKGROUND

The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention.

RESULTS

We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database.

CONCLUSIONS

We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

Translational genomics and multi-omics integrated approaches as a useful strategy for crop breeding

Recent next generation sequencing-driven mass production of genomic data and multi-omics-integrated approaches have significantly contributed to broadening and deepening our knowledge on the molecular system of living organisms. Accordingly, translational genomics (TG) approach can play a pivotal role in creating an informational bridge between model systems and relatively less studied plants. This review focuses mainly on addressing recent advancement in omics-related technologies, a diverse array of bioinformatic resources and potential applications of TG for the crop breeding. To accomplish above objectives, information on omics data production, various DBs and high throughput technologies was collected, integrated, and used to analyze current status and future perspectives towards omics-assisted crop breeding. Various omics data and resources have been organized and integrated into the databases and/or bioinformatic infrastructures, and thereby serve as the ome's information center for cross-genome translation of biological data. Although the size of accumulated omics data and availability of reference genomes are different among plant families, translational approaches have been actively progressing to access particular biological characteristics. When multi-layered omics data are integrated in a synthetic manner, it will allow providing a stereoscopic view of dynamic molecular behavior and interacting networks of genes occurring in plants. Consequently, TG approach will lead us to broader and deeper insights into target traits for the plant breeding. Furthermore, such systems approach will renovate conventional breeding programs and accelerate precision crop breeding in the future.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a 'SyntenyViewer' to view genome synteny between different cucurbit species and an 'RNA-Seq' module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

15 years of GDR: New data and functionality in the Genome Database for Rosaceae

The Genome Database for Rosaceae (GDR, https://www.rosaceae.org) is an integrated web-based community database resource providing access to publicly available genomics, genetics and breeding data and data-mining tools to facilitate basic, translational and applied research in Rosaceae. The volume of data in GDR has increased greatly over the last 5 years. The GDR now houses multiple versions of whole genome assembly and annotation data from 14 species, made available by recent advances in sequencing technology. Annotated and searchable reference transcriptomes, RefTrans, combining peer-reviewed published RNA-Seq as well as EST datasets, are newly available for major crop species. Significantly more quantitative trait loci, genetic maps and markers are available in MapViewer, a new visualization tool that better integrates with other pages in GDR. Pathways can be accessed through the new GDR Cyc Pathways databases, and synteny among the newest genome assemblies from eight species can be viewed through the new synteny browser, SynView. Collated single-nucleotide polymorphism diversity data and phenotypic data from publicly available breeding datasets are integrated with other relevant data. Also, the new Breeding Information Management System allows breeders to upload, manage and analyze their private breeding data within the secure GDR server with an option to release data publicly.

Tripal MapViewer: A tool for interactive visualization and comparison of genetic maps

Tripal is an open-source, resource-efficient toolkit for construction of genomic, genetic and breeding databases. It facilitates development of biological websites by providing tools to integrate and display biological data using the generic database schema, Chado, together with Drupal, a popular website creation and content management system. Tripal MapViewer is a new interactive tool for visualizing genetic map data. Developed as a Tripal replacement for Comparative Map Viewer (CMap), it enables visualization of entire maps or linkage groups and features such as molecular markers, quantitative trait loci (QTLs) and heritable phenotypic markers. It also provides graphical comparison of maps sharing the same markers as well as dot plot and correspondence matrices. MapViewer integrates directly with the Tripal application programming interface framework, improving data searching capability and providing a more seamless experience for site visitors. The Tripal MapViewer interface can be integrated in any Tripal map page and linked from any Tripal page for markers, QTLs, heritable morphological markers or genes. Configuration of the display is available through a control panel and the administration interface. The administration interface also allows configuration of the custom database query for building materialized views, providing better performance and flexibility in the way data is stored in the Chado database schema. MapViewer is implemented with the D3.js technology and is currently being used at the Genome Database for Rosaceae (https://www.rosaceae.org), CottonGen (https://www.cottongen.org), Citrus Genome Database (https://citrusgenomedb.org), Vaccinium Genome Database (https://www.vaccinium.org) and Cool Season Food Legume Database (https://www.coolseasonfoodlegume.org). It is also currently in development on the Hardwood Genomics Web (https://hardwoodgenomics.org) and TreeGenes (https://treegenesdb.org). Database URL: https://gitlab.com/mainlabwsu/tripal_map.

Genetic characterization of worldwide Prunus domestica (plum) germplasm using sequence-based genotyping

Prunus domestica commonly known as European plum is a hexaploid fruit tree species cultivated around the world. Locally it is used for fresh consumption, in jams or jellies, and the production of spirits while commercially the fruit is primarily sold dried (prunes). Despite its agricultural importance and long history of cultivation, many questions remain about the origin of this species, the relationships among its many pomological types, and its underlying genetics. Here, we used a sequence-based genotyping approach to characterize worldwide plum germplasm including the potential progenitor Eurasian plum species. Analysis of 405 DNA samples established a set of four clades consistent with the pomological groups Greengages, Mirabelles, European plums, and d'Agen (French) prune plums. A number of cultivars from each clade were identified as likely clonal selections, particularly among the "French" type prune germplasm that is widely cultivated today. Overall, there was relatively low genetic diversity across all cultivated plums suggesting they have been largely inbred and/or derived from a limited number of founders. The results agree with P. domestica having originated as an interspecific hybrid of a diploid P. cerasifera and a tetraploid P. spinosa that itself may have been an interspecific hybrid of P. cerasifera and an unknown Eurasian plum species. The low genetic diversity and lack of true wild-types coupled with the known cultivation history of Eurasian plums imply that P. domestica may have been a product of inter-specific cross breeding and artificial selection by early agrarian Eurasian societies.

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.

Genome sequences of horticultural plants: past, present, and future

Horticultural plants play various and critical roles for humans by providing fruits, vegetables, materials for beverages, and herbal medicines and by acting as ornamentals. They have also shaped human art, culture, and environments and thereby have influenced the lifestyles of humans. With the advent of sequencing technologies, there has been a dramatic increase in the number of sequenced genomes of horticultural plant species in the past decade. The genomes of horticultural plants are highly diverse and complex, often with a high degree of heterozygosity and a high ploidy due to their long and complex history of evolution and domestication. Here we summarize the advances in the genome sequencing of horticultural plants, the reconstruction of pan-genomes, and the development of horticultural genome databases. We also discuss past, present, and future studies related to genome sequencing, data storage, data quality, data sharing, and data visualization to provide practical guidance for genomic studies of horticultural plants. Finally, we propose a horticultural plant genome project as well as the roadmap and technical details toward three goals of the project.

Identification of QTLs for powdery mildew (Podosphaera aphanis; syn. Sphaerotheca macularis f. sp. fragariae) susceptibility in cultivated strawberry (Fragaria ×ananassa)

Strawberry powdery mildew (Podosphaera aphanis Wallr.) is a pathogen which infects the leaves, fruit, stolon and flowers of the cultivated strawberry (Fragaria ×ananassa), causing major yield losses, primarily through unmarketable fruit. The primary commercial control of the disease is the application of fungicidal sprays. However, as the use of key active ingredients of commercial fungicides is becoming increasingly restricted, interest in developing novel strawberry cultivars exhibiting resistance to the pathogen is growing rapidly. In this study, a mapping population derived from a cross between two commercial strawberry cultivars ('Sonata' and 'Babette') was genotyped with single nucleotide polymorphism (SNP) markers from the Axiom iStraw90k genotyping array and phenotyped for powdery mildew susceptibility in both glasshouse and field environments. Three distinct, significant QTLs for powdery mildew resistance were identified across the two experiments. Through comparison with previous studies and scrutiny of the F. vesca genome sequence, candidate genes underlying the genetic control of this trait were identified.

Cross-Kingdom Commonality of a Novel Insertion Signature of RTE-Related Short Retroposons

In multicellular organisms, such as vertebrates and flowering plants, horizontal transfer (HT) of genetic information is thought to be a rare event. However, recent findings unveiled unexpectedly frequent HT of RTE-clade LINEs. To elucidate the molecular footprints of the genomic integration machinery of RTE-related retroposons, the sequence patterns surrounding the insertion sites of plant Au-like SINE families were analyzed in the genomes of a wide variety of flowering plants. A novel and remarkable finding regarding target site duplications (TSDs) for SINEs was they start with thymine approximately one helical pitch (ten nucleotides) downstream of a thymine stretch. This TSD pattern was found in RTE-clade LINEs, which share the 3'-end sequence of these SINEs, in the genome of leguminous plants. These results demonstrably show that Au-like SINEs were mobilized by the enzymatic machinery of RTE-clade LINEs. Further, we discovered the same TSD pattern in animal SINEs from lizard and mammals, in which the RTE-clade LINEs sharing the 3'-end sequence with these animal SINEs showed a distinct TSD pattern. Moreover, a significant correlation was observed between the first nucleotide of TSDs and microsatellite-like sequences found at the 3'-ends of SINEs and LINEs. We propose that RTE-encoded protein could preferentially bind to a DNA region that contains a thymine stretch to cleave a phosphodiester bond downstream of the stretch. Further, determination of cleavage sites and/or efficiency of primer sites for reverse transcription may depend on microsatellite-like repeats in the RNA template. Such a unique mechanism may have enabled retroposons to successfully expand in frontier genomes after HT.

Identification of a strawberry NPR-like gene involved in negative regulation of the salicylic acid-mediated defense pathway

Hormonal modulation plays a central role in triggering various resistant responses to biotic and abiotic stresses in plants. In cultivated strawberry (Fragaria x ananassa), the salicylic acid (SA)-dependent defense pathway has been associated with resistance to Colletotrichum spp. and the other pathogens. To better understand the SA-mediated defense mechanisms in strawberry, we analyzed two strawberry cultivars treated with SA for their resistance to anthracnose and gene expression profiles at 6, 12, 24, and 48 hr post-treatment. Strawberry genes related to SA biosynthesis, perception, and signaling were identified from SA-responsive transcriptomes of the two cultivars, and the induction of 17 candidate genes upon SA treatment was confirmed by qRT-PCR. Given the pivotal role of the non-expressor of pathogenesis-related (NPR) family in controlling the SA-mediated defense signaling pathway, we then analyzed NPR orthologous genes in strawberry. From the expression profile, FaNPRL-1 [ortholog of FvNPRL-1 (gene20070 in F. vesca)] was identified as an NPR-like gene significantly induced after SA treatment in both cultivars. With a conserved BTB/POZ domain, ankyrin repeat domain, and nuclear localization signal, FvNPRL-1 was found phylogenetically closer to NPR3/NPR4 than NPR1 in Arabidopsis. Ectopic expression of FvNPRL-1 in the Arabidopsis thaliana wild type suppressed the SA-mediated PR1 expression and the resistance to Pseudomonas syringae pv. tomato DC3000. Transient expression of FvNPRL-1 fused with green fluorescent protein in Arabidopsis protoplasts showed that SA affected nuclear translocation of FvNPRL-1. FvNPRL-1 likely functions similar to Arabidopsis NPR3/NPR4 as a negative regulator of the SA-mediated defense.

Repeated translocation of a gene cassette drives sex-chromosome turnover in strawberries

Turnovers of sex-determining systems represent important diversifying forces across eukaryotes. Shifts in sex chromosomes—but conservation of the master sex-determining genes—characterize distantly related animal lineages. Yet in plants, in which separate sexes have evolved repeatedly and sex chromosomes are typically homomorphic, we do not know whether such translocations drive sex-chromosome turnovers within closely related taxonomic groups. This phenomenon can only be demonstrated by identifying sex-associated nucleotide sequences, still largely unknown in plants. The wild North American octoploid strawberries (Fragaria) exhibit separate sexes (dioecy) with homomorphic, female heterogametic (ZW) inheritance, yet sex maps to three different chromosomes in different taxa. To characterize these turnovers, we identified sequences unique to females and assembled their reads into contigs. For most octoploid Fragaria taxa, a short (13 kb) sequence was observed in all females and never in males, implicating it as the sex-determining region (SDR). This female-specific “SDR cassette” contains both a gene with a known role in fruit and pollen production and a novel retrogene absent on Z and autosomal chromosomes. Phylogenetic comparison of SDR cassettes revealed three clades and a history of repeated translocation. Remarkably, the translocations can be ordered temporally due to the capture of adjacent sequence with each successive move. The accumulation of the “souvenir” sequence—and the resultant expansion of the hemizygous SDR over time—could have been adaptive by locking genes into linkage with sex. Terminal inverted repeats at the insertion borders suggest a means of movement. To our knowledge, this is the first plant SDR shown to be translocated, and it suggests a new mechanism (“move-lock-grow”) for expansion and diversification of incipient sex chromosomes.

Sex chromosomes frequently restructure themselves during organismal evolution, often becoming highly differentiated. This dynamic process is poorly understood for most taxa, especially during the early stages typical of many dioecious flowering plants. We show that in wild strawberries, a female-specific region of DNA is associated with sex and has repeatedly changed its genomic location, each time increasing the size of the hemizygous female-specific sequence on the W sex chromosome. This observation shows, for the first time to our knowledge, that plant sex regions can “jump” and suggests that this phenomenon may be adaptive by gathering and locking new genes into linkage with sex. This conserved and presumed causal sex-determining sequence, which varies in both genomic location and degree of differentiation, will facilitate future studies to understand how sex chromosomes first begin to differentiate.