The draft genome sequence of European pear (Pyrus communis L. 'Bartlett')

Genetic and molecular regulation of chilling requirements in pear: breeding for climate change resilience

Pear (Pyrus spp.) is a deciduous fruit tree that requires exposure to sufficient chilling hours during the winter to establish dormancy, followed by favorable heat conditions during the spring for normal vegetative and floral budbreak. In contrast to most temperate woody species, apples and pears of the Rosaceae family are insensitive to photoperiod, and low temperature is the major factor that induces growth cessation and dormancy. Most European pear (Pyrus Communis L.) cultivars need to be grown in regions with high chilling unit (CU) accumulation to ensure early vegetative budbreak. Adequate vegetative budbreak time will ensure suitable metabolite accumulation, such as sugars, to support fruit set and vegetative development, providing the necessary metabolites for optimal fruit set and development. Many regions that were suitable for pear production suffer from a reduction in CU accumulation. According to climate prediction models, many temperate regions currently suitable for pear cultivation will experience a similar accumulation of CUs as observed in Mediterranean regions. Consequently, the Mediterranean region can serve as a suitable location for conducting pear breeding trials aimed at developing cultivars that will thrive in temperate regions in the decades to come. Due to recent climatic changes, bud dormancy attracts more attention, and several studies have been carried out aiming to discover the genetic and physiological factors associated with dormancy in deciduous fruit trees, including pears, along with their related biosynthetic pathways. In this review, current knowledge of the genetic mechanisms associated with bud dormancy in European pear and other Pyrus species is summarized, along with metabolites and physiological factors affecting dormancy establishment and release and chilling requirement determination. The genetic and physiological insights gained into the factors regulating pear dormancy phase transition and determining chilling requirements can accelerate the development of new pear cultivars better suited to both current and predicted future climatic conditions.

A chromosome-scale assembly for 'd'Anjou' pear

Cultivated pear consists of several Pyrus species with Pyrus communis (European pear) representing a large fraction of worldwide production. As a relatively recently domesticated crop and perennial tree, pear can benefit from genome-assisted breeding. Additionally, comparative genomics within Rosaceae promises greater understanding of evolution within this economically important family. Here, we generate a fully phased chromosome-scale genome assembly of P. communis 'd'Anjou.' Using PacBio HiFi and Dovetail Omni-C reads, the genome is resolved into the expected 17 chromosomes, with each haplotype totaling nearly 540 Megabases and a contig N50 of nearly 14 Mb. Both haplotypes are highly syntenic to each other and to the Malus domestica 'Honeycrisp' apple genome. Nearly 45,000 genes were annotated in each haplotype, over 90% of which have direct RNA-seq expression evidence. We detect signatures of the known whole-genome duplication shared between apple and pear, and we estimate 57% of d'Anjou genes are retained in duplicate derived from this event. This genome highlights the value of generating phased diploid assemblies for recovering the full allelic complement in highly heterozygous crop species.

Introduction of a diverse genetic background of Pyrus into Malus through intergeneric hybridization

Wide hybridizations across species and genera have been employed to enhance agriculturally important traits in crops. Within the tribe Maleae of the Rosaceae family, different genera and species exhibit several traits useful for increasing diversity and gene pool through hybridization. This study aimed to develop and characterize intergeneric hybrid individuals between Malus and Pyrus. Through seed germination, shoot multiplication, and rooting in vitro, acclimatized seedlings showing vegetative growth on their own roots were obtained from crosses of Malus × domestica pollinated by Pyrus communis, P. bretschneideri, and the Pyrus interspecific hybrid (P. communis × P. pyrifolia). Comparative analysis of leaf morphology, flow cytometry, and molecular genotyping confirmed the hybrid status of the individuals. Genome-wide genotyping revealed that all the hybrid individuals inherited genomic fragments symmetrically from the Malus and Pyrus parents. To the best of our knowledge, this is the first report on the development of intergeneric hybrid seedlings between Malus × domestica and P. bretschneideri. Furthermore, the Pyrus interspecific hybrid individual served as a bridge plant for introducing the genetic background of P. pyrifolia into Malus × domestica. The results of this study provided a crucial foundation for breeding through intergeneric hybridization between Malus and Pyrus, facilitating the incorporation of valuable traits from diverse gene pools.

PbrWRKY62-PbrADC1 module involves in superficial scald development of Pyrus bretschneideri Rehd.fruit via regulating putrescine biosynthesis

Putrescine plays a role in superficial scald development during the cold storage of pear fruit. However, the molecular mechanism behind this phenomenon has not been un-fully clarified until recently. In this study, a conjoint analysis of metabolites and gene expression profiles in the putrescine-metabolic pathway of P. bretschneideri Rehd. fruit followed by experimental validation revealed that PbrADC1, forming a homodimer in the chloroplast, was involved in putrescine biosynthesis and thus fruit chilling resistance. Additionally, the substrate-binding residue Cys546 in PbrADC1, whose activity was modified by H2O2, played a crucial role in arginine decarboxylation into agmatine. Through a combined analysis of the distribution of cis-acting elements in the PbrADC1 promoter as well as the expression profiles of related transcription factors (TFs), several TFs were identified as upstream regulators of PbrADC1 gene. Further investigation revealed that the nuclear PbrWRKY62 could directly bind to the W-box elements in the PbrADC1 promoter, activate its expression, enhance putrescine accumulation, and thus increase fruit chilling tolerance. In conclusion, our results suggest that the PbrWRKY62-PbrADC1 module is involved in the development of superficial scald in P. bretschneideri Rehd. fruit via regulating putrescine biosynthesis. Consequently, these findings could serve as valuable genetic resources for breeding scald-resistant pear fruit.

Telomere-to-telomere pear (Pyrus pyrifolia) reference genome reveals segmental and whole genome duplication driving genome evolution

Previously released pear genomes contain a plethora of gaps and unanchored genetic regions. Here, we report a telomere-to-telomere (T2T) gap-free genome for the red-skinned pear, 'Yunhong No. 1' (YH1; Pyrus pyrifolia), which is mainly cultivated in Yunnan Province (southwest China), the pear's primary region of origin. The YH1 genome is 501.20 Mb long with a contig N50 length of 29.26 Mb. All 17 chromosomes were assembled to the T2T level with 34 characterized telomeres. The 17 centromeres were predicted and mainly consist of centromeric-specific monomers (CEN198) and long terminal repeat (LTR) Gypsy elements (≥74.73%). By filling all unclosed gaps, the integrity of YH1 is markedly improved over previous P. pyrifolia genomes ('Cuiguan' and 'Nijisseiki'). A total of 1531 segmental duplication (SD) driven duplicated genes were identified and enriched in stress response pathways. Intrachromosomal SDs drove the expansion of disease resistance genes, suggesting the potential of SDs in adaptive pear evolution. A large proportion of duplicated gene pairs exhibit dosage effects or sub-/neo-functionalization, which may affect agronomic traits like stone cell content, sugar content, and fruit skin russet. Furthermore, as core regulators of anthocyanin biosynthesis, we found that MYB10 and MYB114 underwent various gene duplication events. Multiple copies of MYB10 and MYB114 displayed obvious dosage effects, indicating role differentiation in the formation of red-skinned pear fruit. In summary, the T2T gap-free pear genome provides invaluable resources for genome evolution and functional genomics.

Genome-wide identification, comparative analysis and functional roles in flavonoid biosynthesis of cytochrome P450 superfamily in pear (Pyrus spp.)

BACKGROUND

The cytochrome P450 (CYP) superfamily is the largest enzyme metabolism family in plants identified to date, and it is involved in many biological processes, including secondary metabolite biosynthesis, hormone metabolism and stress resistance. However, the P450 gene superfamily has not been well studied in pear (Pyrus spp.).

RESULTS

Here, the comprehensive identification and a comparative analysis of P450 superfamily members were conducted in cultivated and wild pear genomes. In total, 338, 299 and 419 P450 genes were identified in Chinese white pear, European pear and the wild pear, respectively. Based on the phylogenetic analyses, pear P450 genes were divided into ten clans, comprising 48 families. The motif and gene structure analyses further supported this classification. The expansion of the pear P450 gene family was attributed to whole-genome and single-gene duplication events. Several P450 gene clusters were detected, which have resulted from tandem and proximal duplications. Purifying selection was the major force imposed on the long-term evolution of P450 genes. Gene dosage balance, subfunctionalization and neofunctionalization jointly drove the retention and functional diversification of P450 gene pairs. Based on the association analysis between transcriptome expression profiles and flavonoid content during fruit development, three candidate genes were identified as being closely associated with the flavonoid biosynthesis, and the expression of one gene was further verified using qRT-PCR and its function was validated through transient transformation in pear fruit.

CONCLUSIONS

The study results provide insights into the evolution and biological functions of P450 genes in pear.

Stone cell formation in the pedicel of pears and apples

MAIN CONCLUSION

For the first time, stone cells in pear and apple pedicel were studied. The lignification of the pedicel outer part was correlated with flesh, and the secondary cell wall biosynthesis genes were activated. Fruit pedicels act as bridges between the fruit and the shoot. They have secondary thickened cell walls that presumably function in mechanical support, water and nutrient transport. Stone cells are cells with a secondary cell wall thickening. In pears, yet not in apples, the stone cells affect the flesh texture. There have been few reports on stone cell formation in pear and apple pedicels; therefore, we studied these cells for the first time. The apple pedicel had few stone cells in the cortex. The formation of stone cells in pear continued until seven weeks after flowering (WAF), and the density was significantly higher than in apple. The stone cell formation degree (SFD) of pear was 3.6-7.1 times higher than that of apple. Total lignin and lignin non-condensed structure (G and S units) content in the pear pedicle outer part was 1.5-2.7 times higher than that of the apple at harvest. The SFD of the pedicel outer part had a positive correlation with the G and S units content of the flesh. The total lignin and G and S units content between flesh and the pedicel outer part were positively correlated. Correlation analysis revealed a positive relationship between fruit and pedicel formation of the stone cells. The WGCNA showed that NST3 was linked to NAC028, MYB46, CESA, POD, LAC, and VSR6. These genes were highly expressed in the outer part of the pear pedicel, while they were suppressed in that issue of the apple at 4 WAF.

Identification and evolution of nsLTPs in the root nodule nitrogen fixation clade and molecular response of Frankia to AgLTP24

Non-specific lipid transfer proteins (nsLTPs) are antimicrobial peptides, involved in several plant biological processes including root nodule nitrogen fixation (RNF). Nodulating plants belonging to the RNF clade establish symbiosis with the nitrogen-fixing bacteria rhizobia (legumes symbiosis model) and Frankia (actinorhizal symbiosis model) leading to root nodule formation. nsLTPs are involved in processes active in early step of symbiosis and functional nodule in both models. In legumes, nsLTPs have been shown to regulate symbiont entry, promote root cortex infection, membrane biosynthesis, and improve symbiosis efficiency. More recently, a nsLTP, AgLTP24 has been described in the context of actinorhizal symbiosis between Alnus glutinosa and Frankia alni ACN14a. AgLTP24 is secreted at an early step of symbiosis on the deformed root hairs and targets the symbiont in the nitrogen-fixing vesicles in functional nodules. nsLTPs are involved in RNF, but their functions and evolutionary history are still largely unknown. Numerous putative nsLTPs were found up-regulated in functional nodules compared to non-infected roots in different lineages within the RNF clade. Here, results highlight that nodulating plants that are co-evolving with their nitrogen-fixing symbionts appear to have independently specialized nsLTPs for this interaction, suggesting a possible convergence of function, which opens perspectives to investigate nsLTPs functions in RNF.

A MADS-box gene-induced early flowering pear (Pyrus communis L.) for accelerated pear breeding

There have been a considerable number of studies that have successfully sped up the flowering cycle in woody perennial horticultural species. One particularly successful study in apple (Malus domestica) accelerated flowering using a silver birch (Betula pendula) APETALA1/FRUITFULL MADS-box gene BpMADS4, which yielded a good balance of vegetative growth to support subsequent flower and fruit development. In this study, BpMADS4 was constitutively expressed in European pear (Pyrus communis) to establish whether this could be used as a tool in a rapid pear breeding program. Transformed pear lines flowered within 6-18 months after grafting onto a quince (Cydonia oblonga) rootstock. Unlike the spindly habit of early flowering apples, the early flowering pear lines displayed a normal tree-like habit. Like apple, the flower appearance was normal, and the flowers were fertile, producing fruit and seed upon pollination. Seed from these transformed lines were germinated and 50% of the progeny flowered within 3 months of sowing, demonstrating a use for these in a fast breeding program.

The pear genomics database (PGDB): a comprehensive multi-omics research platform for Pyrus spp

BACKGROUND

Pears are among the most important temperate fruit trees in the world, with significant research efforts increasing over the last years. However, available omics data for pear cannot be easily and quickly retrieved to enable further studies using these biological data.

DESCRIPTION

Here, we present a publicly accessible multi-omics pear resource platform, the Pear Genomics Database (PGDB). We collected and collated data on genomic sequences, genome structure, functional annotation, transcription factor predictions, comparative genomics, and transcriptomics. We provide user-friendly functional modules to facilitate querying, browsing and usage of these data. The platform also includes basic and useful tools, including JBrowse, BLAST, phylogenetic tree building, and additional resources providing the possibility for bulk data download and quick usage guide services.

CONCLUSIONS

The Pear Genomics Database (PGDB, http://pyrusgdb.sdau.edu.cn ) is an online data analysis and query resource that integrates comprehensive multi-omics data for pear. This database is equipped with user-friendly interactive functional modules and data visualization tools, and constitutes a convenient platform for integrated research on pear.

The East Asian wild apples, Malus baccata (L.) Borkh and Malus hupehensis (Pamp.) Rehder., are additional contributors to the genomes of cultivated European and Chinese varieties

The domestication process in long-lived plant perennials differs dramatically from that of annuals, with a huge amount of genetic exchange between crop and wild populations. Though apple is a major fruit crop grown worldwide, the contribution of wild apple species to the genetic makeup of the cultivated apple genome remains a topic of intense study. We used population genomics approaches to investigate the contributions of several wild apple species to European and Chinese rootstock and dessert genomes, with a focus on the extent of wild-crop gene flow. Population genetic structure inferences revealed that the East Asian wild apples, Malus baccata (L.) Borkh and M. hupehensis (Pamp.), form a single panmictic group, and that the European dessert and rootstock apples form a specific gene pool whereas the Chinese dessert and rootstock apples were a mixture of three wild gene pools, suggesting different evolutionary histories of European and Chinese apple varieties. Coalescent-based inferences and gene flow estimates indicated that M. baccata - M. hupehensis contributed to the genome of both European and Chinese cultivated apples through wild-to-crop introgressions, and not as an initial contributor as previously supposed. We also confirmed the contribution through wild-to-crop introgressions of Malus sylvestris Mill. to the cultivated apple genome. Apple tree domestication is therefore one example in woody perennials that involved gene flow from several wild species from multiple geographical areas. This study provides an example of a complex protracted process of domestication in long-lived plant perennials, and is a starting point for apple breeding programmes.

PearMODB: a multiomics database for pear (Pyrus) genomics, genetics and breeding study

Pear (Pyrus ssp.) belongs to Rosaceae and is an important fruit tree widely cultivated around the world. Currently, challenges to cope with the burgeoning sets of multiomics data are rapidly increasing. Here, we constructed the Pear Multiomics Database (PearMODB) by integrating genome, transcriptome, epigenome and population variation data, and aimed to provide a portal for accessing and analyzing pear multiomics data. A variety of online tools were built including gene search, BLAST, JBrowse, expression heatmap, synteny analysis and primer design. The information of DNA methylation sites and single-nucleotide polymorphisms can be retrieved through the custom JBrowse, providing an opportunity to explore the genetic polymorphisms linked to phenotype variation. Moreover, different gene families involving transcription factors, transcription regulators and disease resistance (nucleotide-binding site leucine-rich repeat) were identified and compiled for quick search. In particular, biosynthetic gene clusters (BGCs) were identified in pear genomes, and specialized webpages were set up to show detailed information of BGCs, laying a foundation for studying metabolic diversity among different pear varieties. Overall, PearMODB provides an important platform for pear genomics, genetics and breeding studies. Database URL http://pearomics.njau.edu.cn.

Resequencing of Rosa rugosa accessions revealed the history of population dynamics, breed origin, and domestication pathways

BACKGROUND

Rosa rugosa is a shrub that originated in China and has economic and ecological value. However, during the development of R. rugosa, the genetic background was chaotic, and the genetic structure among different wild populations was unclear, as well as wild and cultivated accessions. Here, we report whole-genome resequencing of wild and cultivated R. rugosa accessions.

RESULTS

A total of 19,041,284 SNPs were identified in 188 R. rugosa accessions and 3 R. chinensis accessions by resequencing. Population genetic analysis revealed that cultivated and wild groups were separated very early. All R. rugosa accessions were divided into 8 categories based on genetic structure: (1) Weihai, Yantai, and Liaoning category, (2) Jilin category, and (3) Hammonasset category (above three are wild); (4) traditional varieties, (5) hybrids between R. rugosa and R. chinensis, (6) Zizhi Rose, (7) Kushui Rose, (8) hybrids between R. rugosa and R. multiflora. We found that the heterozygosity and genetic diversity of wild accessions were generally lower than those of cultivated individuals. The genes that were selected during cultivation were identified, and it was found that these genes were mainly related to environmental adaptation and growth.

CONCLUSIONS

The Jilin population was the oldest population and later migrated to Liaoning and then migrated to Yantai and Weihai by sea regression in the Bohai Basin. The Hammonasset naturalized population probably originated from the Jilin population and then experienced separate differentiation. The long-term asexual reproduction pattern of R. rugosa decreased genetic diversity in the wild population. During R. rugosa cultivation, the ancestors of the Jilin population were involved in breeding traditional varieties, after which almost no wild individuals were engaged in breeding. However, in recent decades, cross breeding of R. rugosa started the utilization of wild germplasms. In comparison, some other species play important roles in variety formation. Few genes related to economic traits were selected, suggesting no directional domestication in the R. rugosa cultivation process.

The role of invasive plant species in drought resilience in agriculture: the case of sweet briar (Rosa rubiginosa L.)

Sweet briar (Rosa rubiginosa) belongs to the group of wild roses. Under natural conditions it grows throughout Europe, and was introduced also into the southern hemisphere, where it has efficiently adapted to dry lands. This review focuses on the high adaptation potential of sweet briar to soil drought in the context of global climatic changes, especially considering steppe formation and desertification of agricultural, orchard, and horticultural areas. We provide a comprehensive overview of current knowledge on sweet briar traits associated with drought tolerance and particularly water use efficiency, sugar accumulation, accumulation of CO2 in intercellular spaces, stomatal conductance, gibberellin level, effective electron transport between photosystem II and photosystem I, and protein content. We discuss the genetics and potential applications in plant breeding and suggest future directions of study concerning invasive populations of R. rubiginosa. Finally, we point out that sweet briar can provide new genes for breeding in the context of depleting gene pools of the crop plants.

A large-scale proteogenomic atlas of pear

Pear is an important fruit tree that is widely distributed around the world. The first pear genome map was reported from our laboratory approximately 10 years ago. To further study global protein expression patterns in pear, we generated pear proteome data based on 24 major tissues. The tissue-resolved profiles provided evidence of the expression of 17 953 proteins. We identified 4294 new coding events and improved the pear genome annotation via the proteogenomic strategy based on 18 090 peptide spectra with peptide spectrum matches >1. Among the eight randomly selected new short coding open reading frames that were expressed in the style, four promoted and one inhibited the growth of pear pollen tubes. Based on gene coexpression module analysis, we explored the key genes associated with important agronomic traits, such as stone cell formation in fruits. The network regulating the synthesis of lignin, a major component of stone cells, was reconstructed, and receptor-like kinases were implicated as core factors in this regulatory network. Moreover, we constructed the online database PearEXP (http://www.peardb.org.cn) to enable access to the pear proteogenomic resources. This study provides a paradigm for in-depth proteogenomic studies of woody plants.

Advances in genomics for diversity studies and trait improvement in temperate fruit and nut crops under changing climatic scenarios

Genetic improvement of temperate fruit and nut crops through conventional breeding methods is not sufficient alone due to its extreme time-consuming, cost-intensive, and hard-to-handle approach. Again, few other constraints that are associated with these species, viz., their long juvenile period, high heterozygosity, sterility, presence of sexual incompatibility, polyploidy, etc., make their selection and improvement process more complicated. Therefore, to promote precise and accurate selection of plants based on their genotypes, supplement of advanced biotechnological tools, viz., molecular marker approaches along with traditional breeding methods, is highly required in these species. Different markers, especially the molecular ones, enable direct selection of genomic regions governing the trait of interest such as high quality, yield, and resistance to abiotic and biotic stresses instead of the trait itself, thus saving the overall time and space and helping screen fruit quality and other related desired traits at early stages. The availability of molecular markers like SNP (single-nucleotide polymorphism), DArT (Diversity Arrays Technology) markers, and dense molecular genetic maps in crop plants, including fruit and nut crops, led to a revelation of facts from genetic markers, thus assisting in precise line selection. This review highlighted several aspects of the molecular marker approach that opens up tremendous possibilities to reveal valuable information about genetic diversity and phylogeny to boost the efficacy of selection in temperate fruit crops through genome sequencing and thus cultivar improvement with respect to adaptability and biotic and abiotic stress resistance in temperate fruit and nut species.

The Role of Italy in the Use of Advanced Plant Genomic Techniques on Fruit Trees: State of the Art and Future Perspectives

Climate change is deeply impacting the food chain production, lowering quality and yield. In this context, the international scientific community has dedicated many efforts to enhancing resilience and sustainability in agriculture. Italy is among the main European producers of several fruit trees; therefore, national research centers and universities undertook several initiatives to maintain the specificity of the 'Made in Italy' label. Despite their importance, fruit crops are suffering from difficulties associated with the conventional breeding approaches, especially in terms of financial commitment, land resources availability, and long generation times. The 'new genomic techniques' (NGTs), renamed in Italy as 'technologies for assisted evolution' (TEAs), reduce the time required to obtain genetically improved cultivars while precisely targeting specific DNA sequences. This review aims to illustrate the role of the Italian scientific community in the use of NGTs, with a specific focus on Citrus, grapevine, apple, pear, chestnut, strawberry, peach, and kiwifruit. For each crop, the key genes and traits on which the scientific community is working, as well as the technological improvements and advancements on the regeneration of local varieties, are presented. Lastly, a focus is placed on the legal aspects in the European and in Italian contexts.

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.

Building a foundation for gene family analysis in Rosaceae genomes with a novel workflow: A case study in Pyrus architecture genes

The rapid development of sequencing technologies has led to a deeper understanding of plant genomes. However, direct experimental evidence connecting genes to important agronomic traits is still lacking in most non-model plants. For instance, the genetic mechanisms underlying plant architecture are poorly understood in pome fruit trees, creating a major hurdle in developing new cultivars with desirable architecture, such as dwarfing rootstocks in European pear (Pyrus communis). An efficient way to identify genetic factors for important traits in non-model organisms can be to transfer knowledge across genomes. However, major obstacles exist, including complex evolutionary histories and variable quality and content of publicly available plant genomes. As researchers aim to link genes to traits of interest, these challenges can impede the transfer of experimental evidence across plant species, namely in the curation of high-quality, high-confidence gene models in an evolutionary context. Here we present a workflow using a collection of bioinformatic tools for the curation of deeply conserved gene families of interest across plant genomes. To study gene families involved in tree architecture in European pear and other rosaceous species, we used our workflow, plus a draft genome assembly and high-quality annotation of a second P. communis cultivar, 'd'Anjou.' Our comparative gene family approach revealed significant issues with the most recent 'Bartlett' genome - primarily thousands of missing genes due to methodological bias. After correcting assembly errors on a global scale in the 'Bartlett' genome, we used our workflow for targeted improvement of our genes of interest in both P. communis genomes, thus laying the groundwork for future functional studies in pear tree architecture. Further, our global gene family classification of 15 genomes across 6 genera provides a valuable and previously unavailable resource for the Rosaceae research community. With it, orthologs and other gene family members can be easily identified across any of the classified genomes. Importantly, our workflow can be easily adopted for any other plant genomes and gene families of interest.

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.

Genome-wide identification of apple PPI genes and a functional analysis of the response of MxPPI1 to Fe deficiency stress

Iron (Fe) deficiency affects plant growth and development. The proton pump interactor (PPI) in plants responds to multiple abiotic stresses, although it has not been well characterized under Fe deficiency stress. In this study, we systematically identified and analyzed the PPI gene family in apple. Three PPI candidate genes were found, and they contained 318-1349 amino acids and 3-7 introns. Under Fe deficiency stress, we analyzed the expression of all the PPI genes in roots of apple rootstock Malus xiaojinensis. Expression of the gene MD11G1247800, designated PPI1, is obviously induced by Fe deficiency treatment in M. xiaojinensis. We first cloned MxPPI1 from M. xiaojinensis and determined its subcellular localization, which indicated that it is localized in the cell membrane and nucleus in tobacco. We found that the level of expression of the MxPPI1 protein increased significantly under Fe deficiency stress in apple calli. Moreover, overexpressing MxPPI1 in apple calli enhanced the activities of ferric chelate reductase and H⁺-ATPase, H⁺ secretion, MxHA2 gene expression and total Fe content when compared with the wild type calli. We further found that MxPPI1 interacted with MxHA2 using bimolecular fluorescence complementation and luciferase complementation assays. Overall, we demonstrated that MxPPI1 interacts with MxHA2 to enhance the activity of H⁺-ATPase to regulate Fe absorption in M. xiaojinensis.

The mechanism analysis of exogenous melatonin in limiting pear fruit aroma decrease under low temperature storage

Exogenous melatonin (MT) is widely used in fruit preservation, and can increase the storage time and delay the quality deterioration. Firstly, it was found that 150 μM MT was the optimal concentration to treat 'Xinli No.7' under storage at 4 °C for 60 days. MT could significantly improve oxidase activity and inhibit the reduction of physiological indexes, including pulp hardness, weight loss, titratable acid and soluble solid content. MT could also reduce ethylene release and limit the reduction of fruit aroma. The average content of fruit aroma substance increased by 43.53%. A relevant RNA-Seq database was built to further explore the regulation mechanism of MT. A total of 2,761 differentially expressed genes (DEGs) were identified. DEGs were enriched in 64 functional groups and 191 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. DEGs were mainly enriched in alpha-linolenic acid metabolism, fatty acid metabolism and plant hormone signal transduction pathway. The gene pycom09g05270 belonging to long chain acyl-CoA synthetase family and participating in fatty acid metabolism pathway was identified, and its expression level was consistent with fragments per kilobase per million mapped reads (FPKM) values, implying that pycom09g05270 might play a vital role in maintaining quality during the storage process.

PpZAT5 suppresses the expression of a B-box gene PpBBX18 to inhibit anthocyanin biosynthesis in the fruit peel of red pear

BBX (B-box) proteins play a vital role in light-induced anthocyanin biosynthesis. PpBBX18 was an indispensable regulator for the induction of anthocyanin biosynthesis in the peel of red pear fruit (Pyrus pyrifolia Nakai.). However, the upstream regulation of BBX genes has not been well characterized. In this study, PpZAT5, a cysteine2/histidine2-type transcription factor, was discovered as the upstream negative regulator of PpBBX18. The results showed that PpZAT5 functions as a transcriptional repressor and directly binds to the CAAT motif of PpBBX18 and inhibits its expression. PpZAT5 expression was inhibited by light, which is converse to the expression pattern of anthocyanin-related structural genes. In addition, less anthocyanin accumulated in the PpZAT5-overexpressing pear calli than in the wild-type pear calli; on the contrary, more anthocyanin accumulated in PpZAT5-RNAi pear calli. Moreover, the crucial genes involved in light-induced anthocyanin biosynthesis were markedly down-regulated in the transcriptome of PpZAT5 overexpression pear calli compared to wild-type. In conclusion, our study indicates that PpBBX18 is negatively regulated by a C2H2-type transcriptional repressor, PpZAT5, which reduces anthocyanin content in pear. The present results demonstrate an upstream molecular mechanism of PpBBX18 and provide insights into light-induced anthocyanin biosynthesis.

Next-generation sequencing in the biodiversity conservation of endangered medicinal plants

Medicinal plants have been used as traditional herbal medicines in the treatment of various types of diseases. However, the increased demand for these plants highlights the importance of conservation specifically for endangered species. Significant advancements in next-generation sequencing (NGS) technologies have accelerated medicinal plant research while reducing costs and time demands. NGS systems enable high-throughput whole genome sequencing as well as direct RNA sequencing and transcriptome analysis. The sequence data sets created can be used in a variety of areas of study, including biodiversity conservation, comparative genomics, transcriptomic analysis, single cell mining, metagenomics, epigenetics, molecular marker discovery, multi genome sequencing, and so on. Commercial sequencing service providers are constantly working to improve technologies to address bioinformatics problems in NGS data analysis. Several genome sequencing projects on medicinal plants have been completed recently and a few more are in the works. In some medicinal plants, massive NGS-based data has been developed. In the present review, we have attempted to briefly discuss advancements in NGS technology on medicinally essential plants in India. The review will also provide ideas for applying NGS technologies for exploring genomes of various endangered medicinal plants whose genome sequences are not normally available and thus provides valuable insights for the conservation of these vulnerable species.

The roles of WRKY transcription factors in Malus spp. and Pyrus spp

The WRKY transcription factor gene family is known to be involved in plant defense against pathogens and in tolerance to different environmental stresses at different stages of development. The response mechanisms through which these genes act can be influenced by different phytohormones as well as by many trans- and cis-acting elements, making this network an important topic for analysis, but still something complex to fully understand. According to available reports, these genes can also perform important roles in pome species (Malus spp. and Pyrus spp.) metabolism, especially in adaptation of these plants to stressful conditions. Here, we present a quick review of what is known about WRKY genes in Malus and Pyrus genomes offering a simple way to understand what is already known about this topic. We also add information connecting the evolution of these transcription factors with others that can also be found in pomes.

Computational analysis of potential candidate genes involved in the cold stress response of ten Rosaceae members

Background

Plant species from Rosaceae family are economically important. One of the major environmental factors impacting those species is cold stress. Although several Rosaceae plant genomes have recently been sequenced, there have been very few research conducted on cold upregulated genes and their promoter binding sites. In this study, we used computational approaches to identify and analyse potential cold stress response genes across ten Rosaceae family members.

Results

Cold stress upregulated gene data from apple and strawberry were used to identify syntelogs in other Rosaceae species. Gene duplication analysis was carried out to better understand the distribution of these syntelog genes in different Rosaceae members. A total of 11,145 popular abiotic stress transcription factor-binding sites were identified in the upstream region of these potential cold-responsive genes, which were subsequently categorised into distinct transcription factor (TF) classes. MYB classes of transcription factor binding site (TFBS) were abundant, followed by bHLH, WRKY, and AP2/ERF. TFBS patterns in the promoter regions were compared among these species and gene families, found to be quite different even amongst functionally related syntelogs. A case study on important cold stress responsive transcription factor family, AP2/ERF showed less conservation in TFBS patterns in the promoter regions. This indicates that syntelogs from the same group may be comparable at the gene level but not at the level of cis-regulatory elements. Therefore, for such genes from the same family, different repertoire of TFs could be recruited for regulation and expression. Duplication events must have played a significant role in the similarity of TFBS patterns amongst few syntelogs of closely related species.

Conclusions

Our study overall suggests that, despite being from the same gene family, different combinations of TFs may play a role in their regulation and expression. The findings of this study will provide information about potential genes involved in the cold stress response, which will aid future functional research of these gene families involved in many important biological processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08751-x.

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.

Genome-wide clonal variability in European pear "Rocha" using high-throughput sequencing

Pears (Pyrus) are one of the most economically important fruits worldwide. The Pyrus genus is characterized by a high degree of genetic variability between species and interspecific hybrids, and several studies have been performed to assess this variability for both cultivated and wild accessions. These studies have mostly been limited by the resolving power of traditional molecular markers, although in the recent past the availability of reference genome sequences or SNP arrays for pear have enhanced the capability of high-resolution genomics studies. These tools can also be applied to better understand the intra-varietal (or clonal) variability in pear. Here we report the first high resolution genomics analysis of a pear clonal population using whole genome sequencing (WGS). Results showed unique signatures for the accumulation of mutations and transposable element insertions in each clone, which are likely related to their history of propagation and cultivation. The nucleotide diversity remained low in the clonal collection with the exception of few genomic windows, suggesting that balancing selection may be occurring. These windows included mainly genes related to plant fertility. Regions with higher mutational load were partially associated with transcription factors, probably reflecting the distinctive phenotypes in the collection. The annotation of variants also revealed the theoretical disruption of relevant genes in pear. Taken together, the results from this study show that pear clones accumulate mutations differently, and that those mutations can play a role on pear phenotypes, meaning that the study of pear clonal populations can be relevant in genetic studies, mainly when comparing with traditional association studies.

Overexpression of apple Ma12, a mitochondrial pyrophosphatase pump gene, leads to malic acid accumulation and the upregulation of malate dehydrogenase in tomato and apple calli

Acidity is an important factor influencing the organoleptic quality of apple fruits. In this study, an apple pyrophosphate-energized proton pump (PEPP) gene was isolated and designated MdMa12. On the basis of a phylogenetic analysis in Rosaceae species, PEPP genes were divided into three groups, with apple PEPP genes most closely related to pear PEPP genes. Gene expression analysis revealed that high malic acid content was generally accompanied by high MdMa12 expression levels. Moreover, MdMa12 was mainly expressed in the fruit. A subcellular localization analysis suggested that MdMa12 is a mitochondrial protein. The ectopic expression and overexpression of MdMa12 in "Micro-Tom" tomato and apple calli, respectively, increased the malic acid content. One (MDH12) of four malate dehydrogenase genes highly expressed in transgenic apple calli was confirmed to encode a protein localized in mitochondria. The overexpression of MDH12 increased the malate content in apple calli. Furthermore, MdMa12 overexpression increased MdDTC1, MdMa1, and MdMa10 expression levels, which were identified to transport malate. These findings imply that MdMa12 has important functions related to apple fruit acidity. Our study explored the regulatory effects of mitochondria on the complex mechanism underlying apple fruit acidity.

Pear genetics: Recent advances, new prospects, and a roadmap for the future

Pear, belonging to the genus Pyrus, is one of the most economically important temperate fruit crops. Pyrus is an important genus of the Rosaceae family, subfamily Maloideae, and has at least 22 different species with over 5000 accessions maintained or identified worldwide. With the release of draft whole-genome sequences for Pyrus, opportunities for pursuing studies on the evolution, domestication, and molecular breeding of pear, as well as for conducting comparative genomics analyses within the Rosaceae family, have been greatly expanded. In this review, we highlight key advances in pear genetics, genomics, and breeding driven by the availability of whole-genome sequences, including whole-genome resequencing efforts, pear domestication, and evolution. We cover updates on new resources for undertaking gene identification and molecular breeding, as well as for pursuing functional validation of genes associated with desirable economic traits. We also explore future directions for "pear-omics".

PlantTribes2: Tools for comparative gene family analysis in plant genomics

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

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.

Penicillium and Talaromyces Species as Postharvest Pathogens of Pear Fruit (Pyrus communis) in Serbia

Pears are one of the oldest and the third most important fruit species grown in temperate regions. They are consumed because of their nutritional and health benefits, in fresh form or as various processed products. This article resolves the etiology of the Penicillium-like mold symptoms on pear fruits in Serbia. Samples of pear fruits with blue mold and other Penicillium-like mold symptoms were collected in Serbia from 2016 to 2019, from four storages. The recovered isolates were identified and characterized according to a polyphasic approach. Morphological and physiological analyses were performed on three media and five temperatures, respectively. Four loci (internal transcribed spacer, beta-tubulin, calmodulin, and DNA-dependent RNA polymerase II second largest subunit) were used for sequencing, genetic identification, and phylogenetic analyses. The results of the identification by conventional and molecular methods were in agreement, and they revealed that the obtained isolates belong to five species: Penicillium crustosum, P. expansum, P. italicum, Talaromyces minioluteus, and T. rugulosus. In a pathogenicity test, P. crustosum, P. expansum, T. minioluteus, and T. rugulosus produced decay on artificially inoculated pear fruits, and P. italicum induced tissue response lesions. The results of this study are the first reports of T. minioluteus and T. rugulosus as postharvest pear pathogens. Also, these are the first world records of T. minioluteus, T. rugulosus, and P. italicum on fruits of European pear. Furthermore, this is the first finding of P. crustosum, P. expansum, P. italicum, T. minioluteus, and T. rugulosus on pear fruit in Serbia.

Genome-Wide Identification, Characterization and Expression Profiling of Aluminum-Activated Malate Transporters in Eriobotrya japonica Lindl

Aluminum-activated malate transporters (ALMTs) have multiple potential roles in plant metabolism such as regulation of organic acids in fruits, movement of guard cells and inducing tolerance against aluminum stress. However, the systematic characterization of ALMT genes in loquat is yet to be performed. In the current study, 24 putative ALMT genes were identified in the genome of Eriobotrya japonica Lindl. To further investigate the role of those ALMT genes, comprehensive bioinformatics and expression analysis were performed. In bioinformatics analysis, the physiochemical properties, conserved domains, gene structure, conserved motif, phylogenetic and syntenic analysis of EjALMT genes were conducted. The result revealed that the ALMT superfamily domain was conserved in all EjALMT proteins. EjALMT proteins were predicted to be localized in the plasma membrane. Genomic structural and motif analysis showed that the exon and motif number of each EjALMT gene ranged dramatically, from 5 to 7, and 6 to 10, respectively. Syntenic analysis indicated that the segmental or whole-genome duplication played a vital role in extension of the EjALMT gene family. The Ka and Ks values of duplicated genes depicted that EjALMT genes have undergone a strong purifying selection. Furthermore, the expression analysis of EjALMT genes was performed in the root, mature leaf, stem, full-bloom flower and ripened fruit of loquat. Some genes were expressed differentially in examined loquat tissues, signifying their differential role in plant growth and development. This study provides the first genome-wide identification, characterization, and relative expression of the ALMT gene family in loquat and provides the foundation for further functional analysis.

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

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

New high-quality peach (Prunus persica L. Batsch) genome assembly to analyze the molecular evolutionary mechanism of volatile compounds in peach fruits

Peach (Prunus persica L. Batsch) is an economically important fruit crop worldwide. Although a high-quality peach genome has previously been published, Sanger sequencing was used for its assembly, which generated short contigs. Here, we report a chromosome-level genome assembly and sequence analysis of Chinese Cling, an important founder cultivar for peach breeding programs worldwide. The assembled genome contained 247.33 Mb with a contig N50 of 4.13 Mb and a scaffold N50 of 29.68 Mb, representing 99.8% of the estimated genome. Comparisons between this genome and the recently published one (Lovell peach) uncovered 685 407 single nucleotide polymorphisms, 162 655 insertions and deletions, and 16 248 structural variants. Gene family analysis highlighted the contraction of the gene families involved in flavone, flavonol, flavonoid, and monoterpenoid biosynthesis. Subsequently, the volatile compounds of 256 peach varieties were quantitated in mature fruits in 2015 and 2016 to perform a genome-wide association analysis. A comparison with the identified domestication genomic regions allowed us to identify 25 quantitative trait loci, associated with seven volatile compounds, in the domestication region, which is consistent with the differences in volatile compounds between wild and cultivated peaches. Finally, a gene encoding terpene synthase, located within a previously reported quantitative trait loci region, was identified to be associated with linalool synthesis. Such findings highlight the importance of this new assembly for the analysis of evolutionary mechanisms and gene identification in peach species. Furthermore, this high-quality peach genome provides valuable information for future fruit improvement.

Development of Genome-Wide Functional Markers Using Draft Genome Assembly of Guava (Psidium guajava L.) cv. Allahabad Safeda to Expedite Molecular Breeding

Guava (Psidium guajava L.), a rich source of nutrients, is an important tropical and subtropical fruit of the Myrtaceae family and exhibits magnificent diversity. Genetic diversity analysis is the first step toward the identification of parents for hybridization, genetic mapping, and molecular breeding in any crop species. A diversity analysis based on whole-genome functional markers increases the chances of identifying genetic associations with agronomically important traits. Therefore, here, we sequenced the genome of guava cv. Allahabad Safeda on an Illumina platform and generated a draft assembly of ~304 MB. The assembly of the Allahabad Safeda genome constituted >37.95% repeat sequences, gene prediction with RNA-seq data as evidence identified 14,115 genes, and BLAST n/r, Interproscan, PfamScan, BLAST2GO, and KEGG annotated 13,957 genes. A comparative protein transcript analysis of tree species revealed the close relatedness of guava with Eucalyptus. Comparative transcriptomics-based SSR/InDel/SNP-PCR ready genome-wide markers in greenish-yellow skinned and white fleshed-Allahabad Safeda to four contrasting cultivars viz apple-color-skinned and white-fleshed-Lalima, greenish-yellow-skinned and pink-fleshed-Punjab Pink, purple-black-skinned and purple-fleshed-Purple Local and widely used rootstock-Lucknow-49 were developed. The molecular markers developed here revealed a high level of individual heterozygosity within genotypes in 22 phenotypically diverse guava cultivars. Principal coordinate, STRUCTURE clustering, and neighbor-joining-based genetic diversity analysis identified distinct clusters associated with fruit skin and flesh color. The genome sequencing of guava, functional annotation, comparative transcriptomics-based genome-wide markers, and genetic diversity analysis will expand the knowledge of genomes of climacteric fruits, facilitating trait-based molecular breeding and diversifying the nutritional basket.

Genome-wide association study of individual sugar content in fruit of Japanese pear (Pyrus spp.)

BACKGROUND

Understanding mechanisms of sugar accumulation and composition is essential to determining fruit quality and maintaining a desirable balance of sugars in plant storage organs. The major sugars in mature Rosaceae fruits are sucrose, fructose, glucose, and sorbitol. Among these, sucrose and fructose have high sweetness, whereas glucose and sorbitol have low sweetness. Japanese pear has extensive variation in individual sugar contents in mature fruit. Increasing total sugar content and that of individual high-sweetness sugars is a major target of breeding programs. The objective of this study was to identify quantitative trait loci (QTLs) associated with fruit traits including individual sugar accumulation, to infer the candidate genes underlying the QTLs, and to assess the potential of genomic selection for breeding pear fruit traits.

RESULTS

We evaluated 10 fruit traits and conducted genome-wide association studies (GWAS) for 106 cultivars and 17 breeding populations (1112 F1 individuals) using 3484 tag single-nucleotide polymorphisms (SNPs). By implementing a mixed linear model and a Bayesian multiple-QTL model in GWAS, 56 SNPs associated with fruit traits were identified. In particular, a SNP located close to acid invertase gene PPAIV3 on chromosome 7 and a newly identified SNP on chromosome 11 had quite large effects on accumulation of sucrose and glucose, respectively. We used 'Golden Delicious' doubled haploid 13 (GDDH13), an apple reference genome, to infer the candidate genes for the identified SNPs. In the region flanking the SNP on chromosome 11, there is a tandem repeat of early responsive to dehydration (ERD6)-like sugar transporter genes that might play a role in the phenotypes observed.

CONCLUSIONS

SNPs associated with individual sugar accumulation were newly identified at several loci, and candidate genes underlying QTLs were inferred using advanced apple genome information. The candidate genes for the QTLs are conserved across Pyrinae genomes, which will be useful for further fruit quality studies in Rosaceae. The accuracies of genomic selection for sucrose, fructose, and glucose with genomic best linear unbiased prediction (GBLUP) were relatively high (0.67-0.75), suggesting that it would be possible to select individuals having high-sweetness fruit with high sucrose and fructose contents and low glucose content.

Transcriptome Characterization and Expression Profiles of Disease Defense-Related Genes of Table Grapes in Response to Pichia anomala Induced with Chitosan

Transcriptome analysis (TA) was conducted to characterize the transcriptome changes in postharvest disease-related genes of table grapes following treatment with Pichia anomala induced with chitosan (1% w/v). In the current study, the difference in the gene expression of table grapes after treatment with P. anomala induced with chitosan and that of a control group was compared 72 h post-inoculation. The study revealed that postharvest treatment of table grapes with P. anomala induced with chitosan could up-regulate genes that have a pivotal role in the fruit’s disease defense. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) results also confirmed that GO terms and the KEGG pathways, which have pivotal roles in plant disease resistance, were significantly enriched. The up-regulated genes of the treatment group have a unique function in the fruit’s disease resistance compared to the control group. Generally, most genes in the plant–pathogen interaction pathway; the plant Mitogen-activated protein kinase (MAPK) signaling pathway; the plant hormone signal transduction pathway; the pathway of glutathione metabolism; the pathway of phenylalanine, tyrosine, and tryptophan biosynthesis; and the pathway of flavonoid biosynthesis were all up-regulated. These up-regulations help the fruit to synthesize disease-resistant substances, regulate the reactive oxygen species (ROS), enhance the fruit cell wall, and enrich hormone signal transduction during the pathogen’s attack. This study is useful to overcome the lags in applying transcriptomics technology in postharvest pathology, and will provide insight towards developing other alternative methods to using bio-pesticides to control postharvest diseases of perishables.

Transcriptomics of Differential Ripening in 'd'Anjou' Pear (Pyrus communis L.)

Estimating maturity in pome fruits is a critical task that directs virtually all postharvest supply chain decisions. This is especially important for European pear (Pyrus communis) cultivars because losses due to spoilage and senescence must be minimized while ensuring proper ripening capacity is achieved (in part by satisfying a fruit chilling requirement). Reliable methods are lacking for accurate estimation of pear fruit maturity, and because ripening is maturity dependent it makes predicting ripening capacity a challenge. In this study of the European pear cultivar 'd'Anjou', we sorted fruit at harvest based upon on-tree fruit position to build contrasts of maturity. Our sorting scheme showed clear contrasts of maturity between canopy positions, yet there was substantial overlap in the distribution of values for the index of absorbance difference (I _AD ), a non-destructive spectroscopic measurement that has been used as a proxy for pome fruit maturity. This presented an opportunity to explore a contrast of maturity that was more subtle than I _AD could differentiate, and thus guided our subsequent transcriptome analysis of tissue samples taken at harvest and during storage. Using a novel approach that tests for condition-specific differences of co-expressed genes, we discovered genes with a phased character that mirrored our sorting scheme. The expression patterns of these genes are associated with fruit quality and ripening differences across the experiment. Functional profiles of these co-expressed genes are concordant with previous findings, and also offer new clues, and thus hypotheses, about genes involved in pear fruit quality, maturity, and ripening. This work may lead to new tools for enhanced postharvest management based on activity of gene co-expression modules, rather than individual genes. Further, our results indicate that modules may have utility within specific windows of time during postharvest management of 'd'Anjou' pear.

Study on the differences of gene expression between pear and apple wild cultivation materials based on RNA-seq technique

BACKGROUND

Pears and apples are both perennial deciduous trees of the Rosaceae family, and both are important economic fruit trees worldwide. The emergence of many varieties in the market has been mostly domesticated from wild to cultivated and regulated by the differential expression of genes. However, the molecular process and pathways underlying this phenomenon remain unclear. Four typical wild and cultivar pear and apple trees at three developmental stages were used in our study to investigate the molecular process at the transcriptome level.

RESULT

Physiological observations indicated the obvious differences of size, weight, sugar acid content and peel color in wild and cultivar fruit among each developmental stage. Using next-generation sequencing based RNA-seq expression profiling technology, we produced a transcriptome in procession of a large fraction of annotated pear and apple genes, and provided a molecular basis underlying the phenomenon of wild and cultivar fruit tree differences. 5921 and 5744 differential expression genes were identified in pear and apple at three developmental stages respectively. We performed temporal and spatial differential gene expression profiling in developing fruits. Several key pathways such as signal transduction, photosynthesis, translation and many metabolisms were identified as involved in the differentiation of wild and cultivar fruits.

CONCLUSION

In this study, we reported on the next-generation sequencing study of the temporal and spatial mRNA expression profiling of pear and apple fruit trees. Also, we demonstrated that the integrated analysis of pear and apple transcriptome, which strongly revealed the consistent process of domestication in Rosaceae fruit trees. The results will be great influence to the improvement of cultivar species and the utilization of wild resources.

Characterization of the pectin methylesterase inhibitor gene family in Rosaceae and role of PbrPMEI23/39/41 in methylesterified pectin distribution in pear pollen tube

Pectin methylesterase inhibitor gene family in the seven Rosaceae species (including three pear cultivars) is characterized and three pectin methylesterase inhibitor genes are identified to regulate pollen tube growth in pear. Pectin methylesterase inhibitor (PMEI) participates in a variety of biological processes in plants. However, the information and function of PMEI genes in Rosaceae are largely unknown. In this study, a total of 423 PMEI genes are identified in the genomes of seven Rosaceae species. The PMEI genes in pear are categorized into five subfamilies based on structural analysis and evolutionary analysis. WGD and TD are the main duplication events in the PMEI gene family of pear. Quantitative real-time PCR analysis indicates that PbrPMEI23, PbrPMEI39, and PbrPMEI41 are increasingly expressed during pear pollen tube growth. Under the treatment of recombinant proteins PbrPMEI23, PbrPMEI39 or PbrPMEI41, the content of methylesterified pectin at the region 5-20 μm from the pollen tube tip significantly increases, and the growth of pear pollen tubes is promoted. These results indicate that PMEI regulates the growth of pollen tubes by changing the distribution of methylesterified pectin in the apex.

An Atlas of Genomic Resources for Studying Rosaceae Fruits and Ornamentals

Rosaceae, a large plant family of more than 3,000 species, consists of many economically important fruit and ornamental crops, including peach, apple, strawberry, raspberry, cherry, and rose. These horticultural crops are not only important economic drivers in many regions of the world, but also major sources of human nutrition. Additionally, due to the diversity of fruit types in Rosaceae, this plant family offers excellent opportunities for investigations into fleshy fruit diversity, evolution, and development. With the development of high-throughput sequencing technologies and computational tools, an increasing number of high-quality genomes and transcriptomes of Rosaceae species have become available and will greatly facilitate Rosaceae research and breeding. This review summarizes major genomic resources and genome research progress in Rosaceae, highlights important databases, and suggests areas for further improvement. The availability of these big data resources will greatly accelerate research progress and enhance the agricultural productivity of Rosaceae.

Genome-wide characterization of PEBP family genes in nine Rosaceae tree species and their expression analysis in P. mume

BACKGROUND

Phosphatidylethanolamine-binding proteins (PEBPs) constitute a common gene family found among animals, plants and microbes. Plant PEBP proteins play an important role in regulating flowering time, plant architecture as well as seed dormancy. Though PEBP family genes have been well studied in Arabidopsis and other model species, less is known about these genes in perennial trees.

RESULTS

To understand the evolution of PEBP genes and their functional roles in flowering control, we identified 56 PEBP members belonging to three gene clades (MFT-like, FT-like, and TFL1-like) and five lineages (FT, BFT, CEN, TFL1, and MFT) across nine Rosaceae perennial species. Structural analysis revealed highly conserved gene structure and protein motifs among Rosaceae PEBP proteins. Codon usage analysis showed slightly biased codon usage across five gene lineages. With selection pressure analysis, we detected strong purifying selection constraining divergence within most lineages, while positive selection driving the divergence of FT-like and TFL1-like genes from the MFT-like gene clade. Spatial and temporal expression analyses revealed the essential role of FT in regulating floral bud breaking and blooming in P. mume. By employing a weighted gene co-expression network approach, we inferred a putative FT regulatory module required for dormancy release and blooming in P. mume.

CONCLUSIONS

We have characterized the PEBP family genes in nine Rosaceae species and examined their phylogeny, genomic syntenic relationship, duplication pattern, and expression profiles during flowering process. These results revealed the evolutionary history of PEBP genes and their functions in regulating floral bud development and blooming among Rosaceae tree species.

De novo assembly and characterization of the first draft genome of quince (Cydonia oblonga Mill.)

Quince (Cydonia oblonga Mill.) is the sole member of the genus Cydonia in the Rosacea family and closely related to the major pome fruits, apple (Malus domestica Borkh.) and pear (Pyrus communis L.). In the present work, whole genome shotgun paired-end sequencing was employed in order to assemble the first draft genome of quince. A genome assembly that spans 488.4 Mb of sequence corresponding to 71.2% of the estimated genome size (686 Mb) was produced in the study. Gene predictions via ab initio and homology-based sequence annotation strategies resulted in the identification of 25,428 and 30,684 unique putative protein coding genes, respectively. 97.4 and 95.6% of putative homologs of Arabidopsis and rice transcription factors were identified in the ab initio predicted genic sequences. Different machine learning algorithms were tested for classifying pre-miRNA (precursor microRNA) coding sequences, identifying Support Vector Machine (SVM) as the best performing classifier. SVM classification predicted 600 putative pre-miRNA coding loci. Repetitive DNA content of the assembly was also characterized. The first draft assembly of the quince genome produced in this work would constitute a foundation for functional genomic research in quince toward dissecting the genetic basis of important traits and performing genomics-assisted breeding.

Adaptive evolution driving the young duplications in six Rosaceae species

Background

In plant genomes, high proportions of duplicate copies reveals that gene duplications play an important role in the evolutionary processes of plant species. A series of gene families under positive selection after recent duplication events in plant genomes indicated the evolution of duplicates driven by adaptive evolution. However, the genome-wide evolutionary features of young duplicate genes among closely related species are rarely reported.

Results

In this study, we conducted a systematic survey of young duplicate genes at genome-wide levels among six Rosaceae species, whose whole-genome sequencing data were successively released in recent years. A total of 35,936 gene families were detected among the six species, in which 60.25% were generated by young duplications. The 21,650 young duplicate gene families could be divided into two expansion types based on their duplication patterns, species-specific and lineage-specific expansions. Our results showed the species-specific expansions advantaging over the lineage-specific expansions. In the two types of expansions, high-frequency duplicate domains exhibited functional preference in response to environmental stresses.

Conclusions

The functional preference of the young duplicate genes in both the expansion types showed that they were inclined to respond to abiotic or biotic stimuli. Moreover, young duplicate genes under positive selection in both species-specific and lineage-specific expansions suggested that they were generated to adapt to the environmental factors in Rosaceae species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07422-7.

New Strategies to Overcome Present CRISPR/Cas9 Limitations in Apple and Pear: Efficient Dechimerization and Base Editing

Despite recent progress, the application of CRISPR/Cas9 in perennial plants still has many obstacles to overcome. Our previous results with CRISPR/Cas9 in apple and pear indicated the frequent production of phenotypic and genotypic chimeras, after editing of the phytoene desaturase (PDS) gene conferring albino phenotype. Therefore, our first objective was to determine if adding an adventitious regeneration step from leaves of the primary transgenic plants (T0) would allow a reduction in chimerism. Among hundreds of adventitious buds regenerated from a variegated T0 line, 89% were homogeneous albino. Furthermore, the analysis of the target zone sequences of twelve of these regenerated lines (RT0 for “regenerated T0” lines) indicated that 99% of the RT0 alleles were predicted to produce a truncated target protein and that 67% of RT0 plants had less heterogeneous editing profiles than the T0. Base editors are CRISPR/Cas9-derived new genome-editing tools that allow precise nucleotide substitutions without double-stranded breaks. Hence, our second goal was to demonstrate the feasibility of CRISPR/Cas9 base editing in apple and pear using two easily scorable genes: acetolactate synthase—ALS (conferring resistance to chlorsulfuron) and PDS. The two guide RNAs under MdU3 and MdU6 promoters were coupled into a cytidine base editor harboring a cytidine deaminase fused to a nickase Cas9. Using this vector; we induced C-to-T DNA substitutions in the target genes; leading to discrete variation in the amino-acid sequence and generating new alleles. By co-editing ALS and PDS genes; we successfully obtained chlorsulfuron resistant and albino lines in pear. Overall; our work indicates that a regeneration step can efficiently reduce the initial chimerism and could be coupled with the application of base editing to create accurate genome edits in perennial plants.

Chromosome-level draft genome of a diploid plum (Prunus salicina)

BACKGROUND

Plums are one of the most economically important Rosaceae fruit crops and comprise dozens of species distributed across the world. Until now, only limited genomic information has been available for the genetic studies and breeding programs of plums. Prunus salicina, an important diploid plum species, plays a predominant role in modern commercial plum production. Here we selected P. salicina for whole-genome sequencing and present a chromosome-level genome assembly through the combination of Pacific Biosciences sequencing, Illumina sequencing, and Hi-C technology.

FINDINGS

The assembly had a total size of 284.2 Mb, with contig N50 of 1.78 Mb and scaffold N50 of 32.32 Mb. A total of 96.56% of the assembled sequences were anchored onto 8 pseudochromosomes, and 24,448 protein-coding genes were identified. Phylogenetic analysis showed that P. salicina had a close relationship with Prunus mume and Prunus armeniaca, with P. salicina diverging from their common ancestor ∼9.05 million years ago. During P. salicina evolution 146 gene families were expanded, and some cell wall-related GO terms were significantly enriched. It was noteworthy that members of the DUF579 family, a new class involved in xylan biosynthesis, were significantly expanded in P. salicina, which provided new insight into the xylan metabolism in plums.

CONCLUSIONS

We constructed the first high-quality chromosome-level plum genome using Pacific Biosciences, Illumina, and Hi-C technologies. This work provides a valuable resource for facilitating plum breeding programs and studying the genetic diversity mechanisms of plums and Prunus species.

Systematic analysis of the OFP genes in six Rosaceae genomes and their roles in stress response in Chinese pear (Pyrus bretschneideri)

OVATE family proteins (OFPs) are the plant-specific transcription factors, and have significant functions in regulating plant growth, development and resistance. The OFP genes have been investigated in several plants, but they still lack a systematic analysis of OFP genes in Chinese pear and some other five Rosaceae genomes. Here, 28 PbrOFPs were identified within Chinese pear and compared them with those of other five Rosaceae genomes. Evolutionary tree revealed that all OFP genes from six Rosaceae genomes were divided into eight groups. Seventeen conserved microsynteny regions were detected in Chinese pear genome, suggested that these PbrOFP genes might be considered to have originated from the large-scale duplication events., indicating these PbrOFP genes might contain specialized regulatory mechanisms in these tissues, such as flower, ovary and fruit. Remarkably, two PbrOFP genes (Pbr010426.1 and Pbr010427.1) were up-regulated under Venturia nashicola treatment, and five PbrOFP genes were up-regulated under PEG treatment, suggesting that these genes might play crucial roles in defence to environmental stresses. Our data presented a systematic analysis and might aid in the selection of appropriate PbrOFPs for further functional studies in Chinese pear, especially in response to the mechanism of biotic and abiotic stresses.

Reconstruction of the Largest Pedigree Network for Pear Cultivars and Evaluation of the Genetic Diversity of the USDA-ARS National Pyrus Collection

The USDA-ARS National Clonal Germplasm Repository (NCGR) in Corvallis, Oregon, maintains one of the world's largest and most diverse living Pyrus collection. A thorough genetic characterization of this germplasm will provide relevant information to optimize the conservation strategy of pear biodiversity, support the use of this germplasm in breeding, and increase our knowledge of Pyrus taxonomy, evolution, and domestication. In the last two decades simple sequence repeat (SSR) markers have been used at the NCGR for cultivar identification and small population structure analysis. However, the recent development of the Applied Biosystems Axiom Pear 70K Genotyping Array has allowed high-density single nucleotide polymorphism (SNP)-based genotyping of almost the entire collection. In this study, we have analyzed this rich dataset to discover new synonyms and mutants, identify putative labeling errors in the collection, reconstruct the largest pear cultivar pedigree and further elucidate the genetic diversity of Pyrus.